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ABSTRACT 


Current  procedures  to  assess  the  strength  and 
stiffness  of  a  steel  deck  diaphragm  under  in-plane  loads 
are  reviewed.  A  new  method  to  calculate  the  stiffness  of 
a  steel  deck  diaphragm  is  developed. 

The  influence  of  different  parameters  that  affect  the 
behaviour  of  a  steel  deck  is  analyzed. 

The  behaviour  of  one  storey  rectangular  buildings 
under  horizontal  loads  is  studied.  A  particular  solution 
based  on  the  beam  model  is  applied  for  buildings  having  a 
limited  number  of  bays.  A  general  solution  based  on  a 
beam  on  an  elastic  foundation  is  used  for  any  number  of 
bays.  This  procedure  is  applied  for  symmetrical  as  well 
as  nonsymmetr ical  buildings. 

Practical  examples  of  analysis  for  lateral  load 
resistance  of  buildings  using  the  method  developed  are 
presented . 
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CHAPTER  I 


INTRODUCTION 

1.  GENERAL  REMARKS 

When  a  structure  is  subjected  to  lateral  loads ,  a 
complex  distribution  of  forces  results  within  the  structure.  The 
first  question  that  arises  is  how  the  lateral  forces  are 
transferred  to  the  vertical  elements.  This  distribution  of 
forces  can  become  complicated  depending  on  the  geometry  of  the 
structure  and  the  stiffness  of  the  horizontal  diaphragm  and  the 
vertical  resisting  elements. 

1 . 1  Diaphragm  Action 

To  show  the  influence  of  diaphragm  action  let  us 
consider  the  one-story  building  shown  in  Figure  1.1.  As  shown, 
the  wind  excerts  a  positive  pressure  on  the  windward  side  of  the 
building.  The  wall  cladding  ADKI  is  assumed  to  have  enough 
strength  to  work  as  a  loaded  girder  spanning  vertically  between 
the  wall  foundation  and  the  perimeter  roof  beams  AB,  BC  and  CD. 
The  roof  surface  ADEH  is  loaded  in  its  own  plane  along  its  edge 
AD.  The  roof  system  transmits  this  horizontal  edge  load  from  its 
line  of  application  to  the  frames  and  shear  walls  and  hence  to 
the  building  foundations. 

In  the  case  of  an  earthquake  loading  inertia  forces  are 
important.  These  forces  depend  on  the  mass  of  the  element 
involved.  The  equivalent  static  force  which  must  be  applied  to  a 
structure  in  design  is  generally  determined  by  building  codes  as 
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Figure  M  Horizontal  Forces  on  a  Building 
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a  percentage  of  the  dead  weight  of  the  structure..  This  force  is 
assumed  to  act  at  the  roof  line  of  a.  one-storey  structure  or  at 
each  floor  line  in  the  case  of  multiple  stories.  The  roof  or 
floor  is  designed  as  a  diaphragm  to  transmit  the  force  to  the 
resisting  cross  walls. 

In  conventional  earthquake  and  wind  calculations, 
reinforced  concrete  horizontal  diaphragms  are  assumed  to  be 
infinitely  rigid  relative  to  their  vertical  supports.  The 
lateral  forces  applied  to  the  diaphragm  are  assumed  to  be 
distributed  to  the  vertical  supports  in  proportion  to  the 
relative  stiffness  of  such  supports.  However,  in  a  structure 
with  widely  spaced  shear  walls  or  in  a  structure  where  steel 
decks  or  concrete  filled  steel  decks  are  used  as  horizontal 
diaphragms,  the  deformation  of  the  diaphragm  can  be 
significant.  In  these  cases,  the  distribution  of  seismic  forces 
in  the  vertical-resistant  elements  is  markedly  different  from  the 
distribution  that  would  result  if  the  horizontal  diaphragms  were 
assumed  to  be  rigid.  If  the  diaphragms  are  very  flexible  the 
vertical  supports  are  assumed  to  be  rigid  relative  to  the 
diaphragm  and  lateral  forces  are  assumed  to  be  distributed  to  the 
vertical  supports  in  proportion  to  the  contributing  diaphragms 
areas.  For  diaphragms  intermediate  between  rigid  and  flexible, 
one  or  the  other  of  the  above  methods  is  usually  applied, 
depending,  on  which  idealized  situation  is  more  closely 
approximated.  Therefore  there  is  a  need  to  know  the  distribution 
of  forces  to  the  vertical  resisting  elements  when  the  diaphragms 
are  neither  infinitely  rigid  nor  infinitely  flexible. 
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1.2  OBJECTIVE  AND  SCOPE 

The  purpose  of  this  study  is  to  analyze  the  influence  of 
the  diaphragm  stiffness  on  the  distribution  of  lateral  forces  to 
the  structural  resisting  elements  in  buildings.  The  study  is 
limited  to  one  storey  symmetrical  and  nonsymmetr ical  rectangular 
buildings . 

Chapter  II  reviews  and  analyzes  different  methods  to 
calculate  strength  and  stiffness  for  steel  decks.  Then  a  method 
is  developed  by  the  writer  to  assess  the  stiffness  of  steel  deck 
diaphragms.  Finally  the  parameters  that  influence  the  stiffness 
and  strength  of  steel  deck  diaphragms  are  studied. 

Chapter  III  is  devoted  to  a  study  of  the  influence  of 
diaphragm  stiffness  on  the  distribution  of  lateral  forces  to  the 
vertical  resisting  elements  for  symmetrical  buildings. 

Chapter  IV  presents  the  analysis  of  the  influence  of 
diaphragm  stiffness  on  lateral  force  distribution  for 
nonsymmetr ical  buildings. 

Chapter  V  is  devoted  to  the  application  of  the  theory 
developed  in  Chapter  III  to  practical  examples. 


CHAPTER  II 


CRITICAL  REVIEW  OF  CURRENT  METHODS 
TO  EVALUATE  STRENGTH  AND  STIFFNESS 
OF  STEEL  DECK  DIAPHRAGMS 

2.1  INTRODUCTION 

The  intent  of  this  chapter  is  to  analyze  and  compare 
some  procedures,  developed  in  Great  Britain  and  North  America, 
which  are  currently  used  in  engineering  practice,  to  calculate 
the  strength  and  stiffness  of  steel  deck  diaphragms.  A  new 
approach  to  assess  the  stiffness  of  a  steel  deck  diaphragm  is 
developed  and  on  the  basis  of  tests  performed  in  England  is 
compared  to  the  above  mentioned  procedures.  At  the  end  of  the 
chapter,  the  parameters  that  influence  the  stiffness  and  strength 
of  the  steel  diaphragms  are  studied,  so  that,  when  designing  a 
diaphragm  a  swift  decision  may  be  taken. 

2.2  EVALUATION  OF  STRENGTH  AND  STIFFNESS  OF  STEEL  DECK 

DIAPHRAGMS 

The  strength  and  stiffness  of  decks  can  be  evaluated 
using  several  methods  as  detailed  below. 

Testing .  This  is  the  most  accurate  way  of  assessing  the 
strength  and  stiffness  of  diaphragms,  but  at  the  same  time,  the 
most  expensive  and  time  consuming.  Thus,  it  has  not  been  used 
extensively  in  current  engineering  practice. 

Finite  Element  Method.  The  diaphragm  is  visualized  as 
an  assemblage  of  structural  elements  interconnected  at  a  discrete 
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number  of  nodal  points.  The  corrugated  steel  sheeting  is 
represented  by  an  orthotropic  plate  element,  purlins  and  edge 
beams  are  defined  as  beam  elements  and  seam  fasteners,  end 
fasteners  and  shear  connectors  are  represented  by  spring 
elements.  This  method  yields  results  that  are  accurate  and 
additionally  provides  information  regarding  internal  force 
distributions.  Furthermore,  it  may  be  applied  to  diaphragms  with 
openings,  which  is  not  possible  currently  for  other  methods  with 
the  exception  of  testing.  The  finite  element  method  has  the 
following  shortcomings: 

(  i)  a  computer  is  needed  which  is  not  always  available  to 
designers ; 

(ii)  The  preparation  of  the  input  data  is  time  consuming; 
and , 

(iii)  numerical  errors  can  occur  because  of  the  large  amount 
of  input  data  required. 

Analytical  Method.  This  approach  uses  either  the  Energy 
or  Equilibrium  methods  or  both.  Although  these  methods  are  not 
exact  they  provide  readily  applicable  solutions  to  practically 
any  type  of  steel  diaphragm. 

Empirical  Method.  The  Empirical  method  is  based  on  the 
result  of  testing,  so  that,  it  may  be  applied  to  decks  having 
similar  characteristics.  The  results  are  given  in  tables  usually 
supplied  by  the  firms  that  manufacture  deck  material.  Empirical 
equations  have  been  provided  by  the  U.S.  Department  of  Army,  Navy 
and  Air  Force  (31)  and  by  Luttrell  (21)  whose  research  in  West 
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Virginia  University  has  made  available  information  regarding  the 
behaviour  of  steel  deck  diaphragms.  Some  of  the  companies  that 
make  use  of  these  equations  are  Westeel  Rosco,  Bowman-Peterson, 
Westland  Metals,  etc. 

Generally  what  is  required  for  design  purposes  is  a 
readily  available  design  procedure  which  provides  a  reliable 
assessment  of  the  strength  and  stiffness  of  diaphragms.  There 
are  only  two,  of  the  above  mentioned  procedures  which  satisfy 
that  requirement,  i.e.  the  Analytical  and  the  Empirical 
methods.  For  steel  deck  diaphragms  both  methods  are  available, 
for  concrete  deck  diaphragms  only  Empirical  Equations  have  been 
developed.  The  Analytical  methods  are  discussed  in  detail  in  the 
following  sections. 

2.3  BEHAVIOUR  OF  A  STEEL  DECK  DIAPHRAGM  UNDER  IN-PLANE  LOADS 

A  typical  isolated  diaphragm  is  shown  in  Figure  2.1.  It 
is  made  up  of  a  number  of  component  elements,  that  is,  corrugated 
steel  sheeting,  purlins,  rafters,  shear  connectors,  end  fasteners 
and  seam  fasteners.  Different  alternate  arrangements  are 
possible.  The  sheeting  may  be  fixed  to  the  supporting  structure 
on  all  four  sides  or  only  on  two  sides.  The  corrugations  on  the 
individual  panels  may  be  aligned  parallel  or  perpendicular  to  the 
applied  force. 

When  a  force  is  applied  to  a  diaphragm  a  distribution  of 
forces  among  the  different  elements  that  make  up  the  diaphragm 
takes  place.  By  means  of  fasteners  part  of  the  force  will  be 
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Hgure  2-1  Arrangement  of  typical  shear  panel 


9 

taken  by  end  and  intermediate  purlins  and  the  remainder  is 
carried  by  the  subpanels. 

Under  these  assumptions,  two  different  modes  of  failure 
of  the  diaphragm  may  occur.  First,  if  there  are  many  fasteners 
and  with  enough  strength  not  to  fail  before  buckling  of  the 
sheeting  occurs,  this  failure  that  is  called  "Buckling  as  an 
orthotropic  plate".  This  kind  of  failure  is  uncommon  and 
represents  the  upper  bound  load  that  a  steel  diaphragm  could 
take.  Figure  2-2  illustrates  such  a  buckling  failure.  Second, 
if  the  number  or  the  strength  of  the  fasteners  is  low,  the  steel 
deck  fails  at  loads  below  buckling  by  failure  at  the  fasteners. 
This  failure  may  occur  by  tearing  of  the  sheeting  at  the  fastener 
locations  or  by  simple  shear  failure  of  the  fasteners.  This  kind 
of  failure  is  called  "Failure  at  the  fasteners". 

2.4  ANALYTICAL  METHODS 

2.4.1  Failure  as  an  orthotropic  plate 

There  are  several  published  approaches  that  assess  the 
buckling  failure  of  a  diaphragm;  mainly  Bergmann  and  Reissner's 
(4)  (Appendix  A,  Eq.  A.l),  Hlavacek's  (21)  (Appendix  A,  Eq.  A. 4), 
Easley-McFarland ' s  (14)  (Appendix  A,  Eq.  A. 6),  Easley's  (15) 
(Appendix  A,  Eq.  A. 6).  All  of  them  give  similar  solutions  for 
the  critical  buckling  load.  Easley's  (15)  equations  (Appendix  A, 
Eq.  A. 6)  are  the  most  rigorous,  but  yield  solutions  only  when  the 
parameter  e  is  known;  e  is  only  available  for  simple  and  clamped 
supports.  Hlavacek's  equations  (20)  (Appendix  A,  eg.  A. 4)  are 
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Figure  2-2  Buckled  Pattern 


applicable  for  a  completely  different  kind  of  deck,  that  is  a 
flat  sheet  which  is  stiffened  by  perpendicular  plates  distributed 
symmetrically  to  both  faces  of  the  sheet. 

Bergmann  and  Reissner  (4)  were  the  first  investigating 
buckling  load  and  according  to  Easley  (15)  the  derivation  of 
their  formula  is  the  most  rigorous  of  all. 

2.4.2  Failure  at  the  fasteners 

A  number  of  researchers  have  provided  solutions  for  the 
failure  that  the  fasteners;  among  them  can  be  mentioned  Bryan  (7) 
and  Davies  (11)  in  England  being  the  pioneers,  Easley  (13)  in  the 
United  States  and  Fazio  (17)  et  al.  in  Canada.  The  most  relevant 
contributions  have  been  made  by  Bryan  (7),  Davies  (11)  and  Easley 
(13),  while  Fazio  (17)  et  al.  have  refined  Easley's  (13) 
method.  In  the  following  sections  a  description  of  the  solutions 
proposed  by  the  forement ioned  authors  is  presented  together  with 
the  advantages  and  disadvantages  of  their  methods.  Finally  a 
refined  approach  is  presented. 

2.4. 2.1  Bryan's  Method 

Bryan's  approach  (7)  has  been  accepted  as  a  feasible 
solution  for  steel  decks  in  Great  Britain.  His  method  is  based 
on  the  Energy  Method  and  Equilibrium.  Details  of  his  method  are 
presented  in  Appendix  B. 

Bryan  considers  the  flexibility  of  a  diaphragm  instead 
of  the  stiffness.  He  recognizes  two  types  of  steel  deck  and 
assumes  a  distribution  of  internal  forces  as  shown  in  Figure  2- 
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(b)  Sheet  Fasteners  on  two  edges  only 


Figure  2-3  Sheet  Fasteners 
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3.  If  shear  connectors  are  present  (Figure  2-3(a))  there  will  be 
slip  only  at  the  end  fasteners;  a  mode  of  deformation  that  he 
called  "Direct  Shear  Transfer".  Otherwise,  Figure  2-3(b) 

applies,  in  which  there  is  a  linear  distribution  of  forces  at  the 
end  fastener  locations.  Bryan  referred  to  this  as  the  "Indirect 
Shear  Transfer". 

Bryan  states  that  the  overall  flexibility  of  a  steel 
deck  is  the  result  of  the  summation  of  the  flexibilities  of  the 
components  that  make-up  a  steel  deck.  Consequently  he  analyzes 
the  sheet  deformation,  shear  strain  of  the  sheet,  axial  strain  in 
the  purlins  and  sheet  connection  deformations.  By  studying  the 
effects  separately,  the  influence  of  each  factor  in  the  overall 
flexibility  can  be  assessed  and  the  less  important  factors  may  be 
discarded. 

Bryan's  solution  cannot  be  considered  exact  because  it 
requires  the  determination  of  six  different  parameters  by 
testing.  Rigorously  his  solution  should  only  be  applicable  for 
the  kind  of  deck  and  connections  produced  by  the  different 
fabricators  in  Great  Britain;  when  applying  his  method  to  decks 
manufactured  in  North  America  testing  should  be  required  to 
obtain  new  parameters. 

Bryan's  assumed  distribution  of  internal  forces  is  not 
real  according  to  observations  made  on  tested  panels  (tilting  of 
the  sub-panels).  This  affects  the  flexibility  contribution  of 
the  connections  to  the  overall  flexibility  of  the  panel  but  does 
not  affect  the  results  for  the  axial  forces  in  the  purlins,  the 
sheet  deformation  and  the  shear  strain  of  the  sheet. 
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Bryan  did  not  take  into  account  the  vertical 
distribution  of  forces  that  exists  at  the  purlin  connections; 
this  is  an  item  that  may  be  incorporated  in  his  approach  to 
improve  his  method. 

The  advantage  of  his  method  is  that  the  expressions  for 
stiffness  are  not  difficult  to  calculate. 

Bryan  considers  different  possible  modes  of  failure 
assuming  an  internal  force  distribution  which  is  not  realistic. 
The  ultimate  strength  of  the  deck  is  obtained  by  considering 
failure  at  the  seam  connection,  or,  failure  at  the  shear 
connectors,  or,  failure  at  the  sheet-purlin  connection,  or, 
failure  at  the  end  purlins  by  buckling  whichever  is  least. 


2. 4. 2. 2  Davies'  Method 

Davies'  (11)  method  is  based  on  Bryan's  (7)  work.  He 
assumes  on  the  basis  of  Finite  element  analysis,  an  improved 
distribution  of  internal  forces.  Details  of  his  method  are 
presented  in  Appendix  C. 

Davies'  assumptions  can  be  discussed  best  with  reference 
to  Figures  2.4  and  2.5.  The  first  assumption  is  that  the  seam 
slip  at  each  internal  seam  is  equal  and  denoted  by  2A.  The 
second  assumes  that  the  total  horizontal  force  Qa/b  in  the  two 
edge  purlins  is  equally  shared  between  the  sub-panels.  Using 
these  two  assumptions  vertical  and  horizontal  equilibrium  are 
satisfied  for  an  interior  sub-panel.  When  the  same  principle  is 
applied  to  the  end  panel,  vertical  equilibrium  is  not  satisfied, 
so  that  he  makes  a  further  assumption:  the  slip  at  the  end 
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Figure  2-4  Forces  and  Deformations  Assumed  for  Typical  Internal  Sub-Panel 
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Figure  2-5  Forces  and  Deformations  Assumed  for  End  Sub-Panel 
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connectors  is  Ae  +  A.  This  Ae  is  obtained  from  equilibrium 
considerations  violating  moment  equilibrium  for  the  end 
sub-pane 1 . 

2. 4. 2. 3  Easley's  Method 

Easley's  method  (13)  is  basically  an  equilibrium 
approach  that  requires  some  parameters  that  are  obtained  by 
testing.  His  main  assumption  is  that  the  sub-panel  in  a 
diaphragm  tends  to  tilt  in  the  frame,  as  shown  in  Figure  2-6(a); 
this  was  also  confirmed  by  test  observations.  Davies  (11) 
applied  the  same  concept  to  his  theoretical  approach,  and 
MacGregor  (23)  testing  independently  at  the  University  of  Alberta 
made  the  same  observation.  Easley's  (13)  assumptions  are 
discussed  in  detail  in  the  following  paragraphs. 

The  assumption  that  sub-panels  deform  in  the  same  way  is 
a  simplifying  one  that  leads  to  simple  equilibrium  equations  for 
the  interior  sub-panels,  however  the  connections  used  are  the 
same  for  all  the  sub-panels;  this  is  an  idealization  not 
applicable  to  real  steel  deck  diaphragms. 

The  assumption  that  sub-panel  end  fastener  forces  have 
components  parallel  and  perpendicular  to  the  panel  ends  may  be 
easily  understood  with  the  aid  of  Figure  2.7,  which  shows  that 
panel  rotation  occurs  when  a  lateral  force  is  applied.  The 
magnitude  of  the  perpendicular  component  of  force  increase 
proportionally  to  the  distance  from  the  center  of  rotation.  This 
is  similar  to  an  assemblage  of  bolts  acted  upon  by  a  moment  where 
the  force  taken  by  any  bolt  is  proportional  to  the  bolt  distance 
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Figure  2-6  Fastener  Forces 


19 


from  the  center  of  rotation. 

Easley  assumes  that  the  distribution  of  forces  for  the 
sub-panel  side  fasteners  are  uniform  in  magnitude  and  parallel  to 
the  sides  of  the  sub-panel.  Theoretically  this  assumption  is  not 
correct.  When  the  diaphargm  starts  to  tilt  (Figure  2.7)  some 
forces  perpendicular  to  the  side  of  the  sub-panels  should  be 
developed.  Also,  as  it  has  been  shown  by  Finite  Element  Analysis 
by  Nilson  (28)  these  side  fastener  forces  are  distributed 
parabol ical ly . 

Easley  (13)  also  observed  in  his  testing  that  failure  of 
the  diaphragm  could  occur  because  of  local  bearing  failure  or 
tear  out  of  the  sub-panels  at  the  fasteners  either  those  which 
hold  sub-panels  together  or  those  connecting  sub-panels  to  the 
frames,  or  failure  in  the  fasteners  themselves.  He  recognized 
two  different  types  of  connection  failure,  either  failure  at  the 
end  fastener  located  in  line  with  the  seam  fastener  or  failure  at 
the  seam  fasteners.  The  first  assumption  of  failure  is 
theoretically  correct  because  it  corresponds  to  the  fastener 
taking  the  maximum  vertical  load  (see  Figure  2.6)  according  to 
the  linear  force  distribution.  Testing  has  confirmed  this 
assumption.  The  second  mode  of  failure  theoretically  corresponds 
to  all  the  seam  fasteners  in  a  sub-panel  failing 
simultaneously.  Testing  has  confirmed  this  assumption  too. 
Although  the  actual  distribution  of  forces  is  parabolic,  a 
redistribution  of  forces  takes  place  before  failure,  confirming 
that  the  behaviour  of  a  diaphragm  is  influenced  by  plastic 
deformations.  The  main  shortcoming  of  this  approach  is  that  it 


Figure  2-7  Shear  Panel  Rotation 
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is  not  suitable  to  the  type  of  steel  deck  usually  made  in  North 
America  in  which  fasteners  of  different  shear  stiffnesses  are 
used;  however,  his  equations  can  be  modified  to  include  fasteners 
that  have  different  load-deformation  characteristics. 

2. 4. 2.4  Fazio,  Ha,  El-Hakim's  Method 

This  approach  (17)  is  based  on  the  methods  of  Bryan  (7), 
Davies  (11)  and  Easley  (13)  and  it  has  the  same  shortcomings 
already  pointed  out.  Nevertheless,  unlike  Easley's  (13)  method, 
this  method  considers  equilibrium  of  an  external  sub-panel  using 
different  type  of  fasteners  for  the  end,  purlin  and  side 
attachments.  Therefore  the  purpose  of  Fazio,  Ha,  El-Hakim's 
approach  is  to  analyze  steel  decks  that  are  commonly  used  in 
practice.  The  main  assumption  made  is  that  the  stress 
distribution  at  the  end  fasteners  for  the  end  sub-panel  is  not 
symmetrical;  it  is  symmetrical  for  the  internal  sub-panels  (see 
Figure  2.8).  This  leads  to  simple  equations,  thus  making  the 
method  attractive  and  easy  to  apply.  However,  Fazio,  Ha,  El- 
Hakim's  assumption  satisfies  the  equilibrium  of  the  external  sub¬ 
panel  but  not  the  equilibrium  of  the  internal  sub-panels. 

2. 4. 2. 5  Writer's  Approach 

In  this  section  a  solution  to  the  stiffness  of  a  steel 
deck  diaphragm  has  been  developed  by  the  writer  and  is  discussed 
in  detail  in  the  following  paragraphs. 

The  assumptions  on  which  this  method  is  based  are: 

The  seam  slip  for  a  line  of  seam- fasteners  is 


1. 
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Figure  2-8  Deformation  of  the  Panel 
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constant,  but  it  varies  from  one  line  of  seam- 
fasteners  to  the  next  according  to  the  mode  of 
deformation  assumed  in  Figure  2.9. 

2.  The  vertical  component  of  the  shear  connector  force 
is  constant. 

3.  The  horizontal  component  of  the  shear  connector  and 
seam  fastener  force  is  negligible. 

4.  The  horizontal  force  taken  by  the  end  purlins  is 
shared  evenly  be  each  sub-panel 

5.  The  vertical  components  of  the  end  fastener  and  the 
purlin  fastener  forces  vary  linearly  at  every  sub¬ 
panel. 

For  every  sub-panel  there  are  two  unknown 


displacements.  Two  equations  are  available  for  each  sub-panel, 
from  vertical  force  equilibrium  and  moment  equilibirum,  which 
permit  the  determination  of  displacements  in  each  sub-panel. 
These  equations  are  solved  simultaneously  for  the  complete  panel 
to  obtain  the  total  displacement. 

The  necessary  equations  using  equilibrium  and 
compatibility  are  developed  below  with  the  aid  of  Figure  2.9. 

The  vertical  equilibrium  of  the  first  sub-panel  yields 
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Figure  2-9  Writer's  Approach 
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Also,  because  of  linear  distribution  of  stresses  within 
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Rearranging  terms 
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Moment  equilibrium  of  the  first  sub-panel  gives 
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Vertical  equilibrium  of  the  second  sub-panel  yields 
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Because  of  linear  distribution  of  stresses  within  the  sub-panel 
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Moment  equilibrium  of  the  second  sub-panel  yields 
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P 


n 


+  A.( 


nsh 

4'  S  a 
e 


n 


n 


sh 


2 

x .  +  _ 

i  S  a 


^  n  a  Ar.n  a  ^ 

2  .  s  N  ,  5  s  Qa 

x  •  +  c  —  )  +  c-  —  -  -  =  — 
i  Sn,  Sn  n  , 

s  sh  s  sh 


(17) 


Vertical  equilibrium  of  the  third  sub-panel  yields 


2  A, 

nF«-n  F  0  +  — — 
s  s2  sc  sc2  S 


(1 


x  . 

l 


Lo3 


PA, 


)  + 


Xi 

(i  -  — i) 

o3 


0  (18) 


sc2  S 


sc 


(19) 


and 


o3 


n  ,  o3 
sh 


(20) 
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or 


5  a 


o3  Ar  +  A,  n  , 
5  6  sh 


(21) 


Replacing  12 ,  19  and  20  in  18  results 


n  n  2n  Pn  2n  , 

*  s  ,  A  ,  s  ,  e  ,  p  sh 

4b  D  b  b  b  ba 

s  s  e  p  e 


Pn 


x 


sh 


S  a 
P 


xi} 


n 


n 


n  2n  , 

A  +  Sh 


6  v  S 


sc 


S  a 
e 


H  xi 


Pn 


sh 


S  a 
P 


x  i )  =  0 


(22) 


Moment  equilibrium  yields 


n  A 


sc  6  S 


sc  n 


z 


sh 


x .  P  Ac 

(1  -  3TT)xi  "  — 
°3  p 


(1  -  -!t)x.  = 

x  _  1  n  , 

o3  sh 

(23) 


Replacing  21  in  23. 


2n 


Ac( 


sh 


5  S  a 
e 


Z-i 


~  Pn  . 
2  sh 

S  a 
P 


z«,j  - » 


s  ^_xi>  (24) 

e  p 


+  A -i-  + 


6  S  nsh  S 
sc  e 


nsh  \ 

7^  ilxi 


2  +  nsh 

S  a 
P 


2 .  Qa 
xi  >  =  7T7 

sh 
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where : 

FSi 

=  Force  taken  by  a  seam  fastener  at 

line  i 

Fsci 

=  Force  taken  by  a  shear  connector 

at  edge 

ne 

=  number  of 

end  fasteners 

nP 

=  number  of 

purlin  fasteners 

ns 

=  number  of 

seam  fasteners 

n 

P 

S, 


sc 


S 


sc 


‘Ol 


number  of  shear  connectors 
number  of  intermediate  purlins 
flexibility  of  an  end  fastener 
flexibility  of  a  purlin  fastener 
flexibility  of  a  seam  fastener 
flexibility  of  a  shear  connector,, 

distance  from  the  left  edge  of  sub-panel  i  to  the  end 
and  purlin  fastener  as  shown  in  Figure  2-9 
distance  from  the  left  edge  of  sub-panel  i  to  the 
neutral  axis  of  end  and  purlin  fastener  as  shown  in 
Figure  2-9 

displacement  of  left  end  fastener  of  sub-panel  i 


The  same  procedure  may  be  applied  for  more  than  3  sub¬ 
panels. 

Once  the  individual  displacements  are  obtained,  this 
component  of  the  shear  stiffness  can  be  added  to  the  sheet  shear 
def ormat iona 1  stiffness  and  the  stiffness  resulting  from  axial 
strain  in  purlins,  as  derived  by  Bryan  (7). 

The  advantages  of  this  method  over  the  previously 
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mentioned  methods  may  be  summarized  as  follows:  While  Easley 
(13)  considers  only  one  type  of  connections,  this  method  makes 
use  of  different  types,  which  is  more  applicable  to  real  steel 
decks.  Fazio's  (17)  approach  assures  equilibrium  of  the  external 
sub-panel  but  does  not  satisfy  equilibrium  for  the  interior 
ones.  Davies'  approach  violates  moment  equilibrium  for  the 
external  sub-panel.  The  writer's  approach  meets  moment  and  force 
requirements  for  every  sub-panel.  The  only  shortcoming  of  the 
writer's  approach  is  that  it  requires  a  computer  program  to  solve 
the  set  of  simultaneous  equations. 

2.5  THE  EMPIRICAL  METHOD  OF  THE  U.S.  DEPARTMENT  OF  ARMY, 

NAVY  AND  AIR  FORCE 

The  equations  used  by  the  U.S.  Department  of  Army,  Navy 
and  Air  Force  (33)  (Appendix  F)  have  an  empirical  basis.  They 
are  the  result  of  information  gathered  by  the  Triservice  Seismic 
Design  Committee  and  received  from  the  companies  supplying  the 
material.  According  the  the  U.S.  Department  of  Army,  Navy  and 
Air  Force  (33)  the  equations  for  strength  and  stiffness  are 
subjected  to  some  revision  as  new  data  is  obtained,  which  makes 
the  method  not  applicable  for  a  particular  kind  of  deck,  but 
only,  for  decks  having  similar  characteristics  to  the  ones 
tested . 

This  method  is  the  only  one  that  provides  equations  to 
calculate  the  stiffness  and  strength  for  steel  concrete 
diaphragms.  It  is  illustrated  in  one  example  in  Chapter  V. 
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2.6  COMPARISON  OF  ANALYTICAL  METHODS 

Only  a  limited  amount  of  published  information  regarding 
tests  performed  on  steel  decks  is  available  in  the  literature.  A 
series  of  tests  by  Davies  (11)  were  selected  to  assist  in  a 
comparison  and  assessment  of  each  of  the  applicable  analytical 
methods  described  above.  Calculations  for  deck  flexibilities 
were  made  and  are  presented  in  Table  2.1,  together  with  the 
experimental  results  reported  by  Davies  (11).  Parameters  common 
to  all  the  methods  are  shown  in  Table  2.2. 

Figures  2-10  and  2-11  show  the  dimensions  and  general 
arrangement  of  steel  decks  tested  by  Davies  (11). 

Basically  panels  A  and  B  are  similar  using  an  open 
corrugated  profile  (Amascodeck  R15,  Ref.  7). 

Panels  A1  and  A2  are  of  direct  shear  transfer,  that  is, 
with  shear  connectors.  While  A1  has  purlin  fasteners  and  end 
fasteners  in  every  corrugation;  A2  has  fasteners  only  in 
alternate  corrugations. 

Panels  B1 ,  B2  and  B3  are  of  indirect  shear  transfer  that 
is,  without  shear  connectors  and  with  purlin  and  end  fasteners  at 
every  corrugation,  at  alternate  corrugations  and  at  every  third 
corrugation  respectively. 

Panel  Cl  and  C2  are  similar,  without  shear  connectors 
and  with  purlin  and  end  fasteners  at  every  corrugation  and  at 
alternate  corrugations,  respectively. 

From  the  results  summarized  in  Table  2.2,  it  is  observed 
that  the  overall  flexibility  depends  heavily  on  two  parameters, 
sheet  distorsion  and  the  fastener  flexibility.  These  two 
parameters  account  for  at  least  70%  of  the  total  flexibility  of 


the  steel  deck. 


Table  2.1  Summary  of  Panel  Flexibilities  in  inch/kip 
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Figure  2-10  Typical  Shear  Panel  tested  (Types ABB) 


35 


Figure  2-11  Typical  Shear  Panel  Tested(Type  C) 


Table  2.2  Flexibility  Factors 
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The  sheet  distortion  depends  on  an  experimental  constant 
KQ  Eq.  B.l,  Appendix  B  which  should  be  carefully  assessed.  For 
panel  B1  the  test  value  of  flexibility  is  not  predicted  by  any  of 
the  methods;  suggesting  that  the  test  value  is  not  accurate. 

The  results  given  by  the  writer's  method  predict  the 
flexibility  of  the  diaphragm  with  almost  the  same  degree  of 
accuracy  as  that  of  the  other  methods.  Furthermore  not  having 
the  theoretical  shortcomings  of  Fazio,  Ha,  El  Hakim's  (17)  and 
Davies'  (11)  method  regarding  internal  equilibrium  and  moment 
respectively,  makes  the  method  more  acceptable. 

The  distribution  of  forces  for  the  end  fasteners  is  more 
realistic  than  in  Bryan's  (7)  method. 

Bryan's  (7)  approach  for  the  prediction  of  the 
flexibility  of  panels  Cl  and  C2  results  in  calculated 
flexibilities  that  are  greatly  different  from  test  results, 
which  points  out  that  his  theory  for  indirect  shear  transfer  is 
not  completely  correct. 

2.7  PARAMETERS  THAT  INFLUENCE  STRENGTH  AND  STIFFNESS  OF  A 

STEEL  DECK  DIAPHRAGM 

2.7.1  Introduction 

When  designing  a  building  in  Engineering  practice  it  is 
necessary  to  know  the  strength  and  the  stiffness  of  each  of  the 
elements  that  shape  it.  For  a  steel  deck  diaphragm  under  in¬ 
plane  loads,  there  are  several  variable  such  as,  panel  length, 
panel  width,  purlin  spacing  and  number  and  arrangement  of 
connections,  that  affect  the  values  of  the  strength  and 

Thus,  to  obtain  an  optimum  design,  the  influence  of 


stiffness. 
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each  variable  on  the  overall  strength  and  stiffness  should  be 
known. 

The  Fazio,  Ha,  El  Hakims 's  (17)  approach,  being  one  of 
the  most  accurate  procedures  according  to  the  comparisons  made  in 
the  preceding  section  has  been  selected  to  show  how  panel  length, 
panel  width,  purlin  spacing  and  number  and  arrangement  of 
connections  affect  the  behaviour  of  a  steel  panel. 

Two  computer  programs  based  on  Fazio,  Ha,  El  Hakim's 
(17)  equations  were  developed;  the  first  one  for  strength 
calculations  and  the  second  one  for  stiffness  calculations.  As 
will  be  shown  in  Chapter  III,  the  ratio  of  frame  to  slab 
stiffness  needs  to  be  known  in  order  to  calculate  the  amount  of 
force  taken  by  intermediate  frames. 

2.7.2 _ Panel  Length  and  Width 

Deck  T15-24  was  selected  from  reference  35.  Figures 
2~12(a)  and  2-12(b)  show  the  arrangement  of  connections,  purlins 
and  the  profile  selected.  For  the  calculations  performed  the 
connection  properties  were  kept  constant  as  follows: 

Ss  =  0.0613  in/kip  (flexibility  of  seam  fasteners) 

S  =  0.0613  in/kip  (flexibility  of  purlin  fasteners) 

sr 

=  0.02049  in/kip  (flexibility  of  shear  connectors) 
se  =  0.0613  in/kip  (flexibility  of  end  fasteners) 

Fsu  =  0.382  kips  (strength  of  seam  fasteners) 

Fpu  =  0.917  kips  (strength  of  purlin  fasteners) 

Fsc  =  2.751  kips  (strength  of  shear  connectors) 

=  0.917  kips  (strength  of  end  fasteners) 


Feu 
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30'  b 


(b) 


h - a  - -i 

Figure  2-12  Typical  Arrangement 


ko 


St i ffness  -  The  results  of  calculations  are  tabulated  in 
Tables  2.3,  2.4  and  2.5  for  30',  22.5*  and  15s  width  panels, 
respectively  and  are  plotted  in  Figure  2.13. 

The  first  observation  that  may  be  drawn  from 
Figure  2.13  is  that  there  is  a  well  defined  relation  between 
stiffness  and  length  of  a  shear  panel;  a  power  curve  fits  the 
data  points  well.  The  coefficient  of  determination  ranges 
between  0.976  and  0.993.  The  power  curves  that  best  represent 
the  data  points  ares 
For  30*  width 

K  =  421.4  x  2~le006(k/in) 

For  22 . 5 8  width 

K  =  259.7  x  2”1 e  027 ( k/in ) 

For  15"  width 

K  =  111.1  x  2”1 ' 038 ( k/in ) 

A  second  conclusion  that  can  be  drawn  is  related  to 
Tables  2.3,  2.4  and  2.5.  It  is  seen  that  the  two  most  important 
parameters  that  effect  the  flexibility,  or  its  inverse, 
stiffness,  of  a  steel  deck  diaphragm  are  sheet  distortion  ( c ^  ^ , 
Appendix  B,  eq.  B.l)  and  the  flexibility  of  the  connections. 

When  designing,  control  of  these  two  parameters  is  important. 
Therefore,  if  an  accurate  value  of  the  shear  stiffness  is 
required,  care  should  be  taken  to  ensure  that  experimental 
factors,  such  as  shear  stiffness  of  fasteners  and  KQ  factor  for 
sheet  distortion,  involved  in  the  determination  of  shear 
stiffness  are  properly  obtained. 

Strength  -  To  analyze  the  effect  of  different  widths  on 


Table  2.3.  Width  30'  Length's  Influence  in  Stiffness  for  Width  -  30' 


k) 


T=»ble  2.4.  Width  22.5'  Length's  Influence  in  Stiffness  for  Width 
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Figure  2~I3  Stiffness  of  Diaphragms 
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Figure  2-14  Strength  of  Diaphragms 
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Table  2.6  Width's  Influence  in  Strength 


Panel  Width  (Feet) 

30 

22.5 

15 

Failure  at  Seam  (kips) 

12.05 

9.26 

6.45 

Failure  at  Side  (kips) 

26.99 

20.49 

13.96 

Failure  at  End  (kips) 

22.55 

17.06 

11.55 

Failure  Load  (kips) 

12.05 

9.26 

6.45 

strength  the  same  panels  discussed  previously  were  used.  The 
results  are  shown  in  Table  2.6  where  the  lowest  value  for  each 


47 


panel  is  the  failure  load.  If  these  values  are  plotted  Figure 
2.14  is  obtained,  which  shows  a  linear  increase  in  strength  with 
width. 

2.7.3  Purlin  Spacing 

The  same  panel  used  in  the  previous  section  was  selected 
for  the  purpose  of  analyzing  the  effect  of  purlin  spacing  on 
strength  and  stiffness.  The  dimensions  of  the  panel,  profile, 
purlin  section,  number  of  seam  and  end  fasteners,  and,  shear 
connectors  were  kept  constant  throughout  the  analysis.  In  Figure 
2-15  the  different  arrangements  selected  are  shown.  The  only 
difference  between  any  two  of  them  is  the  number  of  purlins  and 
purlin  fasteners. 

Stiffness  -  The  results  are  shown  in  Figure  2-16.  It  is 
observed  that  stiffness  and  number  of  purlins  are  linearly 
related.  Another  observation  is  that  the  bigger  the  ratio  a/b 
(Length/width)  the  less  important  it  is  to  add  more  purlins  to 
the  panel.  From  Figure  2-16,  the  slope  of  the  straight  lines 
decreases  when  a/b  increases. 

Strength  -  As  for  stiffness,  the  strength  increases  as 
the  number  of  purlins  increased  as  shown  in  Figure  2-17. 

2.7.4  Number  of  Seam  Connectors 

The  effect  of  increasing  the  number  of  seam  connectors 
is  analyzed  for  shear  panels  of  different  sizes.  The  result  of 
the  calculations  are  shown  in  Figure  2-18  and  2-19. 
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Figure  2-15  Influence  of  purlin  spacing 
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Figure  2-16  Influence  of  Purlin  spacing  on  stiffness 
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Figure  2-17  Influence  of  Purlin  spacing  on  strength 
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Figure  2~\Q  Influence  of  the  number  of  Seam  fasteners  on  strength 
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Figure  2-19  Seam  Connection's  Influence 
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Figure  2-18  shows  the  influence  of  the  number  of  seam 
connectors  on  strength.  It  is  observed  that  a  straight  line  fits 
the  data  points  very  well,  and  that  the  effect  of  increasing  the 
number  of  seam  connectors  is  the  same  for  different  shear  panel 
lengths . 

Figure  2-19  shows  the  influence  of  the  number  of  seam 
fasteners  on  stiffness.  A  small  increase  in  the  number  of  seam 
fasteners  implies  a  big  increase  in  stiffness  and  this  influence 
is  less  when  the  length  of  the  shear  panel  increases. 

2.7.5  Number  of  Shear  Connectors 

The  influence  of  increasing  the  number  of  shear 
connectors  is  studied  for  shear  panels  of  different  dimensions. 
The  results  are  shown  in  Figure  2-20,  2-21  and  2-22. 

Figure  2-20  shows  the  influence  of  the  number  of  shear 
connectors  on  strength.  It  is  seen  that  only  a  few  connectors 
are  needed  in  order  to  reach  a  strength  value  above  which  any 
increase  in  the  number  of  connectors  is  useless  from  the  strength 
point  of  view. 

Figures  2-21  and  2-22  show  the  influence  of  the  number 
of  shear  connectors  on  stiffness.  The  same  observation  as  that 
obtained  for  seam  connectors  may  be  pointed  out,  that  is,  a  small 
increase  in  the  number  of  shear  connectors  implies  a  big  increase 
in  stiffness  and  this  influence  is  less  when  the  length  of  the 
shear  panel  increases. 
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Figure  2-20  Shear  Connector's  Influence  on  Strength 
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Figure  2-21  Shear  Connector's  Influence  on  Stiffness 
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Figure  2-22  Shear  Connector's  Influence  on  Stiffness 


CHAPTER  III. 


INTERACTION  FRAME-DIAPHRAGM  FOR  SYMMETRICAL 
ONE-STOREY  BUILDINGS. 

3.1  INTRODUCTION 

This  chapter  is  devoted  to  the  analysis  of  rectangular, 
symmetrical,  one  storey  buildings  subjected  to  uniformly 
distributed  horizontal  loads  such  as  wind  or  earthquake. 

Typical  buildings  analyzed  are  shown  in  Figure  3-1 
(Types  1,  2  and  3).  They  all  have  end  shear  walls  and  parallel 
intermediate  frames.  The  only  difference  between  two  of  them  is 
the  number  of  intermediate  frames.  Theoretically  the  number  of 
intermediate  frames  may  vary  from  zero  to  infinity. 

A  typical  building  is  shown  in  an  isometrical  view  in 
Figure  3-2  and  is  comprised  of  several  elements:  a  front  wall 
ADKI,  a  horizontal  diaphragm  AHDE,  two  interior  frames  BGMN  and 
CFOP  and  two  end  walls  AHKL  and  DEIJ. 

When  horizontal  forces  act  on  a  building  such  as  the  one 
shown  in  Fig.  3-2  a  complex  distribution  of  forces  among  the 
structural  elements  that  make-up  the  building  takes  place.  The 
lateral  horizontal  load  is  taken  by  the  front  wall  ADKI  (Figure 
3-2);  if  this  wall  is  structurally  strong  it  transfers  part  of 
the  force  to  the  horizontal  diaphragm  AHDE,  which  in  turn  in 
combination  with  shear  walls  AHKL  and  DEIJ  and  interior  frames 
BGMN  and  CFOP  finally  transfers  the  lateral  loads  to  the 
foundation  of  the  building. 
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View  of  a  Building 


Figure  3-2  Isometric 
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3.2  Method  of  Analysis 

The  structural  behaviour  of  a  building  may  be  studied  by 
means  of  models.  To  represent  the  horizontal  diaphragm,  two 
models  have  been  selected:  the  beam  model  and  the  truss  model. 

To  model  frames  and  walls,  springs  of  different  constants  are 
used. 

For  buildings  up  to  five  bays,  the  beam  model  yields 
solutions  which  are  not  difficult  to  obtain.  For  more  than  five 
bays  a  solution  based  on  the  theory  of  beams  on  elastic 
foundations  which  was  first  used  by  Muto  (25),  is  developed. 

3.2.1  Beam  Model 

Generally,  the  horizontal  diaphragm  will  behave  as  a 
deep  girder,  so  that  both  bending  and  shear  properties  must  be 
considered.  However,  when  the  building  has  a  steel  deck 
diaphragm  the  beam  model  is  adopted  using  the  shear  properties 
only,  because  the  contribution  of  the  edge  girders  and  the  steel 
deck  to  the  bending  properties  is  negligible.  When  the 
horizontal  diaphragm  is  made  of  concrete,  bending  as  well  as 
shear  properties  must  be  considered. 

The  beam  model  procedure  was  applied  to  determine  the 
force  distribution  for  frames  tested  by  Bryan  (7).  Table  3-1 
presents  the  experimental  and  analytical  results  reported  by 
Bryan  together  with  those  obtained  from  the  application  of  the 
beam  model.  In  this  case  only  shear  properties  are  used  because 
the  horizontal  diaphragm  is  a  steel  deck.  As  shown  in  Table  3.1, 
the  beam  model  predicts  the  force  taken  by  the  frames  well. 
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Table  3 . 1 
STIFFNESSES 


Tests 

kips/in 

Beam  Model 
kips/in 

% 

of 

Tests 

Bryan 

kips/in 

% 

of 

Tests 

TEST  I 

No  Gable 

FI 

10 . 7 

9.1 

85 

8.7 

81.3 

Bracing 
at  ends 

F  2 

11.8 

10.5 

89 

10.2 

86.4 

TEST  I 

Gable 

Bracing 

FI 

6.91 

5.81 

84.1 

5.79 

83.8 

at  Ends 

F2 

8.08 

7.58 

93.8 

7.57 

93.7 

TEST  II 
Gable 

FI 

2.76 

2.74 

99.3 

Bracing 

At  Ends 

F2 

3.55 

3.62 

102 

3.60  101.4 
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3.2.1. 1  Solution  for  a  Limited  Number  of  Bays 

The  solution  for  the  force  taken  by  an  interior  frame, 
on  a  three  bay  building,  (Type  1,  as  shown  in  Figure  3-1)  is  used 
to  explain  some  of  the  concepts  involved  in  the  beam  model. 
Derivations  are  detailed  in  Appendix  J. 

The  force  taken  by  a  frame  is: 


V 


R 


5 

K 


K, 
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K 


+ 


sh 


V 


K, 


K, 


VT  KF 

2*725  kT+  3  FT 

B  sh 


+ 


kw//kf 


+  1 


where : 


(3.1) 


VR  =  Force  taken  by  an  interior  frame 

VT  =  Total  horizontal  force  acting  on  the  building 

Kp  =  Stiffness  of  the  frame 

Kq  =  Bending  stiffness  of  the  horizontal  diaphragm 
Ks^  =  Shearing  stiffness  of  the  horizontal  diaphragm 
K w  =  Stiffness  of  the  wall. 

Simplifying  further,  the  following  quantities  are 
replaced  in  Equation  3.1. 

Kp/Ks  =  x 

Kp/Ksh  =  y 

Kw/Kp  =  c  (in  general  c  =  EK^/EKp) 


63 


Thus,  the  solution  becomes: 


V  x  +  y  +  V2  c 

^  - - T -  (3.2) 

c 

Using  Fig.  3.3,  Equation  (3.2)  represents  for  every  c 
value  a  space  surface  whose  intersection  with  the  planes  y  =  0 
and  x  =  0  yield  the  curves  for  the  pure  bending  solution  and  the 
pure  shear  solution,  respectively. 

Since  factor  c  may  vary  from  »  to  0  representing  the 
lower  and  upper  bounds  respectively,  the  limits  for  the  force 
taken  by  a  frame  would  be: 

Lower  bound  limit 


VR  x  +  y 

VT  2.725x  +  3y  +  1 

c  ->•  00 

which  is  a  space  surface 
of  this  surface  with  the 
curves : 


(3.3) 

shown  in  Figure  3.3.  The  intersection 
bending  and  shear  planes  yields  the 


v 


V 

V 


R 


X 

2.725x  +  1 

y 

3y+l 


and 


The  Upper  bound  limit  is: 


Lim 


0.5 


(3.4) 


c  -►  0 
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Figure  3-3  Schematic  Interaction 
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which  is  a  plane  in  the  space  as  shown  in  Figure  3.3.  Equations 
(3.3)  and  (3.4)  are  the  limits  for  the  complete  set  of  solutions. 

The  plane  x  =  0  is  the  pure  shear  plane  in  that  only 
shear  properties  are  present,  similarly  the  plane  y  =  o  is  the 
pure  bending  plane  since  only  bending  properties  are  present. 

Any  plane  passing  through  the  origin  and  perpendicular  to  the 
plane  z  =  0  yields  a  combined  solution  in  which  bending  as  well 
as  shear  properties  are  present  in  a  constant  ratio. 

The  intersection  of  the  space  surface  (3.2)  with  the 
plane  y  =  x  is  an  example  of  a  combined  solution  as  shown  in 
Figure  3.3  and  given  by  Eqn.  3.5  below. 


V 


R 


2x  +  V? 


VT  5.725x  +  -  +  1 

c 


(3.5) 


The  lower  bound  limit  in  this  plane  is 


Lim 


V 

V„ 


R 


2x 


5 . 7  25x  +  1 


(3.6) 


c  •+■  °° 


which  coincides  with  equation  (3.3)  when  y  is  replaced  by  x. 

Similarly  the  upper  bound  limit  in  this  plane  is: 


V 

L  ^5  =  0.5  (3.7) 

im  VT 

c-*0 


Because  the  pure  shear  plane  and  the  pure  bending  plane 
are  the  limits  where  either  only  shear  properties  or  bending 
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properties  are  represented,  it  is  important  to  have  the  solution 
for  the  force  taken  by  an  interior  frame  for  these  planes. 

Figures  3.4  and  3.5  show  solutions  for  the  force  taken 
by  interior  frames  for  a  three  and  a  four  bay  building 
respectively.  In  these  figures  the  pure  shear  solution  and  the 
pure  bending  solution  are  presented  together.  To  obtain  a  wide 
range  of  stiffness  the  x  axis  is  necessarily  logarithmic.  The  y 
axis  represents  the  percentage  of  the  total  horizontal  force 
taken  by  an  interior  frame.  The  parameter  c  has  already  been 
defined  and  the  parameter  n  =  1  identifies  symmetrical  buildings 
as  will  be  seen  in  the  following  chapter. 

As  it  is  shown  in  figures  3.4  and  3.5,  the  curves 
representing  the  fraction  of  the  total  horizontal  force  taken  by 
an  interior  frame,  range  between  two  extremes  which  correspond  to 
a  relatively  infinitely  stiff  diaphragm  and  an  infinitely 
flexible  diaphragm. 

When  the  horizontal  diaphragm  is  stiff  relative  to  the 
frame,  the  ratio  Kp/Kg  (stiffness  of  the  frame/stiffness  of  the 
slab)  tends  to  zero  and  all  the  curves  flatten  out;  this  means 
that  for  any  c  the  displacement  of  the  horizontal  diaphragm  is 
constant  along  the  building  and  the  force  taken  by  an  interior 
frame  is  proportional  to  its  stiffness.  As  the  ratio  Kp/Kg 
increases,  the  fraction  of  the  total  horizontal  force  taken  by  an 
interior  frame  increases,  so  that  if  it  is  assumed  that  the 
horizontal  diaphragm  is  infinitely  stiff  the  real  force  taken  by 
a  frame  could  be  underevaluated.  Otherwise,  if  it  is  assumed 
that  the  horizontal  diaphragm  is  infinitely  flexible,  this  leads 
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Figure  3-4  Three  Bays  Symmetrical 
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Figure  3-5  Four  Bays  Symmetrical  Loads 


69 


to  an  overevaluation  of  the  force  taken  by  the  frame. 

3. 2. 1.2  Truss  Model 

This  model  is  only  applicable  to  buildings  having  steel 
deck  diaphragms  and  is  equivalent  to  the  pure  shear  solution 
given  by  the  beam  model  as  it  is  explained  in  Appendix  H,  for  the 
case  of  a  two  bay  building. 

As  shown  in  Figure  3.6,  for  a  type  1  building,  the 
horizontal  diaphragm  is  modelled  by  a  diagonal  bar  AB  in  which 
the  shear  deflection  is  the  same  as  in  the  diaphragm. 

The  vertical  and  horizontal  members  that  shape  the  truss 
are  taken  as  infinitely  rigid,  walls  and  intermediate  frames  are 
replaced  by  springs. 

STIFFNESS  OF  THE  HORIZONTAL  DIAPHRAGM  —  Figure  3.6  shows  3  types 
of  buildings.  The  steel  deck  diaphragm  of  a  typical  bay  of  any 
of  those  buildings  is  replaced  by  a  diagonal  bar  AB  that  gives  an 
equivalent  shear  deflection. 

A 

Figure  3.7  shows  the  same  typical  bay  under  the  action 
of  a  lateral  force  K-^ 


K1  =  P1  cos  0 


but  cos  0  =  b/h} 

so  that  P-|_  =  h^/b 

hl 

Triangles  ACB  and  BDE  are  similar  thus 
so  that 


1 

5 
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Type  1 


Type  3 


Figure  3-6  Truss  Model 
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Figure  3~7  Stiffness  of  the  Diaphragm 
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also 


<5  = 


Plhl 

AE 


Finally 


,  Ki  h7 

b  11 

h^  AEb 


thus  the  force  that  causes  an  unitary  displacement  is 


K 


b2AE 


1  h3 

hl 


(3.8) 


and  the  area  of  the  equivalent  bar  is 


A  = 


Kihi 

b2E 


(3.9) 


Three  bays  buildings 

The  layout  of  the  building  is  shown  in  Figure  3.7  (Type 
1).  It  consists  of  three  bays,  two  end  shear  walls  and  two 
interior  frames.  This  distribution  may  be  interchanged  for  two 
end  frames  and  two  interior  shear  walls. 

The  variables  involved  in  the  solution  of  the  problem 
were  already  defined  for  the  beam  model  and  they  are: 


s 


73 


which  is  the  independent  variable. 


The  dependent  variable  is 


Z 


and  the  parameter  c  is  defined  as: 


c  =  2  KW/£KF 

These  above  mentioned  variables  are  related  by  the  set 
of  curves  shown  in  Figure  3.8. 

It  is  observed  from  Figure  3.8  that  there  is  a  well 
defined  relation  between  the  relative  stiffness  Kp/Ks  and  the 
percentage  of  horizontal  force  taken  by  the  frame.  Depending  on 
the  value  of  the  c  constant,  it  can  be  observed  that  there  is  an 
upper  boundary  when  the  external  walls  are  frames  represented  by 
c  =  1.  When  the  parameter  c  increases  the  curves  get  closer  and 
most  of  them  gather  in  a  narrow  band  space,  so  that,  from  a 
practical  viewpoint  there  is  no  difference  for  the  force  taken  by 
the  frame  when  the  outerwalls  are  either  concrete  or  masonry. 

An  important  observation  that  may  be  obtained  from 
Figure  3.8  is  that  the  assumption  of  an  infinitely  rigid 
horizontal  diaphragm  may  lead  to  erroneous  results  for  the  force 
carried  by  the  frame.  A  value  of  the  relative  stiffness  can  be 
chosen  such  that  the  assumption  of  a  rigid  horizontal  diaphragm 
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is  acceptable,  for  example 


V 

V 


R 

Too 


0.20 


This  condition  is  satisfied  for  a  value  of  K^/k 

r  S 

approximately  equal  to  10-2,  thus  Figure  3.8  becomes  Figure  3.9. 

It  is  also  observed  that  when  the  stiffness  of  the  frame 
relative  to  the  stiffness  of  the  horizontal  diaphragm  is  large, 
all  the  curves  gather  in  a  narrow  space  and  finally  for  Kp/Ks  > 

40  the  curves  reach  a  constant  value  equal  to  0.352  for 
infinitely  flexible  diaphragms  which  is  higher  than  the  area 
distribution  method  which  yields  a  value  0.333. 

Four  and  Five  Bays  One  Storey  Buildings 

The  features  of  these  buildings  are  shown  in  Figure  3.6 
(Type  3  and  Type  2  for  four  and  five  bays,  respectively). 

The  same  model  is  used  as  for  three  bays,  but  in  these 
two  cases  for  each  c  value  two  curves  are  obtained  which 
represent  the  fraction  of  the  total  horizontal  force  taken  by 
each  of  the  interior  frames.  (See  Figures  3.10,  3.11  and  3.12) 
(only  one  curve  is  shown  for  the  case  of  four  bays) 

From  these  figures  is  observed  again  that  there  is  a 
well  defined  relation  between  the  relative  stiffness  and  the 
percentage  of  the  total  horizontal  force  taken  by  an  interior 


frame 


77 


78 


79 


PZO 


000 


> 


ac 


<p 

e 


0*0 


E20 


91  0 


800 


Figure  3- i 2  Four  Bays  Truss  Model  Percentage  of  force  token  by  f 
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3.2.2  General  Solution 

The  general  solution  is  based  on  the  theory  of  beams  on 
elastic  foundation  as  first  used  by  Muto  (25). 

As  shown  in  Figure  3.13,  the  stiffness  of  the  frames  is 
converted  into  resisting  elements  distributed  per  unit  length  of 
the  building.  Development  of  the  equations  is  presented  in 
Appendix  N. 

The  differential  equation  that  represents  this  model  is: 


W 
E  I 


s  s 


(3.10) 


0 


3n< 

(--) 

2G  A 
s  s 


where : 


^s 

gF 

G„ 


W 


Shear  area  of  the  horizontal  diaphragm 
Young's  Modulus  of  the  horizontal  diaphragm 
Uniformly  distributed  resisting  force 
Shearing  Modulus  of  the  horizontal  diaphragm 
Moment  of  Inertia  of  the  horizontal  diaphragm 
Uniformly  distributed  horizontal  load 
Displacement  of  the  horizontal  diaphragm 
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Figure  3-13  Beam  on  an  Elastic  Foundation 
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3q  =  Coefficient  of  inelastic  shear  deformation 
<  =  Shear  deformation  coefficient. 

Solutions  for  the  differential  equation  are  dependent  upon  the 
factors  «  and  3. 
if  «  >  3 

y  =  edx(Acosfx  x  Bsinfx)  +  e  dx(Ccosfx  +  Dsinfx)  +  W/gp 

(3.11) 


where : 

d 

f 

If  *  =  3 

y  =  A  edx  +  Bxedx  +  C  e”dx  +  Dxe~dx  +  W/gp  (3.12) 

If  «  <  3 

.  (d  +  f ) x  ^  _  - ( d+f ) x  (d-f)x  -(d-f)x  . 

y  =  Ae  +  Be  +  Ce  +  De  +  w/gp 

(3.13) 


2  2 
p.-tB.:) 


v2 


2  .2 


V2 


The  Boundary  conditions  are 


y 


°w//sw ' 


x  =  1/2, 
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where 

Ms  =  Bending  moment  of  the  horizontal  diaphragm 
Qw  =  Force  taken  by  an  end  wall 
Sw  =  Stiffness  of  an  end  wall 

The  solution  is  developed  in  Appendix  No  A  computer 
solution  is  also  given  in  Appendix  T. 

This  procedure  may  be  applied  for  different  boundary 
conditions,  that  is,  end  shear  walls,  shear  wall  at  the  middle  of 
the  building  or  eccentrically  arranged  shear  walls. 

3.2.3  Comparison  of  the  General  Solution  with  a  Standard 

Program 

Several  computer  runs  were  made  in  order  to  evaluate  the 
results  given  by  the  general  solution  and  a  standard  program 
(SAPIV)  (31)  for  a  3  bay  building  and  a  5  bay  building.  When 
using  the  standard  program,  the  horizontal  diaphragm  is 
represented  by  a  beam  element,  while  the  vertical  resisting 
elements  are  replaced  by  spring  elements.  For  the  general 
solution  the  intermediate  frames  are  replaced  by  a  continuous 
resisting  element. 

The  results  of  these  computer  solutions  are  shown  in 
Tables  3.2  an  3.3  The  percentage  of  difference,  between  both 
solutions  for  the  percentage  of  the  total  horizontal  force  taken 
by  the  interior  frames  ranges  between  0%  and  18%  for  three 
bays.  For  five  bays  differences  range  between  0.2%  and  16%. 

The  main  reason  for  these  differences  is  related  to  the 
different  models  used  to  represent  the  frames.  For  longer 


buildings 

with  a  larger 

number  of  bays 

the  two  models 

should  give 

similar  results. 

/ 
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Table  3.2 

THREE  BAY  BUILDING 


Case 

Kp/Ks 

vF/vT 

% 

Di f  f erence 

MUTO  1 

108.1 

0.00106 

0.00530 

5.36 

SAP  1 

0.00560 

MUTO  2 

105.6 

0.00109 

0 . 00546 

5.04 

SAP  2 

0.00575 

MUTO  3 

1.0 

0.00109 

0.25018 

0.00 

SAP  3 

0 . 25017 

MUTO  4 

108 . 1 

4 . 2500 

0 . 36553 

7 . 09 

SAP  4 

0.33961 

MUTO  5 

105.6 

0.0920 

0.24840 

2.86 

SAP  5 

0 . 25571 

MUTO  6 

108 . 1 

0 .1290 

0.08014 

18 . 20 

SPA  6 

0.09795 

MUTO  7 

108 . 1 

0.0270 

0.00612 

7.80 

SAP  7 

0.00664 

MUTO  8 

108 . 1 

0 . 1700 

0 . 09886 

16 . 40 

SAP  8 

0.11826 
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TABLE  3.3 
FIVE  BAY  BUILDING 


Case 

2Kw 

3kf 

kf 

Ksl 

kf 

Ks2 

V  .  * 

T  Difference 

VF2 

V  % 

T  Difference 

MUTO  1 

0.5 

0.00339 

0.00195 

0.1659 

0.80 

0.1728 

1.62 

SAP  1 

0 .1673 

0 . 1700 

MUTO  2 

72.08 

0 . 00339 

0.00195 

0 .00635 

16 . 00 

0.00897 

4.60 

SAP  2 

0.00756 

0.00942 

MUTO  3 

0.95 

0.00339 

0.00195 

0.12843 

0.26 

0.12977 

0.20 

0.12877 

0 .12951 

MU TO  4  9.15  0.00339 


0.00195  0.02530  3.98 

0.02635 


0 .02778 
0.02807 


1.03 


CHAPTER  IV 


INTERACTION  FRAME-DIAPHRAGM  FOR  NON SYMMETRICAL  ONE-STOREY 
BUILDINGS 

4.1  INTRODUCTION 

This  chapter  is  devoted  to  the  analysis  of  rectangular, 
nonsymmetr ical  one  storey  buildings  subjected  to  uniformly 
distributed  horizontal  loads  such  as  wind  or  earthquake. 

Typical  buildings  studies  are  shown  in  Figure  4.1. 

Figure  4-la  and  Figure  4-lb  show  buildings  where  the  end  walls 
have  different  stiffnesses.  These  different  stiffnesses 
influence  the  displacement  of  the  building  in  such  a  way  that  it 
is  no  longer  only  translational  but  is  also  rotational. 

When  the  stiffness  of  the  horizontal  diaphragm  is 
infinite  as  shown  in  Figure  4-2a,  the  total  acting  horizontal 
load  can  be  concentrated  at  the  geometric  center  of  the  building, 
which  is  defined  as  the  center  of  gravity  of  the  plan  area  of  the 
building.  In  this  particular  case  the  building  will  rotate 
around  the  center  of  rotation  which  is  defined  as  the  center  of 
gravity  of  the  stiffnesses  of  the  vertical  elements  of  the 
building.  Once  the  center  of  rotation  is  known  the  force 
resisted  by  any  frame  can  be  easily  determined. 

When  the  stiffness  of  the  horizontal  diaphragm  is  not 
infinite  the  total  load  can  not  be  concentrated  at  the  Geometric 
center  and  the  load  acting  in  any  frame  depends  on  the  stiffness 
of  the  frame  relative  to  the  stiffness  of  the  horizontal 
diaphragm  and  also  on  the  location  of  the  frame.  Instead  a 
distributed  load  is  considered  to  act  on  the  horizontal  diaphragm 
which  is  modelled  as  a  beam  with  internal  shearing  force  and 
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Figure  4  "1  Typical  Arrangements 
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Figure  4-2  Torsion  on  Buildings 
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bending  moment  (see  figure  4-2b) . 

When  the  number  of  bays  is  less  than  four  the  beam  model 
yields  results  which  are  readily  calculated.  For  more  than  three 
bays  a  solution  based  on  the  theory  of  beams  on  elastic 
foundations  is  developed. 

4.2  METHOD  OF  ANALYSIS 

As  in  Chapter  III  the  beam  model  has  been  selected  to 
represent  the  horizontal  diaphragm.  To  model  frames  and  walls, 
springs  of  different  constants  are  used. 

For  buildings  having  up  to  three  bays  the  beam  model 
yields  solutions  easily  obtainable.  For  more  than  three  bays  a 
solution  based  on  the  theory  of  beams  on  elastic  foundations  is 
developed. 

4.2.1  BEAM  MODEL 

The  horizontal  diaphragm  will  behave  as  a  deep  plate 
girder,  so  that,  bending  and  shear  must  be  considered.  When  the 
building  has  a  steel  deck  diaphragm  only  shear  properties  are 
considered  since  the  bending  contribution  is  slight.  When  the 
horizontal  diaphragm  is  made  of  concrete,  bending  as  well  as 
shear  properties  must  be  considered. 

4. 2. 1.1  SOLUTION  FOR  A  LIMITED  NUMBER  OF  BAYS 

The  solution  for  the  force  taken  by  an  interior  frame  of 
a  three  bay  building,  shown  in  Figure  4-lb,  will  be  used  to 
explain  how  the  nonsymmetr ical  arrangement  of  vertical  elements 
affects  the  fraction  of  the  total  horizontal  load  taken  by  an 
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in  Appendix  M. 
in  Figure  4-lb  is 

(4.1) 

in  which: 

A  and  B  are  shown  in  Appendix  M. 


interior  frame.  Derivations  are  detailed 
The  force  taken  by  frame  1  shown 


V 

V, 


R 


A 

B 


VR  =  Force  taken  by  the  first  interior  frame 

VT  =  Total  horizontal  force  acting  on  the  building 


As  in  the  symmetrical  case  equation  (4.1)  represents  for 
every  c  value  a  space  surface  whose  intersection  with  the  planes 
y  -  0  and  x  =  0  yield  the  curves  for  the  pure  bending  solution 
and  the  pure  shear  solution,  respectively. 

The  limits  for  the  force  taken  by  a  frame  would  be 


Lim 

c 


0 , 1 8 17x  2  +  x  +  1  .  1 8  2  9  xy  +  y  +  y2 _ 

0.495x2  +  2.907x  +  3.27xy  +  4y  +  3y2  +  1 


_ x  +  y _ 

2 . 7  2  5x  +  3y  +  1 

which  is  the  lower  bound  limit. 

The  intersection  of  this  surface  with  the 
shear  planes  yields  the  curves 


(4.2) 

bending  and 


x 

(2.725x  +  1) 


and 
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X 


(  3y  +  1) 


The  upper  bound  limit  is 


V  0 . 05555(n+l ) 2 

Lim  ~  =  - — — — — =  0.5  (4.3) 

T  0.1111(n+l) 
c+0  4n 

which  is  a  plane  in  the  space  as  shown  in  Figure  3-3. 

Equations  (4.2)  and  (4.3)  are  the  limits  for  any 
possible  set  of  solutions. 

As  in  the  symmetrical  case  the  pure  shear  plane  and  the 
pure  bending  planes  are  the  limits  where  either  only  shear 
properties  or  only  bending  properties  are  represented;  it  is 
necessary  to  have  the  solution  for  the  force  taken  by  an  interior 
frame  for  the  pure  shear  and  the  pure  bending  planes  respectively 
as  shown  in  Figure  4-3.  To  obtain  a  wide  range  of  stiffness  the 
x  axis  is  necessarily  logarithmic.  The  y  axis  represents  the 
percentage  of  the  total  horizontal  force  taken  by  an  interior 
frame.  Parameter  c  has  already  been  defined;  n  for  a  3  bay 
building  is  related  to  c  according  to 

=  (n+1)  Kw/2Kp 


c 


(4.4) 


1.5  SHEAR 
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As  in  the  symmetrical  case  the  curves  that  represent  the 
fraction  of  the  total  horizontal  force  taken  by  an  interior 
frame,  range  between  two  extremes.  These  are  an  infinitely  stiff 
diaphragm  and  an  infinitely  flexible  diaphragm,  relative  to  frame 
stiffness.  . 

When  the  horizontal  diaphragm  is  stiff  relative  to  the 
frame,  the  ratio  Kp/Kg  tends  to  zero  and  for  any  c  the 
displacement  of  the  horizontal  diaphragm  is  comprised  of  a 
lateral  displacement  and  a  rotation,  subsequently  the  force  taken 
by  an  interior  frame  is  proportional  to  its  stiffness  and  to  the 
relative  location  with  respect  to  the  center  of  rotation.  When 
the  ratio  Kp/Kg  tends  to  infinity  the  horizontal  diaphragm  is  not 
able  to  rotate  and  the  distribution  of  loads  to  each  intermediate 
frame  is  a  function  of  the  tributary  area  which  is  identical  to 
the  symmetrical  case.  When  the  ratio  Kp/Kg  is  intermediate  the 
diaphragm  will  be  able  rotate  and  this  will  increase  or  decrease 
the  amount  of  horizontal  force  taken  by  an  interior  frame 
depending  on  its  location. 

4. 2. 1.2  RELATIONSHIP  WITH  THE  SYMMETRICAL  CASE 

The  relationship  between  the  symmetrical  and  the 
nonsymmetrical  case  is  best  shown  when  comparing  the  results  for 
the  force  taken  by  an  interior  frame  for  a  three  bay  building. 

For  a  three  bay  building  the  upper  bound  limit  for  the 
force  taken  by  an  interior  frame  is  the  plane  Z  =  0.5  which  is 
the  same  limit  as  for  symmetrical  buildings.  The  lower  bound 


limit  is  the  surface 


x  +  y _ 

2 . 7  2  5  x  +  3y  +  1 
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z  = 


which  is  the  same  limit  for  symmetrical  buildings. 

When  applying  the  general  solution,  the  solution  for  the 
differential  equation  is  the  same  for  the  symmetrical  and  the 
nonsymmetr ical  cases.  The  only  differences  are  the  boundary 
conditions  imposed  on  the  solution. 


4.2.2  GENERAL  SOLUTION 

The  general  solution  for  a  nonsymme t r ical  arrangement  of 
walls  is  based  on  the  theory  of  a  beam  on  an  elastic  foundation. 

As  shown  in  figure  4-4  the  stiffness  of  the  frames  is 
converted  into  resisting  elements  distributed  per  unit  length. 

Development  of  the  equations  is  explained  in  Appendix  Q. 
Solutions  for  the  differential  equation  are  dependent 
upon  the  factors  3  and  a. 

If  3  <  a 

y^  =  e  1  (  A^  cosf  x  ^ +B^  sinf  x^  )  +  e  1  (  C^cosfx-^+D^  sinfx^  )  +  W/^ff) 


dx  ”^x2 

Y  2  =  e  2 ( A2cosfx2+B2sinfx2 )  +  e  ( C2cosf x2+D2sinf x2 )+W/gF 


If  3  >  a 

( d+f ) x  - ( d+f ) x,  (d-f)x 

y^  =  e  +B^e  +C^e  +^le 


-  ( d-f  )  x 
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Figure  4-4  Nonsymmetrical  Arrangement  of  Walls 
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(d+£)x  ~(d+f)x  (d-f)x  -  ( d-f ) x 

y 2  =  A2 e  +B2e  +C2e  +D2e  +g¥ 


If  3  =  a,  the  solutions  are 


dx,  dx,  -dx  -dx, 

y i  =  Aie  +Bixie  +5  e  +?  xie  v 


dX2  dX2  -dx  2  ~dx 

y0  =  A0e  +B  x  e  Z+C„e  Z+D„x  „e  Z+— 
J  2  2  22  2  22  gF 

in  which 


y-^  =  displacement  of  the  right  part  of  the  building 
Yy  =  displacement  of  the  left  part  of  the  building 

The  boundary  conditions  are 


xi  =  h 


yi  = 


0wi/Kw 


xi  =  1i 


Ms =  0 


x  =  0 


y,  =  Qw/K 


w 


x  _  =  0 


0w2/KW 


X1  =  x2  =  0 


QX  G  A  S  , 
s  s  1 


_f£lZ 

ldx 


O  K 


0S) 

G  A  0 
s  s  2 
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0 


M  =  M 
Si 


y2  =  ^W2/KW 


where 


=  Displacement  of  the  left  part  of  the  building 

y 2  =  Displacement  of  the  right  part  of  the  building 

=  Force  taken  by  the  right  wall 

Kt,  =  Stiffness  of  the  right  wall 

Qw  =  Force  taken  by  the  middle  wall 

Ktt  =  Stiffness  of  the  middle  wall 
W 

0  =  Force  taken  by  the  left  wall 

W2 

=  Stiffness  of  the  left  wall 
W 


CHAPTER  V 


APPLICATION  OF  THE  METHOD 


5.1  Introduction 

This  chapter  is  devoted  to  the  application  of  the 
procedure  discussed  in  the  first  four  chapters  to  two  practical 
examples.  Details  of  calculations  are  shown  in  Appendix  R. 

Both  buildings  analyzed  have  three  bays,  two  end  shear 
walls  and  two  steel  intermediate  frames  as  shown  in  Figure  5,1 
The  horizontal  deck  for  both  cases  is  a  steel  deck;  in  the  first 
case  it  is  filled  with  concrete  whereas  in  the  second  case  there 
is  no  concrete  fill.  To  obtain  the  stiffness  of  the  horizontal 

diaphragm  the  procedure  developed  by  the  U.S.  Department  of  the 

/ 

Army,  Navy  and  Air  Force  (33)  is  applied.  The  horizontal  steel 
deck  characteristics  are  those  given  by  the  manufacturer  (35). 

5.2  THREE  BAY  BUILDING  WITH  CONCRETE  DECK 

The  building  shown  in  Figure  5.1  is  first  designed  for 
vertical  loads,  then  making  use  of  the  beam  model  to  represent 
the  horizontal  diaphragm,  the  building  is  analyzed  for  horizontal 
loads.  To  calculate  the  stiffness  of  the  horizontal  diaphragm 
three  constants  are  assessed: 

The  first  one  considers  the  contribution  of  the  edge 
girders  to  the  overall  moment  of  inertia. 

The  second  one  accounts  for  the  contribution  of  the 
concrete  and  the  steel  deck  to  the  overall  moment  of  inertia. 
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Figure  5-1  One  Storey  Symmetrical  Buildings. 
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Figure 5-2  Steel  Deck  with  Concrete  Fill 
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The  third  coefficient  accounts  for  the  shear  properties 
of  the  horizontal  diaphragm.  All  calculations  are  presented  in 
Appendix  R. 

5.3  THREE  BAY  BUILDING  WITH  STEEL  DECK 

The  same  procedure  as  explained  in  section  5.2  is 
applied,  but  in  this  case  to  calculate  the  stiffness  of  the 
horizontal  diaphragm  only  two  constants  are  assessed: 

The  first  one  considers  the  contribution  of  the  edge 
girders  to  the  overall  moment  of  inertia. 

The  second  one  accounts  for  the  shear  properties  of  the 
horizontal  diaphragm.  The  contribution  of  the  steel  deck  to  the 
overall  moment  of  inertia  is  negligible.  Details  are  given  in 
Apendix  R. 

5.4  DISCUSSION 

The  importance  of  considering  the  stiffness  of  the 
horizontal  disphragm  is  appreciated  when  displacements  and 
horizontal  forces  on  the  vertical  elements  are  calculated. 

Several  variations  of  end  walls  are  considered  in 
Appendix  R  and  the  results  are  summarized  in  Tables  5.1  and  5.2 
respectively  for  concrete  horizontal  diaphragm  and  steel 
horizontal  diaphragm.  Three  different  end  walls  are  considered 
and  are  shown  as  cases  1,  2  and  3  respectively. 

Table  5.1  presents  results  for  a  building  with  the 
horizontal  diaphragm  very  stiff  since  it  is  filled  with 
concrete.  The  first  column  of  Table  5.1  shows  the  stiffness  of 


the  frame  relative  to  the  stiffness  of  the  slab. 


This  value  is 


103 


extremely  small  and  according  to  the  study  performed  in  Chapter 
III  the  horizontal  diaphragm  can  be  assumed  infinitely  stiff. 

The  second  column  of  Table  5.1  shows  the  stiffness  of  the  end 
wall  relative  to  the  stiffness  of  the  interior  frame.  As  the 
stiffness  of  the  end  wall  decreases  the  force  taken  by  an 
interior  frame  increases. 

In  the  second  case,  summarized  in  Table  5.2,  the 
horizontal  diaphragm  is  a  steel  deck  with  no  concrete;  this  makes 
the  horizontal  diaphragm  very  flexible.  The  assumption  that  such 
a  horizontal  diaphragm  is  infinitely  stiff  is  no  longer  valid  as 
it  is  seen  when  comparing  results  of  the  force  taken  by  an 
interior  frame.  When  the  stiffness  of  the  horizontal  diaphragm 
is  calculated  the  force  taken  by  an  interior  frame  may  be  greatly 
increased.  This  increase  is  best  shown  when  comparing  the  first 
two  lines  of  Table  5.2.  There  is  an  increase  of  the  force  taken 
by  an  interior  frame  of  more  than  two  hundred  per  cent.  When  the 
stiffness  of  the  wall  relative  to  the  stiffness  of  the  frame  is 
1.0  the  assumption  that  the  horizontal  diaphragm  is  infinitely 
stiff  is  not  bad,  since  according  to  study  performed  in  Chapter 
III,  this  case  approximates  the  upper  bound  solution  for  the 
whole  set  of  curves. 


Forces  in  frame  for  concrete  Diaphragm 
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CHAPTER  VI 


SUMMARY  CONCLUSIONS  AND  RECOMMENDATIONS 


6 . 1  SUMMARY 

An  investigation  of  the  interaction  between  frame  and 
diaphragm  for  horizontal  loads  acting  on  one  storey  rectangular 
buildings  has  been  presented. 

Chapter  II  includes  a  critical  review  of  different 
approaches  to  calculate  the  strength  and  stiffness  of  steel  deck 
diaphragms.  A  new  approach  to  calculate  the  stiffness  of  steel 
diaphragms  has  been  developed  and  presented.  A  parametric  study 
of  the  different  factors  that  govern  the  strength  and  stiffness 
of  a  steel  deck  diaphragm  has  been  presented. 

In  Chapter  III  a  particular  solution  for  the  force  taken 
by  an  interior  frame  in  a  symmetrical  one  storey  rectangular 
building  has  been  developed.  This  solution  may  be  applied  to 
buildings  that  have  up  to  five  bays.  A  general  solution  based  on 
the  theory  of  a  beam  on  an  elastic  foundation  has  been  presented 
as  well.  This  solution  applies  for  any  number  of  bays. 

In  Chapter  IV  a  particular  solution  for  the  force  taken 
by  an  interior  frame  for  rectangular  nonsymmetr ical  one  storey 
buildings  under  horizontal  loads  has  been  presented.  This 
solution  may  be  applied  to  buildings  having  up  to  three  bays.  A 
general  solution  based  on  the  theory  of  a  beam  on  an  elastic 
foundation  has  been  developed  which  provides  the  solution  for  the 
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force  taken  by  an  interior  frame  for  buildings  having  any  number 
of  bays. 

Chapter  V  presents  two  practical  examples  of  the  theory 
developed  in  Chapter  III  with  variations  in  stiffness  of  end 
walls  and  interior  frames.  Effects  of  variation  in  diaphragm 
stiffness  relative  to  end  wall  and  frame  stiffness  are  discussed. 

6.2  CONCLUSIONS 

1.  The  most  accurate  analytical  method  to  assess  the  stiffness 
of  a  steel  deck  diaphragm  under  in-plane  loads  is  the  method 
developed  by  Davies  (11)  in  Britain. 

2.  The  most  important  factors  that  affect  the  strength  and 
stiffness  of  a  steel  deck  diaphragm  are:  panel  length, 
purlin  spacing,  number  of  seam  connections  and  number  of 
shear  connectors. 

3.  The  parameter  Kp/Ks  (stiffness  of  the  frame/stiffness  of  the 
horizontal  diaphragm)  is  the  first  most  important  parameter 
that  affects  the  amount  of  horizontal  load  taken  by  an 
internal  frame 

4.  The  parameter  c  =  EKri/EK„  (summation  of  the  stiffness  of 
the  outerwalls/summation  of  the  stiffness  of  the  internal 
frames)  is  the  second  most  important  parameter  that  affects 
the  amount  of  horizontal  load  taken  by  an  internal  frame. 

5.  The  curve  that  represents  the  percentage  of  horizontal  force 
taken  by  an  internal  frame  is  dependent  upon  the  bending  and 
shearing  properties  of  the  horizontal  diaphragm.  The 
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solution  lies  in  a  plane  which  is  bounded  by  the  pure  shear 
plane  and  the  pure  bending  plane. 

6.  The  solution  in  the  pure  shear  plane  is  a  valid  solution  for 
buildings  having  a  steel  deck  horizontal  diaphragm  when 
using  the  beam  model. 

7.  The  solution  for  the  force  taken  by  an  internal  frame  is 
always  bounded  by  a  lower  bound  space  surface  when  c  +  °°  and 
an  upper  bound  space  surface  when  c+0.  This  two  limits  are 
dependent  upon  the  number  of  internal  frames. 

6.3  LIMITATIONS  OF  THE  METHOD  AND  RECOMMENDATIONS 

The  solution  for  the  stiffness  of  a  steel  deck  diaphragm 
under  in-plane  loads  in  Chapter  II  is  an  elastic  one;  non  linear 
behaviour  of  connections  and  redistribution  of  stresses  takes 
place  before  the  failure  of  the  diaphragm  occurs.  It  is 
recommended  that  nonlinear  behaviour  be  incorporated  in  the 
analysis  of  the  diaphragm. 

The  solution  discussed  in  Chapters  III  and  IV  is  based 
in  the  Beam  model  so  that  the  ratio  between  length  and  width  of 
the  building  must  be  larger  than  2.0.  It  is  also  a  linear 
solution  since  both  shear  and  bending  deflections  are 
superimposed.  Therefore,  there  is  a  need  to  study  this  matter 
more  thoroughly.  A  finite  element  analysis  that  incorporates  the 
Geometrical  and  Physical  properties  of  the  horizontal  diaphragm 
may  be  applied.  The  study  may  be  widened  for  buildings  of  more 
than  one  storey. 
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APPENDIX  A 


PLATE  BUCKLING  FAILURE 


Bergmann,  V.S.  and  Reissner 's  Method 

Bergmann  and  Reissner  (4)  calculated  buckling  loads  of 
rectangular  corrugated  orthotropic  plates  with  the  aid  of  Eg. 

A.  1 . 
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Buckling  load 
panel  width 

Bending  stiffness  per  unit  length  of  diaphragm  in  the 
x  direction 

Bending  stiffness  per  unit  length  of  diaphragm  in  the 
y  direction 

Multiplier  dependent  upon  0  and  <j> 
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(A. 3) 


where 

a  =  Panel  length 

The  dependence  of  X  on  0  and  <j)  is  given  in  the  form  of 
curves  in  Reference  30. 
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Hlavacek's  Method 
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Hlavacek  (21)  calculated 


Equation  A. 4. 
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buckling  loads  with  the  aid  of 


(A. 4) 


Easley  -  McFarland's  Method 

Easley  and  McFarland  (14)  calculated  buckling  loads  using 
Equation  A. 5. 
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(A. 5) 


Easley's  Method 

Easley  (15)  calculated  buckling  loads  with  the  aid  of  Eq. 

A.  6 . 
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cr 


(A.  6) 


Where  £  is  a  parameter  that  represents  boundary 
conditions.  The  lower  and  upper  boundaries  for  £  are  1.0  and 
1.9,  but  the  true  variation  of  £  with  diaphragm  end  restraint  is 


unknown. 


APPENDIX  B 


BRYAN'S  METHOD 


Stiffness 

Bryan  (7)  states  that  the  overall  shear  flexibility 
shear  panel  can  be  evaluated  by  considering  the  effect  of 
separate  components.  These  components  can  be  summarized 
follows: 

Sheet  Distortion 


0.144  a  d 


0 


fl  Ko 


1 . 1 


3  3 
Et  b 


(B.l) 


where 
a  = 


Panel  Length 
Pitch  of  corrugations 
correction  factor 
sheeting  constant 
Modulus  of  Elasticity 
sheet  thickness 
Panel  width 


(Reference  1,  p.  110) 
(Reference  7,  p.  109) 


Shear  strain  in  the  sheet 


2  a  f  ( 1+v) ( l+2h/d) 
C1 . 2~  btE 


(B.2) 


of  a 


as 
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where 

a  =  Panel  Length 

f2  =  correction  factor  (Reference  7 ,  p.  110) 

v  =  Poisson's  ratio 

h  =  height  of  a  corrugation 

d  =  Pitch  of  corrugations 

b  =  Panel  width 

t  =  sheet  thickness 

E  =  Modulus  of  Elasticity; 

Axial  strain  in  the  Purlins 

(B.3) 

(Reference  7 ,  p.  110) 


1.3 

where 


o  3, 
2  a  f . 


3b  AqE 


f-^  =  correction  factor 


Aq  =  cross  area 


Sheet  -  Purlin  Fasteners  (Sheeting  fixed  on  all  four  edges) 


2  a 

1  = 


(B.4) 


where 

q  =  pitch  of  sheet-purlin  fasteners 

Sp  =  flexibility  per  sheet-purlin  fastener 

f  3  = 


correction  factor 
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Sheet-Purlin  Fasteners  (sheeting  fixed  to  the  purlins) 
2S  q 

C2.1  =  (6/P+2  +  a2f3/b2)  (B.5) 


Seam  Fasteners  (sheeting  fixed  on  all  four  edges) 


C2.2  =  (nsh"1)  Ss/ns 


(  B.  6) 


where 

ng^  =  number  of  sub-panels 

Sg  =  flexibility  of  a  seam  fastener 

ng  =  number  of  seam  fasteners 


Seam  Fasteners  (sheeting  fixed  on  two  edges) 


2 . 2 


n  S  /n 
sh  s  s 


(  B  c  7 ) 


Sheet-Connector  Fasteners 

c0  =  2  S  /n 
2 . 3  sc  sc 

where 


(  B0  8) 


SC(,  =  flexibility  of  sheet-connector  fasteners 
n  =  number  of  sheet-connector  fasteners 


APPENDIX  C 


DAVIES'  METHOD 


Stiffness 

In  addition  to  the  components  already  obtained  by  Bryan, 
Davies  (10)  modifies  the  following  coefficients 


Sheet-purlin  fastener 


c2.1 


2  a  S 

P 


q  f3/b2 * * S 


(Cel) 


where 

a  =  panel  length 

S  =  flexibility  per  sheet-purlin  fastener 

■b'' 

q  =  pitch  of  sheet-purlin  fastener 
f^  =  correction  factor 
b  =  panel  width 


Seam  fastener 


2  S  S  (n  ,  -1) 

s  p  sh _ 

2 . 2~  2  n _  S^  +  g,  (p  +  2)  S, 

S  P  -L 


where 

Sc  =  flexibility  of  a  seam  fastener 
ng^  =  number  of  sub-panels 


(C.2) 
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n3  =  number  of  seam  fasteners 
p  =  number  of  intermediate  purlins 


n  /2 

lp 

i  =  l 


(C.3) 


where 

nQ  =  total  number  of  sheet-purlin  fasteners 
per  sub-panel,  n  even 

AT 


P 


n  odd 

hr 


Shear  connector  (direct  shear  transfer) 


2  S 


sc 


'2.3  n 


sc 


(C.5) 


APPENDIX  D 


EASLEY'S  METHOD 


Strength 
End  Fastener 


Easley 

fastener  is 
„  _  Q 


epli  nsh  n 


Fa i lure 
(13)  assumes 
constant. 

e 


that  the  horizontal  force  taken  by  each 

(Del) 


The  vertical  component  is 


F 

epvi 


2x .  n  ,  Qb 
1  sh 

a2  8, 


where 


(  D.  2) 


31  =  ns  + 


41  +  2  P  I 

e_ g 


''n  J 


sh 


where 


le  ■  r  x 

i  =  ’ 

Ig  =  Ip 

=1 


where 

a  =  Length  of  the  panel 

b  =  width  of  the  panel  (parallel  to  the  direction  of  the 

corrugations ) 


122 


123 


ep 


h 


epvi 


R. 

i 


n< 

n{ 

nj 

n 

P 

Q 

x. 


sh 


=  Horizontal  force  component  in  the  fasteners 
=  Vertical  Force  component 
=  Resultant  Force 

=  number  of  end  fasteners  per  sub-panel 
=  number  of  purlin  fasteners 
=  number  of  seam  fasteners 
=  number  of  sub-panels  per  panel 
=  number  of  intermediate  purlins  per  panel 
=  Applied  force  as  shown  in  Figure  2.6a 
=  Distance  from  every  end  fastener  to  the  center  of 
rotation  of  the  sub-panel 

=  Distance  from  every  purlin  fastener  to  the  center  of 


rotation  of  the  sub-panel 

( Fe  ,  2  +  Fe  .  )  ^ 
eph  pvi 


(  D.  3 ) 


Seam  Fastener  Failure 


From  Figure  2-5 

F  = 
s  a  31 

where 


(  D.  4 ) 


F  =  side  fastener  force 
s 


The  lowest  value  from  equations  D. 3  and  D.  4  gives  the 
strength  of  the  steel  deck  diaphragm  under  in  plane  loads. 


STIFFNESS 


Easley's  equation  for  the  stiffness  of  a  diaphragm  is  given 
following  the  AISI  (1)  recommendations  for  shear  diaphragms. 


K  = 


Gtl 


+ 


A-,  b  Gk 
1  sh 

0392  F  a  $ . 

u  1 


(  D.  5 ) 
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In  which 


d 

length  of  one  corrugation  measured  along  the 

corrugation 

Fu 

ultimate  strength  of  the  fastener  connection 

G 

shear  modulus  of  elasticity  of  the  panel  material 

l 

projection  of  the  length,  d,  on  the  plane  of  the  panel 

nsh 

number  of  sub-panels 

t  = 

material  thickness 

X1 

constant 

According 

to  Easley  the  parameters  Fu  and  are  obtained  by 

testing. 

APPENDIX  E 


Fazio's  Et  Al.  Approach 

Strength 

According  to  this  method  (17)  the  steel  deck  may  have  three 
modes  of  failures. 

Failure  at  Seam 


0 


K  ( n  F  +PF  +  2  F  ) 
1  s  su_ pu_ eu 


ult 


n  ,  X 
sh  o 


-  i)( 


-2®  +  +  -1) 
s  s  s  J 

s  d  e 


where 


K1  =  s 


ns 


n  ,  X 
sh  o 


-  i)  +  (g2e  "  gip}  +  ^  (go^  "  gi-} 

e 


S  '  ^  2p 
P 


lp 


and 


ns 


X  = 
o 


2gle  l  P  „ 

S  S  glp 

e  p  c 


n  2  2  P  n 

sc  +  ns  +  ne  +  _ p 


sc 


S 


and 


le 


n 


sh 


n 


r 

i  =  l 


X 


n 


g 


2e  aX 


n 


o  1  =  1 


l  ,  Xi 


n  ,  n 
sh  vp 

9ip  '  a  p=i 


n 


sh 


2p  a  X 


n 


o  p  =  l 
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N 


Failure  At  Side 


°ult 

Ko  (n  F  +  P  F  +  2  F  ) 

■n-  n  °  ^ - —  <E-2) 

sc  p  e 
Failure  At  End 


Qult 

F 

=  eu  (F  SI 

(F  2  +  F  2)  1/2 

eph  epv 

where 


a  = 

length  of  diaphragm  in  direction  perpendicular  to 

corrugations 

Feph~ 

horizontal  constant  force  as  shown  in  Figure  E.l 

F  = 
epv 

maximum  vertical  force  at  the  end  fastener 

F  = 
r  eu 

strength  of  end  fastener 

c 

it 

strength  of  sheet  to  purlin  fastener 

Fsc 

strength  of  side  fastener 

Fsu 

strength  of  seam  fastener 

ne  = 

number  of  end  fasteners  per  panel  end 

np 

number  of  purlin  fasteners  per  panel 

P 

number  of  purlins 

ns 

number  of  seam  fasteners  per  seam  line 

nsc 

number  of  side  fasteners  per  side 

nsh  = 

number  of  sub-panel  per  panel 

se  = 

flexibility  of  end  fastener 

SP  = 

flexibility  of  purlin  fastener 

Ss  = 

flexibility  of  seam  fasteners 

Ssc 

flexibility  of  side  fasteners 

APPENDIX  F 


EMPIRICAL  METHOD  FOR  STRENGTH  AND  STIFFNESS 
U.S.  DEPARTMENTS  OF  THE  ARMY ,  NAVY  AND  AIR  FORCE 

CONCRETE  DIAPHRAGMS 

This  type  of  diaphragm  is  made  of  a  steel  deck  with  a 
superimposed  fill  of  concrete  having  a  minimum  fic  of  2,500 
p.s.i.  at  28  days  and  a  minimum  weight  of  90  p.s.f. 

Galvanized  steel  deck  is  used.  The  method  to  calculate 
shear  stiffness  and  the  working  shear  is  given  empirically  by  the 
Department  of  the  Army,  Air  Force  and  Navy  (31). 

Two  types  of  steel  deck  diagphrams  with  concrete  are 
categorized : 

Type  1.  If  the  diaphragm  is  loaded  without  shear  stresses 
passing  through  the  deck  or  its  attachments  as  shown  in  Figure 

F .  1 .  a . 

Type  2.  If  the  shear  passes  through  the  deck  or  its  attachments 
as  shown  in  Figure  F.l.b. 

Stiffness 

The  factor  F  is  introduced  to  calculate  the  stiffness  of  a 
horizontal  diaphragm  and  is  defined  as  the  average  deflection  in 
microinches  of  the  diaphragm  web,  per  foot  of  span,  stressed  with 
a  shear  of  one  pound  per  foot. 

A  •  106 
„  w 
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8 


Figure  F-l  Steel  Deck  with  Concrete  Fill 
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(a)  Type  A 


Detail  A 


(b)  Type  B 


Figure  F-2  Steel  Deck  Diaphragms 
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where 

]_1  =Distance  in  ft.  between  vertical  resisting  elements  and  the 
point  to  which  the  deflection  is  to  be  determined. 

Qave  ^Average  shear  in  diaphragm  in  pounds  per  ft.  over  length 
L-^  .  r 

Aw  =Web  deflection  in  inches. 

TYPE  1 


Strength 

The  criteria  to  calculate  the  strength  of  a  concrete 
diaphragm  is  given  by  ACI  318-~71  (2)  but  limited  to  the  value  of 


the  equation 

6.0  f'c  .  t 


*ud 


(  F .  1 ) 


1  + 


2  2 

( P+1)  t 


where 

qU(j  =shear  capacity  in  pounds  per  foot 
f'c  =compress ive  strength  of  concrete 
b  =width  of  the  building 
P  =  Number  of  interior  purlins 
t  =Thickness  of  slab  in  inches 

The  flexibility  factor  F  is  calculated  using  the  formula 


F  =  - 10  _____  (  F .  2  ) 

8 . 5  t  y1 *  *  /f'c 

where 

t  =Thickness  of  the  slab  in  inches 
y  =Weight  of  the  concrete  in  pounds  per  cubic  foot, 
f'c  =The  compressive  strength  of  the  concrete  at  28  days  in 
pounds  per  square  inch. 
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TYPE  2 


Strength 


qD  +  ql  =  q6 


where 


qi  = 


92  S  (t1+t2’)  K  (P+1) 
— -  —  — 


in  which  K2  =  1,000 
q6  =  q6’+  qfi" 


where 

V  = 

and 


tf  y  1,5  /f8C 

200 


V  =  z/ 
where 


ka 


Xl(tl+t2)  nsh 


n 


sh 


(F.3) 


(F.4) 


(  F .  5  ) 


(F.6) 


(F.7) 


qD  -  Working  shear  in  pounds  per  foot 
a  =  Length  of  the  diaphragm 

=  Distance  in  feet  between  outermost  puddle  welds  attaching  a 
deck  unit  to  the  supporting  framing  member. 

S  =  Section  modulus  in  feet  of  puddle  weld  group  at  supports 
t^  -  Thickness  of  flat  sheet  elements  in  inches 
t2*=  Effective  thickness  of  fluted  elements  in  inches 
tf  =  Thickness  of  fill  over  top  of  deck  in  inches. 

The  Flexibility  factor  F,  is  defined  using  the  Formula 


F 


20 


STEEL  DECK  DIAPHRAGMS 


(F.8) 


Two  types  of  steel  deck  diaphragms  are  categorized:  Type  A 
and  Type  B  as  shown  in  Figure  F-2.  Definition  of  the  symbols 
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used  are 


P 


w 


■5 


F1  'F2*  * ‘ 


as  follows: 

=  Length  of  the  diaphragm 

Averge  spacing  of  profile  channel  closures,  in  feet. 

=  Center  to  center  spacing  of  seam  welds  in  feet. 

=  Spacing  of  margina  Iwelds  in  feet. 

=  1  for  galvanized  decks;  0.65  for  painted  decks. 

=  1  for  button-punched  seams;  40t  ^^l'w  for  welded 

s 

seams 

=  1  for  button-punch  seams;  150  t  1'  for  welded 

o  w 

seams . 

6 ( P+1 ) 

=  1  for  button-punched  seams ;  for  welded 

seams . 

=  1.2  for  continuous  angle  closure;  1  for  continuous 

zee  closure;  ^  for  profile  channel  closure. 

ci 

P 

=  Distance  in  feet  between  outermost  puddle  welds 
attaching  a  deck  unit  to  the  supporting  framing 
member. 

=  Components  contributing  to  the  flexibility  factor  F 


f  ' 


c 


h 


I 


D 


(F  =  ZF  ) . 
n 

Compressive  strength  of  fill  concrete  at  28  days  in 
pounds  per  square  inch. 

Height  of  fluted  elements  in  inches 

Gross  moment  of  inertia  of  deck  unit  about  vertical 
centerline  axid  through  unit  in  inches  to  the  fourth 
power. 

Gross  moment  of  inertia  of  deck  unit  above  the 
horizontal  neutral  axis  of  the  deck  cross-section 
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l2 

!e 

LR 

Lv 

^w 

!’w 


n 


n 


s 


qD 


ql '  q2  *  *  * 


qave 


R 

S 


t 


1 


t 


2 


t' 


2 


per  foot  of  width  in  inches  to  the  fourth  power. 
Distance  in  feet  between  vertical  resisting  element 
(such  as  shear  wall)  and  the  point  to  which  the 
deflection  is  to  be  determined. 

Average  length  of  each  deck  unit  in  feet. 

Length  of  edge  lip  deck  panel  in  inches 
Distance  in  feet  between  shear  transfer  elements 
Vertical  load  span  of  deck  units  in  feet. 

Minimum  length  in  inches  of  seam  weld. 

Effective  length  in  inches  of  seam  weld. 

Average  number  of  vertical  deck  elements  per  foot 
which  are  laterally  restrained  at  the  bottom  by 
puddle  welds. 

Number  of  seam  attachments  in  span  b  along  a  seam 
Working  shear  in  pounds  per  foot.  The  one-thir 
increase  usually  permitted  on  working  stresses  is 
not  applilcable  to  this  vaue. 

Components  or  limiting  vaues  of  working  shear  in 
pounds  per  foot. 

i 

Average  shear  in  diaphragm  over  length  in  pounds 

per  foot, 
b 

( P+1 ) l2 

Section  modulus  in  feet  of  puddle  weld  group  at 
supports.  (Each  weld  assumed  as  unit  area.) 
Thickness  of  flat  sheet  elements  in  inches. 
Thickness  of  fluted  element  in  inches. 

Effective  thickness  of  fluted  elements  in  inches. 


13^ 


ts 

Y 

TYPE  A 
Strength 


Thickness  of  closure  element  in  inches. 

Thickness  of  fill  over  top  of  deck  in  inches. 
Thickness  in  inches  of  deck  steet  at  seams. 

Unit  weight  of  fill  concrete  in  pounds  per  cubic 
foot . 


qn  =  (q-,  +  q0)  q3/q2«  where  g^/g0  <  c 


LD 
and 


I  x  10 
x 


3  2 


2 


( p  +  1 ) 


or 


(F.9) 


10 


1 . 5/ 


P+ 1 { F 1  + 


F2  + 


f3  +  4 

12 


if  1  <0.5” 

e 


92  S  (t  +  t2')  K2(P+1)  n  h 
ab 


where 


K2=- 


1,000 


0.5 


1  +  S 


(tl  +  t2) 


+  100  m 


0.5 


t2/li  r- 

2  h  *-t 


1  +  fc2 


.  0.5 

n  a  t  crm 

s  s _ 2  c  f 500 

qo  ^  n  1  q ]_  l  t 


0.5  (P+1) 


0.5 


n 


sh 


D 


b  x1S 


(tl  +  fc2)2 


2 
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3600  t  n  C-.  ( p+1 ) 
S  S  j 

k 


The  flexibility  F  is  defined  as 

F  =  Fi  +  F2  +  F3 

where 

F  =  1 
1  12  (tl  +  t2) 

a  b  C4  f500  0.5  (P+1)  >  ql 

^  2  2  l  t  +  2^rr+rr 

( P+1 )  n  ,  160  1 D  b  X,  S  (t,  +  t0)  ql  q2 

sh  112 


F3  = 


R  (P+1) 


( t  ,  + 


2  2  3 
12.5  in  CZ 

- L_f 

h  ’ 


TYPE  B 
Strength 


,  or  q,_  whichever  is  the  lesser  but  limited  to  1050  lb/ft 


0.6  t2  n  11(P+1) 
s  s  w 

b 
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q5  7 


2 

C,t 
5  c 


10 


2h 


0.5 


Flexibility 
F  =  Fi  +  F4  +  F5 


where 


F1  12  (t1  +  t2) 


_  _  3500 

4  =  ~q7" 


„  _  20000 
= 


APPENDIX  G 


STIFFNESS  OF  THE  SLAB 


The  horizontal  diaphragm  is  shown  in  Figure  G-l.  It 
considered  as  a  simple  supported  beam  between  shear  walls. 
Considering  bending  plus  shear  deflections 


A  = 


r 


ii  (1-V  ^ 


Wx' 


•)  dx 


Es  Is 


X 


,  r  /  T  wWlx  Wx  .  , 

1  61(l-x)(— J-  -  -y)  dx 


Es  Is 


+ 


Wl 

(  1-61)  (y  -  W) 
Gs\ 


dx 


Wl 

5i(—  -  Wx)  dx 


GsAs 


The  integration  of  these  quantities  yield 


A 


Wl 


24ESIS 


(  6. 


■2  <S^  +  51  ) 


Wl 


2GsAs 


Therefore  the  stiffness  of  the  slab  is 


K  = 


Wl 

A 


13(  S14-2«13  +  «1  )  1(«1-612) 


24  Esxs 


2  GsAs 


is 
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i 


J  Z3HEZ T'"~T 


W 

r  i  r  i  r 


(b)  Shear 


Figure  G-l  Stiffness  of  the  slab 


Length  of  the  slab  between  end  shear  walls 


Shear  area  of  the  horizontal  diaphragm 
Young's  Modulus  of  the  horizontal  diaphragm 
Shearing  Modulus  of  the  horizontal  diaphragm 
Moment  of  Inertia  of  the  horizontal  diaphragm 
Uniformly  distirbuted  horizontal  load 
Location  coefficient 

Deflection  of  the  horizontal  diaphragm 


APPENDIX  H 


SYMMETRICAL  TWO  BAY  BUILDING 
TRUSS  SOLUTION 

The  structure  is  shown  in  Figure  H-l 
The  matrix  solution  yields? 


4  -  K  U,  =  K  U,  -  K  U„  +  O.U, 
2  wl  si  s2  1 


P1  -  KFU2  “  -KsUl  +  KsU2  -  KsUl 


-  K  U,  =  O.U,  -  K  U-  +  K  U. 
2  wl  1  s2  si 


From  which 


2P  (K  +  V?  K  ) 
U  =  _  3 _ J—H— 

U2  2  K  K  +  K  K  +  K  K 
s  w  F  S  F  w 


V  = 


V 


_F 

2P 


KFU2 

TTT7k 


K, 


^  4k 
F  s 


1  kf 
1  +  c  +  2K 


KF 

Replacing  —  for  x 
K 

s 


VF 

and  ^  for  Z 


Z  = 


i  + 1 

c  4 


1  +  ±  +  2 
c  2 
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P 


i 


Figure  H-l  Two  Boy  Building 
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which 

coincides  with  the  pure  shear  solution  given  by  the  beam 

model 

when  replacing 

K  ,  = 
sh 

8  GsSs 

1 

for 


K  ,  = 

sh^ 

2  GsAs 

1 

APPENDIX  I 


Two  Bay  Building 

Referring  to  Figure  1-1,  the  deflection  at  the  center  of  the 
span  for  a  uniformly  distributed  load  iss 


A  5  Wl4  ,  Wl2  Wl  W1 

A,  =  zr— —  +  — —  +  —  + 


1  384  ESIS  8G  A  4Kw  4  nKw 


The  Deflection  at  the  center  of  the  building  for  a  point  load  is 


2  "Vs  '  4GSAS  '  4Kw  4"  Kw 


Equating  the  deflections 


V 

A1  ”  A2  =  K 


R 


Thus 


5  Wl 


384  EI0 
S  S 


+ 


Wl 


8GsAs  4K« 


Wl  Wl 


V 


V 


4  n  K 


w 


48VS  4GgAs 


V 


R 


V 


R 


V 


R 


4K  4 n  K  K„ 

w  w  F 
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co 

O 


in 
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S£*  0 


Sfi  *0 


oe*o 


si  *o 


oo*o 


Figure  t-l  Two  Bays  Torsion 
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Therefore 


4  2 

Wl  ,  Wl 

T  — - - — 


VR  “ 


384  ESIS  8G  A 

O  u 


W1  ,  W1 

T  ;■  ■■  -  + 


4K  4nK 
w  w 


48EsIs  .  4=sAS 


1  11 

— +  - -  + 

4K  4nK  K_ 
w  w  F 


If  the  Bending  stiffness  of  the  horizontal  diaphragm  is  defined 


as 


KB  =  5 


384  ESIS 


and  the  shear  stiffness  of  the  horizontal  diaphragm  defined  as 


K 


8GgAs 


sh 


The  solution  is 


_i  +  _L_  +  JL  +  _J_ 

K„  K  .  4 K  4nK 

R  _  B  sh _ w _ w_ _ 

VT  l-6  ,  2  +  1  +  _J_  +  _I 

Kn  K  ,  4 K  4nK  K 

B  sh  w  w  F 


Replacing 


x 


K 

c  =  —  (  n+1 ) 


Thus 


x  +  y 


,  (n+l ) 2 
+  ““ 4nc~ 


1 . 6x  +  2y  + 


(n+l)2, 

4nc 


1 


If  n  =  1  the  symmetrical  case  is  obtained 

Is  .  x  +  y  +  ^ 

VT  1.6x  +  2y  +  -  +  1 

o 


v 

R  x  +  y 

Lim  VT  1.6x  +  2y  +  1 
c  00 


L. 


im 


c  o 


1 


These  are  the  limiting  surfaces  for  the  complete  set  of 


surfaces 
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The  intersection  of  the  above  written  expression  wi 
plane  x  =  o  is  represented  by 


1  + 


Z  = 


2y  +  —  +  1 

c 


where  2 


h 

V„ 


which  is  a  pure  shear  solution 


L.  Z 
im 

y  +  o 


1 

c 


I +  1 


im  Z 

y  +  ® 


=  0.5 


The  intersection  with  y  =  o  is 


Z  = 


x  +  - 
c 


1 . 6x  +  —  +  1 
c 


which  is  a  pure  bending  solution 


L.  Z 
im 


_1 

c 


+  1 


th  the 


x  -»■  o 
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This  value  coincides  with  the  value  obtained  for  pure  shear 


L.  Z  =  0.625 
im 

X  ->■  oo 


The  combined  solution  is  represented  by  a  set  of  surfaces 
which  may  be  intersected  by  planes  representing  selected 
combinations  of  shear  and  bending  stiffnesses.  If  the  plane 
x  =  y  is  selected 


Z  = 


2x  +  — 
c 


3 . 6x  + 


i  + 1 
c 


L.  Z 
im 


x  -*o 


L.  Z  =  0.5555 
im 


x  °° 
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If  the  plane  y  =  2x  is  selected 


Z 


5 


3x 


.  6x 


1 


L .  Z 
lm 

x  -»■  o 


1 


c 


c 


L.  Z  =  0.5357 
im 

x  00 


APPENDIX  J 


THREE 

BAYS  SYMMETRICAL 

BUILDING 

Based 

in 

Figure 

J.l  the 

deflections 

A ^  and 

> 

ii 

VR 

i3 

486 

Vs 

A  — 

7 

vi3 

R 

A2 

486 

Vs 

Thus , 

the 

total 

bending 

deflection 

is 

A  —  • 

15 

V3 

ar 

486 

EsIs 

The  shear 

deflection  is 

• 

• 

A  ,  = 
sh 

VR 

1 

3Gs 

5s 

Then 

replacing 

the  end 

supports  by 

yielding 

deflection  is 

ii 

rH 

< 

15 

3 

V  , 

v3  , 

VR 

486 

EsIs 

3GsSs 

Kw 

are  respectively 


supports,  the  total 
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Figure  J-l  Three  Bays  Symmetrical  Building 
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The  Deflection  in  the  opposite  direction  is 


= 


4.35  Wl4  Wl2 


W1 


3,4  Vs  «s»s 


2K 


w 


Equating  the  deflections,  finally 


4.35  Wl 


Wl 


Wl 


v_  = 


9GSAS 


2K 


W 


15  1' 


486  E„I„ 
S  S 


_L__  +  _i  + 

3GgAg  Kw  KF 


( J-l) 


Defining  the  bending  stiffness  as 


K  _  384  Vs 
4.35  13 


and  the  shearing  stiffness  as 


K 


9GsAs 


sh 


Replacing  in  J.l 


K 


+ 


B 


K 


sh 


2K 

w 


2.725 


K 


B 


K 


+ 


sh 


K 


w 


Defining 


c 


(n+1 ) 


if  n 


1 


x 


=  y 


x  +  y 


+ 


1 

2c 


2 „ 725x  +  3y  + 


If  a  solution  based  in  bending  is  sought 


x 


+ 


2c 


2.725x 
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This  curve  lies  in  the  Bending  Plane. 
The  solution  based  in  shear  yields 


1 


L imi ts 


Bending  solution 


Lim 
X  -►  o 


V 

V. 


R 


2c 


1  +  - 
c 


VT? 

L.  ™  =  0.367 
im  VT 

x  00 


Shearing  solution 


VR  2c 


im  VT  i+l 
y  +  o  c 
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L  . 


irn 


y  - 


=  0.333 
VT 


06 


and 


o 


2 


x 


7  25x  + 


1 


0 . 5 


APPENDIX  K 


FOUR  BAYS  SYMMETRICAL  BUILDING 

Solution  for  pure  shear  and  pure  bending  are  obtained  if  the 
following  procedure  is  used  (See  Figure  K-l). 

Equilibrium  yields: 

2  V  +  2  VA  +  V  =  2  F  +  2  FA  +  Fd 
w  B  w  B 

Diaphragm  shear  first  span 


V 

w 


FA  + 


JB 

2 


Diaphragm  distortion 


V 

VA  w 


Thus 


FA  + 


VA 


V 

w 

K~ 

w 


(K.  1) 
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Figure  K-! 


Three  and  Four  Bay  Buildings. 
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Taking  moments  with  respect  to  the  wall 


VA  + 


Diaphragm  shear  second  span 


V 


2 


Diaphragm  distortion 


ad2  =  ab 


V 


D. 


V 


B 


K, 


K, 


V. 


K, 


2K, 


(  F 


B 


V 


From  (K-l) 


h 

Kt 


V 

K 


w 


B 


V, 


V 


B 


w 


K, 


2K, 


K, 


2K, 


A  system  of  three  simultaneous  equations  is  obtained 


+  V, 


B  2K, 
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•V  —  +  V  [~  +  — — .1 
A  Kp  B  LKp  2Ks  J 


B 


2K, 


VA  +  VB  V2  +  Vw 


Wl 

2 


Finally  a  solution  is  obtained 

,  K_  K  K_ 

v  "3c  '1  +  2F^  +  FA  k7U  +  2k7 

A  _  1 _ bl _ ^1 

Wl 


^  +  FB  2K_  [1  3c ^ 


1  +  _i  +  .  ,  2kf  ,  kf2 
c  Kc  3c  K  2K  2 

bl  bl  bl 


If  a  Bending  solution  is  sought 
Wl 

FA  =  1.14  3  — r 


Wl 

Fb  =  0.929  ™ 


Therefore 


1a 

Wl 


1  K  K 

™  +  0.128  —  +  0.4019  - 

3c  K  C  K 

bl  bl 


p  K  2 

F  +  0.143 


i  +  Kp 

c  K„  3c 


+  1  + 


2K  K  2 

F  +  0.5  F 


K, 


1 


Kq  2 
bl 


Kc  2 
bl 


in  which 


c  =  2Kw/3Kp 
Kg  =  107.78  Eglg/l3 
Kc  =  85.33  EgIq/1 

o  o  o 


K. 


K. 


(K.  2) 


0.792 
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Replacing  in  K.  2 


W1 


1  kf  kf  kf2 

-  +  0.612  _  +  0.507  _  +  0.228  _ 

i  =  1.263 
c  K^c 


K  K  2 

+  1  +  2.525  ~  +  0.797  tA, 

Kg  Kg2 


L  . 


im 


o 


_JL 

3c 


1 H — " 
C 


L. 


im 


If  a  shearing  solution  is  sought 
FA  =  0.25  W1 


FB  = 


0.25  W1 


Thus 


V 


K. 


3c 


A 
W1 


6K  c 
bl 


v4 


K. 


K, 


+  1  kf2 


8  K„  2 


_1 

12  K„ 


I  +  ft 

c 


2K 

35  +  1  +  iC 


„  Kp2 
F  +  F 


2Kq  2 
1 


K  =  10.67 

O 


Gs*s 


K. 


GsAs 


5.33 


1 
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Thus 


K. 


Wl 


3c 


3KSC 


2K, 

3K, 


+ 


V 

2Ks2 


K, 


1 . 5cK, 


+  1  + 


4Kp  2Kf2 

V  +  Ks2 


L  . 


im 


->•  o 


Lim  Wl 
KS 


0.25 


-*■  00 


APPENDIX  L 


FIVE  BAYS  SYMMETRICAL  BUILDING 

Making  use  of  the  same  procedure  As  for  four  bays,  solution 
for  the  force  taken  by  the  frame  is  obtained  (see  Figure  L-l). 
Equilibrium  yields 


2V  +  2VA  +  2V_  =  2F  +  2FA  +  2F„ 
w  A  B  w  A  B 

VA  +  VD  +  V  =  F.  +  F_  +  F 
A  B  w  A  B  w 

i 

(Lei) 

Diaphragm  distortion 


Diaphragm  Shear  First  Span 


=  V  -  F 
D^  w  w 


V  =  F  +  FA  +  FB  -  V.  -  VD 
w  w  A  A  B 


D. 


A  A  ~  A 
A  w 


V 


D. 


K, 


VA  V 
A  _  __w 

K0  K 
F  w 


Moments  with  respect  to  wall 


VA  k  +  VB  !  L  +  VB  1  L  +  VA  I  L  +  VWL 


W1 

VA  +  VB  +  Vw  =  “2 
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Diaphragm  shear  second  span 


+  V, 


Diaphragm  distortion 


D, 


=  A 


B 


A 


V 


D. 


K. 


V 

K, 


B 


V. 

K, 


K, 


(F 


B 


V 


From  Lei 


K, 


V 

\ 

T 


w 


A 


B 


V, 


w 


K, 


K, 


K, 


V 


B 


K, 


The  set  of  simultaneous  equations  is  obtained 


+ 


]  +  V, 


B  K. 


w 


VA  +  VB  + 


V 

w 


W1 

2 


vb  [ri  +  ir-1 


vw 

vA 

VB 

- - -  1 

m< 

CD 

h 

Figure  L-l  Five  Bays  Symmetrical  Building 
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V 

+  FAo  + 

bl 


1 


K. 


B  K. 


+  F 


V 

B  KQ  2 
bl 


K. 


B  K, 


kf2  +  1  +  kf  1 
K„  2  c  2c 


_1 

2c 


v  o  K 

,  F  W1  Wl  _JF 

B  K  2  4c  4c  K 

bl  1 


1  + 


3K1 


kf2 

K„ 


-  +  I  +  - 

2  c  K 


K 


F  1 
2c 


APPENDIX  M 


NONSYMMETRICAL  THREE  BAY  BUILDING 


Referring  to  Figure  4*1  the  third  span  deflections  are 


A,  = 


4.35W14  ^  Wl2 
+  — — - — 


Wl  ,  Wl 

i  — - — —  “T 


384E  I  -  3nK  6K 

S  S  9GsAg  w  w 


V  1 2  VI3  VI  VI 

R  p  P  p 

8  1  _2  _6  2 

2  -  486  EsIg  486  Eglg  27  27  - 


^  i  *  |  \ 

+  - T - — - - — “ — B  i  — 


nK 


w 


K 


w 


= 


4, 35W1 

384  ESIs 


+ 


Wl 


9GsAs 


+ 


Wl 


6nK 


w 


Wl 

3K 


w 


A,  = 


VR  1 

7  1 

488  V7 


VR  1 

3  R1 


8 


27 


GSAS 


486 


VR  1 
R2 

EsIs 


VR  1 
6  R2 


27 


GsAs 


V 


R, 


2V, 


2  V 

i  ri  +  -9 


nK 


+ 


1  |  4  Vp 

9  9  2 


K 


w 


w 


KB  =  88 . 28  Eg I g/12  3 

Ksh  =  9  SV1 
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If  the  displacements  are  equalized  a  set  of  two  simultaneous 
equations  is  obtained,  furthermore  if 


x 


K 


B 


K, 


K 


sh 


Kw  (n+1) 
2K^. 


The  force  taken  by  the  first  frame  is 


^  -  0.181  lx2  +  1 . 183xy+y 2  +  <>•  0«039x(n  +  l )  +  0.1612  x  (n+1) 


2cn 


0 . 38  39y ( n+1 )  0.2222y(n+l)  0.167(n+l) 

X  Y  2cn  2c  2c 


0 . 0555 ( n  +  1 )  ^  0 . 333 ( n  +  1 ) 

+  — — - — - — -  *r  — 


n4c' 


2cn 


0 . 49  5x2  +  3 . 26  8xy  +  3y2  +  2.907x  +  4y  +  242x  <  n  +  1 ) 


2cn 


0 . 242x ( n+1 )  0.666y (n+1)  0,666y(n+l) 

2c  2cn  2c 


,  0 . 1111 ( n+1 ) 2  ,  0 . 5555( n+1 )  ,  0.5555(n+l)  ,  ,  _  A 

+  - 71 -  +  - 2, SH  +  25  +  1  ‘  B 

n4c 
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When  y  =  o,  pure  bending  solution 


181x2 


0 .04039x(n  +  l)  ^  0.1612x(n+l)  0.167(n+l) 

2c  2cn  2c 


0 . 0555 ( n  +  1 ) 

.  2 
n4c 


*  0 . 333 ( n  +  1 ) 

+  — —  ~  .  —  +  x 

2cn 


V 


R. 


V, 


0.495x2  +  2 . 9  07x  +  0  -,2,^2xlH±II  + 

2cn  2c 


+ 


0. llll(n+l) 2 

n4c  2 


+ 


0. 5555(n  +  l) 
2cn 


0 . 5555 ( n  +  1 ) 
+  - 2c - 


+  1 


V 


R. 


im  V, 


=  0 . 5 


c  -►  o 


V 


R. 


L  . 


0 . 1 8 1  x  +  x 


x 


im  Vr 


0.495x  +  2.907X  +  1 


2 „  7  25x  +  1 


C  4-  oo 


L 


im 


0.367 


X  4-  co 
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When  x  =  o,  pure  shear  solution 


„2  ,  „  0 . 3889y  (n  +  1  )  0.2222y(n  +  l) 

Y  Y  2cn  2c 


0 . 0555 ( n+1 ) 2 
+  -  - 

n4c 


0 . 333(n  +  l) 
2cn 


V 


3y 


2 


+ 


4  f  0  •  666y ( n  +  1 )  +  0 . 666y ( n  +  1 ) 

Y  2cn  2c 


0.5555 (n+1 )  ^  0.5555(n+l)  ,  . 

i  - - - - — - - - - —  T1  - - - - - - -  +  X 


2cn 


2c 


VT 


L. 


lm  V„ 


=  0.5 


c  >  o 


V 


R. 


L. 


y  +1 


lm  Vr 


c  -►  °° 


3y  +  4y  +  1 


3y  +  1 


V 


R. 


1m  Vr 


=  0.333 


->  00 


0 . 167(n  +  l) 
+  2c 


0 . 1111 (n+1) 2 
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When 


The  force  taken  by  the  second  frame  is: 


0 . 18 12  +  l.isixy  +  y2  +  P  •  1605x(  n  +  1).  +  0 . 3885y ( n+l> 


0 . 2222y ( n+1 )  ,  0.0409x(n+l) 

+  — — - — — - +  — — — — >  +  x  +  y 

2cn  2cn  * 


0. 05555( n+1) 2  0. 16666(n+l)  0.333(n+l) 

A  2  2nc  2c 

4nc 


0 . 493x2  +  3 . 2  6  8  xy  +  3y2  +  2G905x  +  4y 

0 . 2419x( n+1 )  ,  0.2419x(n+l)  ,  0.6666  ,  lix 

+  - ~2E~ - -  +  2cn - +  *  (n  +  1) 

0 . 1 1 1 1 ( n+1 ) 2  0 . 6666y ( n+1 )  0.5555(n+l)  0.5555(n+l) 

„  2  2c  2cn  2c 

4nc 

y  =  o,  pure  bending  solution 


+  1 
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0 . 181X2  +  °/1^Q5x(n  +  1) 

2c 


0 . 0409x( n  +  1 )  , 

- 23H -  + 


X 


+  0. 0555( n  +  1)  2  0 . 1666 ( n  +  1 )  0.333(n  +  l) 

.  2  2nc  2c 

4nc 


0.495x2  +  2 . 9  07x  +  -^---l3-^x^n  +  1I  +  0.2419x(n  +  l) 


2c 


2cn 


0 . 1 1 1 1 ( n+1 )  ,  0 . 5555 ( n+1 )  ,  0.5555(n+l)  ,  , 

“r  - - - - — — - - — - - — —  T  — " — ' — — — - - - —  i  - - — — - - - — —  T  X 


4nc 


2nc 


2c 


L. 


im 


c  -►  o 


0.5 


_^2  =  0.181x2  +  x 

lm  VT  0 . 495x2  +  2 . 907x  +  1 
c  ->•  00 


x 

2.725  +  1 


L. 


im 


0.367 


x:  +  00 
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when  x  =  o,  pure  shear  solution 


,2  „  0 . 3885y ( n+1 )  0.2222y(n+l) 

Y  Y  2c  2cn 


V 


R, 


+ 


0 „ 05555 ( n  +  1 ) 
4  nc  2 


0 . 1666 ( n+1 )  ^  0 . 333 ( n+1 ) 

T  — - — — ■ — — - —  +  — — 


2nc 


2c 


Vr 


3  2  +  4  +  0 . 6666y ( n+1 )  +  0.6666y(n+l) 

Y  -*■  9  r*  n  O  r* 


0e 1111( n+1) 

4  nc  2 


0 . 5555( n+1 )  ,  0.5555(n+l)  ,  , 

+  — — —  -K  —  ■■  — — — —  1 


2cn 


2c 


V 


L  . 
im 


r, 

"v! 


=  0.5 


c  -►  o 


V 


R, 


Lim 


c  00 


V, 


+  y 


3y  +  4y  +  1 


3y  +  1 


L  . 


im 


0.333 


y  + 


00 


APPENDIX  N 


BEAM  ON  ELASTIC  FOUNDATION  MODEL 
The  rigidity  of  the  frames  is  converted  into  resisting 
elements  distributed  in  the  spacing,  a,  therefore  a  rigidity  g^ 
per  unit  length  is  used  (See  Figure  3.13). 

The  uniformly  distributed  load,  will  be  shared  by  the  slab 
and  the  frames. 

Thus 

W  =  W  +  W  (N.l) 

F  S 

in  which 


Wp  =  load 

taken 

by 

the 

frames 

Wg  =  load 

taken 

by 

the 

slab 

For  the  slab 


dZ]  =  fdXl  +  I^L) 
dx  '  ^dx  v  ^-dx '  B 


where 


-r^-)  =  shear  angle 

dx  Jv 


=  bending  angle 

ldx  JB 


Differentiating  once 


,2 

dx  v 


) 

B 
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H)  ■ 


'd 
dx"  B 


M 


ESIS 


and 


,2 

(Hi 

dxz  v 


GS  As 


w 


Substituting 


Bolc 


,2 

d  y 

dx2  G^A 


W 


M 


s 


S  S 


EsIs 


or 


+  M 

,2  El 
dx  S  S 


Bo* 


GSAS 


w 


Differentiating  twice  we  get  consecutively 


dM 


^Vs 


Bo* 


GsAs 


dV 


dxEszs 


®oK  d 


GsAs 


dV 

dx 


d4y 

W  =  E  l  £  +  _ 

S  b  dx  G_A 


3  k  d  W 
o  s  p  T 

2  S  S 


S  S 


dx‘ 


(Nc  2) 


From  N.l 


W  =  W  + 

s  F 


constant 


dW 

dx 


dW 


0 


dx 


+ 


awp 

dx 
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dw  dVJ 

_ _ s  _  _ _ F 

dx  dx 


Therefore 

2  2 
d  W  d  w 

s  F 

2  •  2 

dx  dxz 


(N.3) 


The  force  taken  by  the  unitary  spring  is 

Wp  =  gp-y 


Thus 


dWp 

dx  ^F 


d£ 

dx 


Differentiating  again 


2 

d  WT 


dx‘ 


'F 


,  2 

,  2 
dx 


d2W 


Making  use  of  N.3 


dx 


.2 

d  y 

2  =  -g  F  2 

^  dxz 


Substituting  in  (N.2) 


,  4  3  <EI  ,2 

w  -  F  t  _  o  S  d_jy 

Ws  ■  ESZS  .  4  „  T  gF  .  2 

dx  GsAs  dx 


and  finally: 


d^y  . 

dx2  G„A„  9F  dx2  ESTS 


.4  3  k 

d  y  _  o 


W 


S  S 


(  N.4) 
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where 


$Q  =  coefficient  of  inelastic  shear  deformation 
k  =  shear  deformation  coefficient 
g_  =  unitary  rigidity 
G  =  shearing  Modulus  of  slab 

O 

E  =  Young's  Modulus  of  slab 

O 

I  =  Moment  of  Inertia  of  slab 
A  =  shear  area  of  slab 

o 

W  =  Uniformly  distributed  horizontal  load 


If 


26 


2 


Bo< 


GsAs 


and 


'4 

a 


g, 


V 


s 


Replacing  in  N.4 


dx4 


=  26 


2  d2y 
dx2 


+  a  y  = 


W 


E„  I 


(N.5) 


S  S 


Solutions  for  this  equation  are  dependent  upon  the  value  of  a 


relative  to  6 
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1.  When  a  >  8 


y 


dx 

e 


( Acosfx 


+ 


Bsinf x) 


+ 


-dx 

e 


(Ccosfx  +  Dsinfx) 


where 


d  = 


and 


2  e  when  8  >  a 


y  =  Ae(d+f>X  +  Be-(d+f)x  +  ce(d-f)x  +  De-<d-f,X  + 

gF 


where 


3.  when  8  =  a 


y  =  Aedx  +  Bxedx  +  Ce-dx  +  Dxe-dx  + 


APPENDIX  P 


ONE  STOREY  BUILDING  WITH  END  SHEAR  WALLS 


1.  a  >  3 


The  Boundary  conditions  are 
Q_ 


x  =  1/2 
0 

y  = 


w 


K 


w 


K 


w 


Mg  =  0  (Neglecting  torsional  resistance) 

Symmetrical  conditions 
y  (x)  =  y  (-x) 

y  =  2A  cosfx  coshdx  +  2B  sinfx  sinhdx  +  W/g^, 

=  2  sinfx  coshdx  (Bd-Af)  +  2  cosfx  sinhdx  (Ad+Bf) 
dx 


~  =  2  cosfx  coshdx  ( Ad 2+2Bdf -Af 2 )  +  2  sinfx  sinhdx 
dx2  9  9 

( Bd  ~2Afd”Bf  ^ ) 


, 3  ? 

— =  2  sinfx  coshdx  ( Bd  -3 Af d 

dx  9  99^ 

( Ad  +3Bfd  -3Af  d-Bf  ) 


3Bf  2+Af  ^ )  +  2  cosfx  sinhdx 


Bending  moment  of  slab 

2 


M  (x)  =  -EqIq  "232y  +  23 

s  s  s  dx2  9 


2  W 

F 


X  =  1/2  0s(l/2)  =  -  ysw  =  -0W 
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Shear  of  slab 

°s  ■  -Vs  (“3  -  2s2 


El  ,3 

v  =  _J>  rd  y 

Y(1/2)  Sw  ldx3 


232  &) 


2Acos^4cos-^4  +  2B  sin^-ir  sinh^-i  +  — 
2  2  2  2  g 

r 


Es1s  ro  .  fl  ,  dl  ,  , 3  -  A -  ,2  onr2j,.r3. 

— g—  [2 si n — ^  cosh—  (  Bd  ~3 Afd  -3Bf  d+Af  ) 
w 


+  2  cos— ^  sinh^—  (  Ad 2+3Bfd 2-3Af  2d-Bf  2  ) 


~„2  ,  .  fl  ,  dl  A,,  ^  0  fl  ,  dl 

20  (sm—  cosh— ^  (Bd-Af)  +  2  cos— j  smh— 


Substituting 


P 

0 

M 

N 

R] 

S 

V 

T 


.  fl  ,  dl 

sin—  cosh— 2 

d3  -  3f 2d  -  2g2d 

fl  .  .  dl 
cos—  smh— 

f3  -  3  f d2  +  2  02  f 

d2  -  f2  -  2  $2 

tg-^|  tgh^| 

fl  .  , dl 
sm—  smh— 

fl  udl 

cos—  cosh— 


(Ad+Bf ) ) ) 
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Thus 


K 


w 


r  ~ ,  fl  ,dl  _  „  r±  ,  ai  w  >1 
[2Acos—  cosh— rr  +  2B  sin—  smh— -  +  — J 

h  l  g  ^  2  2  2  2  g 


fl  .  , dl 


W 

[F 


2A  (NP  +  MQ)  +  2B  ( PQ  -  MN ) 


K  fl  HI  K 

A  (  NP  +  MQ  -  cos— r  cosh—)  +  B  ( PQ  -  MN  - 

X  ^  Z  Z  X 


s  s 


s  s 


■p  1  HI  K  ©  W 

sin-=-=-  sinh^i)  =  ,  w 

2  2  29FESIS 


Applying  M  (1/2)  =  0 


_  _  fl  ,dl  ,.,2  2.  ,  _  .  fl  .  ,dl 

0  =  2  cos—  cosh—  (Ad  +2Bfd-Af  )  +  2  sin— ^  smh— 


(  Bd2-2Afd-Bf 2  )  -  2  02  (  2A  cos^~  cosh-XX  +  2B  sin^X 


•  u d  1  ^  W,  ,  a2  W 
sinh-r  +  — )  +  2$  — 

2  gF  g  f 


222  fl  dl  fl 

A(d  -f  -20  )  +  2Bfd  -  2Afd  tg-^  tgh-|  +  B  tg— 


.  ,  dl  .2  ^ 2  0  2 )  _  n 
tgh—  (d  -f  -2  0  =  0 


A  =  B  ( 2 f d  +  SR1)/(2fds  -  Rx) 


B  = 


KT7.W  (  2fds  -  R ) 
W  _  1 

2gpEIs  Hx 
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K  .  W  ,  -  _  , 
w  ( 2  fd+  SR 1 ) 

2gFESISHl 


where 

K  K 

H,  =  (2fd(  PN+MQ4-SPQ-SMN~~i  -  ) 


+  R1(SPN+  SMQ-PQ+MN)) 


The  solution  is 


y 


K 


w 


('1+HiEsI 


( 2  f d+SR,  )  cost x 
S  1 


coshdx 


K 


+ 


w 


uiEszi 


sinfx  sinhdx)  (P.l) 


2. When  3  >  a 


d  = 


2  2 
+  3 


y  =  Ae<d+f)x+Be-(d+f)x+Ce(d-f)x+De-(d-f)X+  ^ 


( 2fdS-R1 ) 


^  =  (d+f)Ae(d+f)x-(d+f)Be-(d+f)x  +(d-f)Ce(d-f)x 


( d-f ) De 


( d-f ) x 


^  =  (d+f)2Ae(d+f)x+(d+f)2Be-(d+f)x+(d-f)2Ce(d"f)x 
dx 

+  (d-f ) 2De~(d+f )X 


=  (d+f ) 3Ae(d+f )x-(d+f ) 3Be  ( d+f } x+ ( d-f ) 3 

dx3 


Ce(d-f)x-(d-f)3De'(d-f)x 


Symmetrical  conditions 


y  ( x )  =  y  ( -  x ) 


Thus 
y  (  x )  = 


2A  cosh ( d+f ) x  +  2B  cosh(d-f)x  +  W/gF 
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Boundary  Conditions 


1.  x  =  1/2 
Ms(l/2)  =  0 

2 

Ms  =  ”ESIS  {^~2  “  +  2^2w/gF) 

dx 

A(d-f)2  cosh  (d+f)  1/2  +  B(d+f)2  cosh(d-f)l/2  =  0 

(P.2) 

2.  x  =  1/2 

y ( 1/2) 


Vs 

Kw 


■diy 

dx' 


2B 


2  dy 
dx 


a( (d+f ) 2sinh(d+f  )  1/2  -  202(d+f)  *  sinh(d  +  f)l/2 


-  cosh(  d  +  f  )  1/2  •  +  B  ( ( d-f  )  2  sinh(d-f)l/2 

ESIS 

?  l  K 

-  20^(d-f)  sinh  (d-f)l/2  -  cosh(d-f)  p— y— ) 

Z  t-glg 
WK 

_  _ w 

2gF  ESIS 

Replacing 


(P.3) 


sinh  (d+f)  1/2  =  P2 

cosh  (d+f)  1/2  =  RM2 

sinh  (d-f)  1/2  =  VI 

cosh  (d-f)  1/2  =  T2 
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d  +  f  =  R2 
d  -  f  =  RN1 

W  K 

A  =  —  .  w  _ 
2gFESIS 


R2T_  (  R3 P  -2  e2R0  P0  -RM  S  /E  l  ) -RN2  RM„  (  RN3V,  - 
Z  Z  2  2  Z  Z  Z  W  b  b  1  Z  11 


2TRN1V1 


T0  K  /E  I 

2  w  S  S 


B  = 


3.  When  3  =  a 

y  =  Aedx  +  Bxedx  +  Ce”dx  +  Dxe-dx 


Symmetrical  conditions 
y ( x)  =  y  ( -x) 


Therefore 


y  =  2A  coshdx  +  2B  sinhdx*x  +  W/gp 


Boundary  Conditions 


1.  x  =  1/2 

Ms ( 1/2 )  =  0 

A(d2cosh~  -  2b2cosh~)  +  B(d2^  sinh™  + 


2d  cosh—  -  2b2  sinh^-)  =  0 


2.  x  —  1/2 

y  (1/2,  -  Wf-  2B2^) 

w  dx 

A(d2sinh-^  -  2  32dsinh-—  -  cosh-^-  ^ - ) 


nf J  1  , dl  t  0,2  .  ,  dl  0„2dl  Kdl 

+  B(d  ~2  cosh— ^  +  3d  sinh— ^  -  23  —  cosh— ^ 


2b2  sinh-^  -  sinh-^-i  — 


K 


w 


sIs 


)  » 


WK 


w 


2<3pBs^S 
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Replacing 


cosh—  =  rm3 


•  i  ci  1 

Slnh—  =  P3 


2 

V  =  2dRM  -  -  p 

2  3  2  3 


T3  =  d  RM 3 


R3  =  P3d3  +  ^3  eTT 


s  s 


d2l  2  K 

RN0  =  -  — ~  RM0  +  d  P0  -  P-3  -=r^r~ 

Z  Z  i  6  3  Ig 


Therefore 


wK  V0 
w  2 


A  = 


2gp  Eglg  (T3RN2  R3V2^ 


B  = 


AT. 

~V, 


APPENDIX  Q 


ECCENTRICALLY  ARRANGED  SHEAR  WALLS 
The  solution  is  divided  in  two  branches  as  it  is  seen 
Figure  4.4C 


1 . 


a  >  3 

dx^  -dx^ 

y,  =  e  1 ( A. cost x. +B, sinf x. )  +  e  ±(C,cosfx» 
1  1111  11 


+  D, sinfx . )  +  — 

1  i  gF 


y 


2 


dx2  ”^x2 

e  1 ( A^cosf x „+B„s infx _ )  +  e  (C0cosfx 
2  2  2  2  2 


+  D^sinfx^) 


The  Boundary  conditions  are; 


in 


187 
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3. 


K 


w 


2 


4  . 


( 


GsAs 


GsAs 


J 


that  is  the  rotation  angle  is  the  same 


6 «  X1  =  x2  =  0 


M 


2 


that  is  the  Moment  is  continuous 


7. 


2 


8 


=  0 
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Thus,  the  constants  to  be  determined  are: 


The  bending  moment  equation  is 


2B2y  + 


The  shear  equation  is 


dl 


A1[cosfl1e  1(d2-f2-2$2) 


dl 


-  2  fde  ^s inf 1 ,  ] 


dl,  oo  9  dl, 

+  B1  [  sinf 1^  e  i(d  -f  -23  )  +  2fde  cosfl^ 


+  [cosfl^  e 


-dl 


12  2  2 
1  ( d  -f  -2  3  1 


-dl 


+  2fde 


s  inf  1 


-dl 


-dl 


+  D-, 


[sinfl1e  1(d2-f2-232  )  -  2fde  1cosfl1  ]  = 


0 
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Making  the  following  substitutions 


dl 


V 


sinf 1 


dl 


cosf  1  ^  e 


=  sinfl^e 


=  cosfl^e 


-dl 


-dl 


2  2  2 
R  =  d  -r  -2 3 


The  above  written  expression  becomes 


A1  [TiR-2fdSi  ]  +  B1  (siR+2fdT1 ]  +  C1 [V 


+  D1  [uiR-2fdV1  ]  =  0 


Another  set  of  substitutions  is 


TXR  -  2fdS][  =  P1 

S1 R  +  2fdT1  =  Q1 

V1R  +  2fdU1  =  R} 

R  -  2fdVx  =  Wx 


1R+2fdU1 ] 
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Thus 


A1P1  +  B1Q1  +  C1R1  +  D1W1  =  0  (2) 


Condition  8  states 


x0  =  M  =  0 

2  2  S2 


A, [cosfl,e  2(d2-f2-2B2 )-2fde  2sln£12] 


d  d  1. 

+  B9[sinfl  e  (d  -f  -2$  )  +  2fde  cosfl9] 


+  C2[cosfl2e  2(d2-f2-2g2) +2fde  2  s i n  f 1 2  ] 


-dl. 


d  1 9  9  9  9  d  1 9 

+  D9  [sinfl?e  id  -f'-2g  )~2fde  cosfl2]  =  0 


Making  substitutions 


■dl 


So  = 


sinf 1 2e 


dl 


To  = 


cost 1 2e 


-dl 


Uo  = 


sinf 1 2e 


•dl 


Vo  = 


cosf 1 2e 


Thus 


A2 [T2R"2fds2 ]+b2 [S2R+2fdT2 ]+C2 [V2R+2fdU 
+  D2 [u2R-2fdV2 ]  =  0 
P2  =  T2R  -  2fdS2 

02  =  S2R  +  2fdT2 

R2  =  V2R  +  2fdu2 

W2  =  U2R  -  2fdV2 

Thus 

A2P2  +  B2°2  +  C2R2  +  °2W2  =  0 

Condition  6  states 
X1  =  x  2~  0 


M 


M 


s 


1 


s 


2 


M 


2(J2y  +  2B2  -^] 

y  p 


F 


M 


s 
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Ax (d2-f 2-2B2)  +  B12fd  +  Cx (d2-f 2-232) 

+  D1(-2fd)  = 

A2(d2™f 2~202)  +  B22fd  +  C  2 ( d  2-f  2- 2  02  ) 

+  D2(-2fd) 

AXR  -  A2R  +  B12fd  -  B22fd  +  C1R  -  C2R  -  2DX fd 
+  2D2 f d  =  0 


AXR  -  A2R  +  B12fd  -  B22fd  +  CjR  -  C2R  -  D^fd 


+  D2  2fd  =  0 


(6) 


Condition  5 


X1  =  x2  =  0 


-  -4-  °s  )  ■ 

ax  (i  a 

S  1  1 


fdy  _ 

ldx 


-4-  Qs  ) 
GSAS  2  2 


dG  CAC  f G pAp  dGQA 

Mx„+  F  t  «— )  -  Bi(yc-  e  i  cBx+«  F  t  J 

S  S  o 


^  °  Vs6o 


e  e  i  ic 

o  S  S 
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fG«Ac  fGcAc  fG„A„ 

Dl<Y°E^W)  =  'A2(xo+EsisBo<)  +  B2(VesIs80^ 


dGcAc  fGeA0 

+  C  0  ( x  +--  1  +  D  f  y  — — — — — —  1 

2  OSoESISkJ  2lY°  WokJ 


In  which 


dGsAs 

X°+  ESISB  o  El 


3  2  2 

xQ=  d  -3df  -20  d 


fGsAs 

y°  VsV  1 


3  2  2 

y Q=  f  »3d  f  +  2Tf 


Thus 


A1E1+A2E1  B1F1  B2F1  C1E1  C2E1  D1  F1  °2  F1  ° 


(5) 


Condition  3 


+  0. 


=  0 


Yt  = 


K 


w 


-A, G , — A0X  -B , y  +B 0y  -C , H . +C0X  +D,y  +D0y 
11  2o  ho  2J  o  11  2o  lJo  2J( 


K  •  W 
w 


gfE  I 
s  s 


(3) 
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in  which 


Ktt 

G  =  x  +  ^ 

1  o  I 


S  S 


K 


H 


w 


3  EsIs  ° 


Condition  4 


x  =  0 


Q  +0 
S1  S2 


K 


w 


KwW 


.  x-A  G  +B  y+B  y+C  *-C2H1+D  y+Djy  =  E  t  <4> 


F  S  S 


in  which 


K 


w 


X  +  E  I 

1  is 


„ 

1  EsIs"  X 


Condition  1 


X1  =  31 


0, 


K 


w 


T3  =  d3T1  +  f3S1-  3d2fS1  -  3df2Tx 


S3  =  d3S1  -  f3T3  +  3d2fT1  -  3df2S3 


V3  =  d3VT  -  f3U1  +  3d2fUL  -  3df2V1 
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U3  = 


d3Ux  +  f3V1  -  f3Vx  -  3df2U  -  3d  2f  V. 


T 


4 


fS 


1 


s4  =  dsx  +  f T  x 


v4  =  dv2  +  fu4 


u4  =  fv4  -  dux 


Also 


?  Tl  K 

T5  *  T3  "  20  T4 


S5  =  S3  -  23  S4  - 


V5  =  V3  -  23  V4  + 


E_ 

T 

S 

s 

Si 

K 

1 

w 

E_ 

T 

S 

s 

V, 

K 

1 

w 

E_ 

T 

S 

S 

U, 

K 

1 

w 

U5  =  U3  +  20  °4  +  E.I 


s  s 


WK 


A,TC  +  B ,  S  c  -  C,Vc  -  D-,  U  (- 
15  15  15  15 


w 


gFES  TS 
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Condition  7 


x 


2 


K 

w 


K 

w 


dl 


dl 


A2cosf 12e 


‘+B2sinf  1 2e 


'+C  2cosf  1  2e 


+ 


-dl 


D2sinf 1 2e 


EsIs 

K 


[A2(d3T2+f 3S2~3d2fS2”3df 2T2) 


+  B2(d3S2-f 3T2+3d2fT2-3df 2S2 
-  C2(d3V2-f3U3+3d2fU2-=3df 2V2) 
-D2(d3U2+f 3V2-3df 2U2~3d2fV2) 
-232A2(dT2-fS2 )-262B2(dS2+fT2 ) 
+  2  $2c  2  ( dv  2+ f  u  2 )  -  2  e  2d  2  (  f  V  2»dU  2 ) 


-dl 


2 


In  which 

Tg  =  d3T2+f3S2-3d2fS2-3df 2T2 
S  =  d3S2-f3T2+3d2fT2-3df2S2 
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V6  =  d3V2-f 3U2+3d2fU2-3df 2V2 

U  =  d3U0+f3V.-3df 2U0-3d2fV^ 
o  2  2  2  2 

T8  =  dT2-fS2 

S8  =  dS2+fT2 

V8  "  dV2+fU2 

U8  -  fv2  -  du2 


Thus 


A2(Tg-2$  Tg-T2  E  t 


— >  +  B2(S6-2b2s8-S2  E~I 


s  s 


s  s 


-  C2(V6-282V3+  -  D2(U6+282V8+U2  -g^ 


s  s 


s  s 


W  K 
_ w 

9FES  IS 


in  which 


?  T?Kw 

T6  -  28  T8  -  =  T9 


S6  "  2b  S8  “  S 


Kw 

2  Vi =  89 
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v6  -  2B%  +  v2  o:  ■  V9 

D  O 


u6  +  2B\  +  U2  =  u9 


W  K 


w 


A2T9  +  B2S9  -  C2V9  -  °2U9  =  = 

rob 


( M.  7 ) 


Thus  a  system  of  8  simultaneous  equations  is  obtained 


W  K. 


w 


(1)  A,T,-  +  B,SC  -  C,Vc  -  D,VC  =  -  _  _ 

15  15  15  15  Eg  Ig 


(2)  Aipi  +  B1Q1  +  C1R1  +  D1W1  =  0 


W  K. 


( 3 )  -A, G , -A0x+B , y+B^y-C , H , +C  0x+D , y+D0y 


w 


11  21^  21  1  1  2  ly  21  o  E  I 

y  F  S  S 


W  K. 


w 


(4)  -A^-AjG^B^+C^-CjH^D^+Djy  =  — E  , 

r  b  b 


(5)  aie1+a2e1-b1f1-b2f1-c1e1-c2e1-d1f1-d2f1  =  0 


(6)  A1R1-A2R+B12fd-B22fd+C1R-C2R-D12fd+D22fd  =  0 


W  K 


w 


(7)  A2T9  +  B2S9  +  C2V9  D2U9  pi 

r  b 


(S)  a2p2  +  b2o2  +  c2r2  +  d2w2  =  0 


Using  any  program  that  solves  a  set  of  simultaneous 
equations ,  the  constants  A^  ,  B^  ,  ,  A2  ,  B2  ,  C2  ,  D2  would  be 

known  and  the  values  incorporated  in  the  solution  for  y-^  and  Y2* 
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2 .  3  >  a 


The 

solution  is 

*1 

=  Ale(d+f)x  +  Ble'(d  +  f)x  +  Cle(d-f)x  +  D^e-^  d~f ) x 

^2 

=  A2e(d+f)x  +  B2e-(d+f)x  +  C2e(d-f)x  +  D2e-(d-f)x 

The 

Boundary  conditions  are 

1.  Xx 

=  1,  y,  =  Q  /K  =  -  Q  /K 

1  1  w  w  sl  w 

rH 

X 

• 

CM 

■  h  %  ■  0 

3 . 

■  0  *i  ■  Sc  ■  +  °s  >/Kw 

w  12 

4.  x2 

-  0  ^2  ■  VKw  ■  <%  +  %)/Kw 

5  .  x  ^ 

=  x2  ■  0  %  =  Ms2 

6 .  x 

-f  0, ) 

bSAS  1  ^SAS  2 

°s 

7*  x2  ^2  ^2  °w/Kw  K 

w 


xs  =  d2 


8 


0 
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Bending  Moment 


Ms  =  -Vs  (^“2  '  2s2y  +  21)2  W/gp) 

dx 


Shear 


°s  “ 


EgT  A  -  282  ^ 
dx 


dx 


) 


X1  11 


M  =  0 
S1 

9  ( d  +  f ) 1 ,  9  -( d+f ) l  (d-f)l 

(d+f)zA1e  x  +  (d+f)zB1e  1+  (d-fPC^e 


- ( d-f ) 1 ,  (d+f ) 1  -( d+f ) 1 1 

+  (d-f)  D1e  -23  (A.^  +B] e  1  + 


(d-f )1 


(d-f )1 


Cf  e 


+  e 


1  +  W/g  ) +  2b2  W/ g  =  0 


(d+f)l1  9  9x 

Axe  1 ( (d+f)  -2$  '  + 


-(d+f)l,  9  9 

+  B,e  1  ( ( d  +  f )  -23  J  + 


(d-f )1 


+  Cj^e 


1 ((d-f)  2 - 2  6  2  ) 


-  ( d-f  )  1  9  9 

+  D,  e  ((d-f) -23  )  =  0 
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8  X2  =  X2  %  =  0 
A2e(d+f } 12 ( (d+f ) 2-  202)  + 

e“(d+f } £2((d+f ) 2  -  232)  = 

C2  e“(d~f  ) 1 2  ( (d-f ) 2~2g2 )  + 

°2  e_(d“f } l2{ (d-f ) 2  -  2 B2 )  =  0 

V 

5  X1  =  x2  =  °  Msx  =  Ms2 

A1  ((d  +  f )2  -  2e2  )  + 

B1((d+f)2  -  232)  + 

D1  ((d-f)2  -  232)  = 

=  A2  ((d  +  f)  2  -  2 3 2  )  + 

B2  ((d+f)2  -  232)  + 

C2  ((d-f)  2  -  2  3 2  )  + 

D2  ((d-f)2  -  232) 
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3  X1  =  0  *1  -  VKw  ■  (QSl+  °s2)/Kw 


-d  y  2  dy 


°s  ESJS  (  3  26  dx 

1  dx 


0=  =  -ESIC  (^4  -  2B2  |£) 


s0  S  S  ^  ,  3 

2  dx 


Qs  =  -EsIs((d+f ) 3A1  -  (d+f ) 3B1+(d-f ) 3C1-(d-f ) 3D1 


2  32  ((d  +  f  )  Ax-(d+f  )B1+(d-f  )C1-(d-f  )D1 ) 


Qs  =  -Eg Ig  ((d+f) 3A2-(d+f ) 3B2+(d“f ) 3C2“(d-f) JD; 
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£ 
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K 
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1  ESZS 


+  A2  [(d+f)3  -  2 B2 ( d  +  f ) ]  +  b2  [-(d+f) 3  +  2B2(d+f )] 
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K  w 
w 


gFESIS 
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R 
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K 
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s  s 
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K  W 
w 
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1  xi  =  1i 


Q, 

l  =  ~  F 


w 


( d  +  f ) 1  - ( d  +  f ) 1 ,  ( d-f ) 1  - ( d-f ) 1  w 
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11  1  i  gt 
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K 


w 
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K 
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+  D1  (- ( d-f ) 3e  (d  f)11 


K 


W  e"(d"f)11  +  23  2(d-f)  e  -(d"f)11) 


Vs 


W  K 

w 


gF  ESIS 


7 
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W  K 


w 


gF  EgIs 
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3<  2b  E  I 
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+  D1  (-(d-f) 


0  <  ,  202E  I 

E0I0  (d-f)J  +  — - eo<  (d-f)) 


GsAs 


s  s 


GsAs 
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GsAs 


GsAs 


3  =  a 


The  same  boundary  conditions  apply 

i 


M  =  0 
S1 


A1d2edll  +  B1d  edll(dl1-2)  +  C1d2e  dll  + 


+  D  d  e  dll  (dl1+2)  =  0 


0 


A2d2edl2  +  B2d  edl2  (dl2~2)  +  C2d2e  dl2 


+  D2d  e  d 1 2  (dl2+2)  =  0 


0 
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Axd2  +  B1  d  (dl  1~2)  +  C  ^d  2  +  D^fdl  +2) 
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3  <  9  e  k 
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APPENDIX  R 


1.1  DESIGN  OF  A  THREE  BAY  ONE  STOREY  BUILDING  WITH  END  SHEAR 
WALLS  AND  CONCRETE  HORIZONTAL  DIAPHRAGM. 

The  building  is  shown  in  Figure  5.1  and  is  comprised  of 
three  bays,  two  end  shear  walls  and  two  interior  frames. 

The  type  of  deck  used  is  selected  from  Reference  33  which  is 
T-15  Hi  Bond  24"  coverage,  deck  thickness  0.030"  and  a 
concrete  cover  of  2.5"  (see  Figure  5.2). 

The  design  loads  are  shown  below: 

Vertical  loads: 

Dead  loads: 

Deck:  35  psf 

Beams:  15  psf 

Live  loads: 

Snow:  50  psf 
Horizontal  Loads 

The  procedure  recommended  by  the  National  Building  Code  of 
Canada  (25)  states: 

V  =  A  S  K2  I  F  W  (  R.  1  ) 

where : 

V  =  Minimum  lateral  seismic  force 
S  =  Seismic  response  factor 

A  =  Acceleration  ratio 

K2  =  Constant  that  represents  damping  and  ductility  of  the 
building 

F  =  Foundation  factor 

W  =  Dead  load  including  25%  of  the  design  snow  load 
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2  14 


W  =  (50  +  0.25  x  50)  x  90  x  15  =  84375  lb. 
Therefore,  replacing  in  (R.l) 

V  =  0.08  x  1.0  x  0.8  x  1.0  x  1.0  x  84375  =  5400  lb. 
Thus,  the  uniformly  distributed  horizontal  load  is 
V/90  =  60  lb.l.f. 


Stiffness  of  the  Horizontal  Diaphragm 

To  calculate  the  stiffness  of  the  horizontal  diaphragm  three 
constants  must  be  assessed. 

1.  Moment  of  Inertia  of  the  perimeter  beams  with  respect  to 
the  centerline  of  the  building. 

2.  Moment  of  Inertia  of  the  transformed  area  of  steel  deck 
with  concrete  fill. 

3.  Shear  area  of  the  horizontal  deck  with  concrete  fill. 


1.  Moment  of  Inertia  of  the  perimeter  beams 
Beams:  W8  x  28  Area  =  8.23  in2 

V  2A 

where 

Is  =  Moment  of  Inertia 

A  =  Cross  sectional  area  of  the  perimeter  beams 
b  =  width  of  the  building 


Thus 
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15  x  12 
2 


2 


Is  =  2  x  8.23 


133,326  in 


4 


Moment  of  Inertia  of  the  transformed  area  of  concrete 


1 . 5 


33  y 


/f  'c 


E 


c 


where 

Y  =  145  pcf 
f 1 c  =  3000  ps i 

Thus 

Ec  =  3,155,924 
and 

n  =  E/E  =  29,000,000/3,155,924  =  9.2 

O 

\ 

For  a  1  ft.  width  the  steel  area  is  =  0.499  in^  and  the 

.  9 

concrete  area  =  36.75  m. 

The  equivalent  steel  depth  of  the  slab  is 


t 


+  0.499  j/12  =  0.375" 


Therefore  the  moment  of  Inertia  of  the  horizontal  diaphragm 


3 


I s=  t  b3  =  —  0.375  x  1 8 0 3  =  182,735  in4 
Shear  area  of  the  steel  deck  with  concrete  fill 
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To  calculate  the  shear  deflection  of  the  deck  the 
flexibility  factor  is  used  (33)  such  that 


V  L  F 

S  2  x  106 
where 

V  =  Average  Shear  along  length  L/2 
L  =  span 

F  =  Flexibility  factor 


(  R.  2) 


but 


V 


wl 

4b 


(R.  3) 


Thus 


A  =  ^L_L _ £ _  (  R.  4  ) 

s  8b  x  106 

The  shear  deflection  for  an  equivalent  beam  is: 


A 

s 


wl ' 


8GsAs 


where 

G  =  shear  modulus  of  the  horizontal  diaphragm 
S 


(  R.  5) 
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shear  area  of  the  horizontal  diaphragm 

Web  area  of  plate  girder 

A  =  b.  t 
s 

where 

b  =  width  of  the  horizontal  diaphragm 
t  =  thickness  of  the  horizontal  diaphragm 

Comparing  R.4  and  R. 5  the  thickness  of  the  eguivalent  slab 
is  obtained  as 

t  =  io6/fgs 

From  Reference  33  the  flexibility  factor  F  is  0.90  for  the 
6'  spacing 


Thus 


1 , 000  ,  000 
11, 200,000  x  0.90 


0.099" 


The  shear  area  is 
A  =  0.099  x  15  x  12  =  17.9  in2 

O 

Frame  Design 

Designing  the  interior  frame  for  the  vertical  loads,  a 
minimum  size  for  beam  and  columns  from  strength  and  deflection 
point  of  view  is  obtained.  Thereafter  any  size  of  beam  and 
columns  may  be  used  with  the  results  shown  in  Table  R-l. 
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Table  R-l 


Stiffness  of  the  Frames 


Beam  IB  Column  Ic  (-^)B  (— K 

Kips/in 


W8  x  48 

184 

W6  x  20 

41.  5 

1 .02 

0.346 

3.574 

W10  x  49 

273 

W8  x  20 

69.4 

1 .52 

0 . 578 

5 .867 

W12  x  58 

476 

W10  x  29 

158 

2.64 

1.317 

12.740 
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Stiffness  of  the  End  Walls 

The  exterior  walls  may  be  made  of  any  material;  for  this 
particular  case  concrete  is  used. 

The  stiffness  of  the  exterior  wall  is  obtained  from  equation 

R-6 . 


K 


1 .  2 
G 


(  R-6) 


where 


e  =  8"  (thickness  of  the  wall) 
b  =  180"  (width  of  the  building) 

H  =  120"  (Height  of  the  building) 
E  =  3,156  k/in  (Young's  Modulus) 
Gg  =  1,372  k/in  (Shear  Modulus) 
v  =  0.15  (Poisson's  ratio) 


Solving  K  =  8346  k/in 


The  force  taken  by  an  interior  frame  is: 


K 


+ 


sh 


2.725 


K 


+ 


B 


Thus 


^R 

VT 


643  1661  28346 

2.725  ,  3,1 

643  17T6T  8T46- 


+ 


1 

3574 
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rr~  =  5.4  kips 
T 


VR  =  0.042  kips,  Vw  =  2.66  kips 


c  =  K  /K  =  8346/3574  =  2335 

w  r 


The  displacements  of  the  frame  and  wall  are  respectively 


A  _  0.040 
1  3,574 


0.0112" 


A  0  =  =  0.0003 

2  8346 


If  the  horizontal  diaphragm  is  considered  infinitely  stiff 


F  =  k  A 


where 


F  =  horizontal  force  applied  to  the  building 
A  =  Uniform  displacement  of  the  building 
A  =  5.4/16699  =  0.000323" 
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The  force  taken  by  a  frame  and  a  wall  are  respectively 

VR  =  3.574  x  0.00032  =  0.00116  kps 
Vw  =  8346  x  0.00032  =  2.699  kps 


and  c  i  =  00 


1.2  THREE  BAY  BUILDING  WITH  END  BRACED  WALLS  AND  CONCRETE 
DIAPHRAGM. 


Stiffness  of  the  End  Walls 

.2 

Kx  =  AE  -2- 
hi 

where 

A  =  Area  of  the  diagonal  element 
b  =  width  of  the  building 
E  =  Modulus  of  Elasticity 
h^  =  Length  of  Diagonal 

2 

K -i  =  2.70  x  29000  x  (15- 250  k/in 

( 18  x  12) 

The  coefficient  C  is 


x  =  250/3.574  =  69.9 


The  force  taken  by  an  interior  frame  is 


V 


1  +  1 


R  643  1999  2  •  250 


VT  2.725  +  _ 3 


1  +  —L 


643 


1999 


+ 


250 


3 . 574 
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Thus 


V 


R 


V. 


=  0.014 


VR  -  0.0757  Kips  =  2.624  kips 


The  displacements  of  the  frame  and  wall  are  respectively 


A,  = 


0.0757 


1  3.574 


=  0.0212" 


= 


2.623 


2  250 


=  0.0105" 


If  the  concrete  diaphragm  is  considered  infinitely  stiff 
displacement  of  the  building  is  uniform: 

4  =  IF1  ■  =  °-01065" 

and  the  force  taken  by  an  interior  fram  and  the  end  wall 
respective ly 

V  =  0.0381  kips 
R 

V  = 
w 


the 


are 


2.662  kips 
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Reducing  the  stiffness  of  the  braced  walls.  Let  use  a  frame 
having  a  stiffness  Kp  =  12.75  K/in  from  Table  R-l.  Also 
increasing  the  stiffness  of  the  interior  frames. 

Kp  =  12.75  K/in 

Thus,  the  force  taken  by  an  interior  frame  is: 


643  1999  2x12.74 

2.725  311 

643  1999  12.74  12.74 


0.254 


Thus 


VR  =  0.254  x  5.4  =  1.37  kips 
Vw  =  1.33  kips 
c  =  1 


The  displacement  of  the  frame  and  wall  are  respectively 


1 . 37 
1  12.74 


0.108" 


h2  =  0.104 


For  a  infinitely  stiff  diaphragm,  the  displacement  of  the 
building  is  uniform: 

5.4 
50.96 


A 


0.106" 
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The  force  taken  by  an  interior  frame  and  an  end  wall  are: 

VR  =  12.74  x  0.106  =  1.35  kips 

V  =  1.35  kips 
w  ^ 

2.  DESIGN  OF  A  THREE  DAY  ONE  STOREY  BUILDING  WITH  END  SHEAR 
WALLS  AND  STEEL  DECK  DIAPHRAGM 

The  same  building  is  designed  as  in  the  first  example. 

The  type  of  deck  used  is  selected  from  Reference  33  which  is 
T-15  24"  coverage,  3  transverse  welds,  span  between  purlins  is  6' 
and  the  deck  thickness  is  0.030". 

Horizontal  Loads 

The  procedure  recommended  by  the  National  Building  Code  of 
Canada  ( 25)  for  wind  states 


P  3  Ce  Cg  Cp 
q  =  10.5  lbsf 

ce  =  1.0 


Thus 

p  =  10.5  x  1.0  x  2.0  x  1.2  =  25.2  lbsf 
V  =  25 . 2  x  5  =  126  lb.  1 .  f . 


126  x  90 


11340  lb 
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Stiffness  of  the  Horizontal  Diaphragm 

Two  constants  must  be  assessed: 

1.  Moment  of  Inertia  of  the  perimeter  beams  with  respect 
to  the  centerline  of  the  building. 

2.  Shear  area  of  the  horizontal  disphragm 


1.  Moment  of  Inertia  of  the  Inertia  of  the  perimeter  beams  is 
calculated  below 

Using  W8  x  28  beams  from  Reference  8. 

The  area  is  8.23  in^  and  the  moment  of  inertia  is: 

2 

I  =  2  x  8.23  (1-|— )  =  133,326  in4 

2.  Shear  area  of  the  horizontal  disphragm. 

The  flexibility  factor  is  used  (33)  such  that 
F  =  23  +  185R. 

where 

R  =  Ration  between  purlin  spacing  and  length  of  deck 
supplied 

In  this  case 

F  =  23  +  185  x  6/18  =  84.7 

and  the  thickness  of  the  eguivalent  slab  is 
t  =  106/FG  =  106/84 . 7  x  11,200,000  = 


0.0013" 
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The  shear  area  is 
Ag=  0.0013  x  15  x  12  =  0.195  in2 


Stiffness  of  an  Interior  Frame 

Table  R-l  provides  values  for  the  stiffness  of  the  frame. 
The  value  selected  from  this  example  is  12.74  k/in. 


Stiffness  of  the  End  Walls 

The  exterior  walls  in  this  case  are  braced  walls. 
The  stiffness  is 


where 


A  = 


b 

E 

h 

A 

K 


1 


1 


Area  of  the  diagonal  element 
width  of  the  building 
Modulus  of  Elasticity 
Length  of  Diagonal 
2.70  in2 


2.70  x 


29000 


( 15x12) 2 


( 18x12) 


3 


250  K/in 
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The  force  taken  by  an  interior  frame  is 


V 

V. 


R 


K  +  %  +  2lT 


KB  = 


2 . 725  3  1  1 

Kr  Kch  K  Kp 
B  sh  w  F 

88.36  ~ 


Iq  =  133,326  in' 


Kb  =  88.36  x  29,000  x  13 3 , 3 26/ ( 9 0x12 ) 3  =  271  K/i 


n 


K  =9 
Ksh  y  l 


Ag  =  0.195  in' 


Ksh  =  9  x  11,200  x  0.165/90x12  =  18.2  K/in 


c  =  KW/KF  =  250/12.74  =  19.6 


271  18.2  2x250 

2.725  3  1 

271  18.2  250 


+ 


1 

12.74 


0 . 2356 


.  *  .  V  =  11340  x  0.2356  =  2672  lb 
R 


2.998  Kips 
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The  displacement  of  frames  and  end  braced  walls  are  respectively 


A 


1 


2.672 
12 .74 


0.210" 


A 


2 


2.998 

250 


0.012" 


If  the  horizontal  diaphragm  is  considered  infinitely  stiff 


F 

A 


KA 

11.34 
5  25.5 


0.022" 


The  forces  taken  by  frames  and  end  braced  walls  are  respectively: 
VR  =  12.74  x  0.022  =  0.275  kips 
Vw  =  5.500  kips 


Concrete  end  walls 


K, .  =  8  346 
w 

Kp  =  12.74 


27T  +  18.2  2x8346 

2.725  ^  3  J  1  ^  1 

271  18.2  8346  12.74 


c  =  655 


0.231 
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Frame  displacement 


A  _  0.231  x  1134 
1  12.74 


0.206" 


Considering  the  horizontal  diaphragm  infinitely  stiff 
F  =  KA 
K  =  16717 

A  =  =  0.0007 


Force  taken  by  the  frame 

F  =  0.007  x  12.74  =  0.009  K 


tt~  =  0.0008 

VT 


Frame  End  Wall 


Kw  =  12.74  K/in 
Kp  =  12.74  K/in 


271  +  18.2  +  2x12.74 

2.725  _ 3 _  1  ,  1 

271  18.2  12.74  12.74 


c  =  1.0 


0.295 
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Frame  displacement 


A  _  0.295  x  11.34 
1  12.74 


0.263" 


Considering  the  horizontal  diaphragm  infinitely  stiff 
F  =  KA 
K  =  50.96 


A 


11.  34 
50 . 96 


0.223" 


Force  taken  by  the  frame 


F  =  12.74  x  0.223  =  2.84  kips 


0.251 
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STIFFNESS 
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1 

2 

3 

4 

S 

6 

7 

6 

9 

1  O 

1  1 

1  2 

1  3 

1  4 

1  5 

1  6 

1  7 

1  8 

1  9 

20 

2  1 

22 

2  3 

24 

25 

26 

27 

2  6 

2  9 

30 

3  1 

3  2 

33 

3  4 

35 

36 

3  7 

3  8 

39 

4  0 

4  1 

4  2 

4  3 

4  4 

4  5 

4  6 

4  7 

4  8 

4  9 

50 

5  1 

52 

5  3 

5  4 

S  5 

56 

5  7 

5  8 

5  9 

60 

6  1 

6  2 

6  3 

6  4 

6  5 

S  6 

6  7 

6  8 

6  9 

70 

7  1 

7  2 

73 

7  4 

75 

7  6 

7  7 


1  O 

1  02 


C 

c 

c 


1  00 

200 

300 

4  OO 


5  00 


IMPLICIT  REH>fiU-S) 

R  E  A  l  *  8  NP.NSH.NS.NSC.SS.SE.SSC.S.A.FSU.FPU.FEU.FSC.X!  100)  ,  X  1  <  1  OO  ' 
INTEGERS  NE.NF.KOUNT.IV 

’21*4,  4R  o  8  .NE.NF.NP, 

’  4  R  *  8  •  ..  SS.SE.SSC.SP) 

'  6  R  *  8  '  .A.B.FSU.FPU.FEU.FSCI 


CALL  F  R  E  A  D  (  5  , 
CALL  F  R  £  A  0  (  5  , 
CALL  FRE AO ( 5 , 

NP  «  NP  -  2  OOO 
K  OUNT  *  NE 
IFINF  CT.KOUNT) 


NSH.NS.NSC) 


K  OUN  T  x  NF 


CALL 

FREAO i 

IS,' 

R  *  8 

DO  1 

0  I  V  C  1 

,  WE 

C  1  E  * 

c  i  e  +  x  <  : 

1  V  > 

C  2  E  * 

C  2  E  ♦  X  (  ] 

!  V  )  * 

*  2 

CONT 

i  nue 

CAL  L 

FREAO i 

IS,' 

R  «  8 

DO  1 

02  1  V  = 

!  .  NF 

GIPs 

G  1  P  ♦  X  1  i 

(IV) 

C  2  P  = 

G  2  P  «  X  1  i 

1  I  V  > 

*  ®  2 

CONTINUE 

G  1  PsNSH/A*  G  1  P 

G1E*NSH/A*G1E 

XO*  <  2  .  OOO*  NS  /  SS  ♦  2  000*G1E/S6«-NP*G1P/SP  */<NSC/SSC*2  000»NS/SS«2  000 
4*NE/SE+NP*NF/SP )*A/NSH 
G2P=NSH/A/XO«G2P 
C2E=NSH/A/XO»G2E 

K  1*2  .  OOO*  NS  *  (A/NSH/KO*  1  .  OOO  ) /  SS  +  2 . 000»  IC2f«C1EI/SC  +  Nf»MG2P*G1Pt/SP 
FAILURE  AT  SEAM 

0ULTlcK1*<NS*FSU*NP*FPU-»2  000*FEU  I  /  (  A/NSM/XO*  1  0001/(2  000«NS/SS* 

8NP/SP+2  ODO/SE > 

FAILURE  AT  SIDE 

OUL  T2*K  1  *  (  NSC»FSC*NP*FPU*2  OOO*  FEU  )  /  (  NSC/SS C *NP /  S P *  2  OOO/SE  I 

FAILURE  AT  END 

I F  (  XO  .  GT  O  .  SOO  )  KOUNT*  1 

FEPVI:  (  1  .  OOO • X  ( K  OUNT  )/X0)/SE/K1 

FEPHs A/NSH/B/NE 


F  R  =  (  FEPH*  *  2  ♦  F  E  P  V  I  *  *2  )  *  *  O . 5*0 
0ULT3sFEU/(FEPH»«2*FEPVI»»2)*«0  5 


OULTsOUL  T  1 

IF10ULT2  LT  OULT)  0ULT*0ULT2 
IF (0ULT3  LT  OULT)  OULT*OULT3 


WRITE' 

6 

.  1  OO  ) 

NE 

.  NF 

,  NP  .  NS  H  , 

NS  .  NS  C 

WR  I T  E  < 

6 

.  200 

} 

SS 

,  S  E 

.  SSC  .  SP 

WR  I TE  ( 

6 

,  300 

) 

A  . 

FSU 

,  FPU  .  FEU 

,  FSC 

FORMAT 

l 

8 

'O', 

T  5 

,  ' 

NE 

',19, 

‘  f 

4 

'O', 

T  5 

.  ' 

NF 

'.19. 

'  •  '/ 

8 

'O'. 

T  5 

NP 

'  o  C20 

8  / 

8 

'O', 

T  5 

,  ' 

NS  H 

'  ,  G20 

8/ 

8 

'O'. 

T  5 

,  ' 

NS 

'  ,  G20 

8/ 

8 

'O'. 

T  5 

NSC 

'  ,  G20 

.  8  ) 

FORMAT 

< 

8 

'O'. 

TS 

,  ' 

SS 

'  ,  G20 

8/ 

8 

'O'. 

T  5 

S  E 

'  .  G20 

8/ 

8 

'O', 

T  5 

ssc 

,  G  2  0 

8  / 

8 

'O'. 

T  5 

S  p 

'  ,  C20 

8  > 

FORMAT 

( 

8 

'O'. 

T  5 

,  ' 

A 

'  .  G20 

8  / 

8 

•O', 

T  5 

,  ' 

FSU 

'  .  G20 

8/ 

8 

'O'. 

T  5 

FPU 

'  ,  G  2  0 

8  / 

8 

'O', 

T  5 

FEU 

'  ,  G  2  0 

8  / 

8 

'O', 

TS 

.  ' 

FSC 

'  ,  C20 

6  ) 

WRITE! 

6 

,400  1 

G  1 

E,G2E,CtP,C2P.XO,K1 

FORMAT 

( 

8 

'O', 

T  5 

,  ' 

G  1  E 

'  .  C20 

8/ 

8 

'O', 

TS 

,  ' 

C  2E 

'  .  G  2  O 

8/ 

8 

'O', 

T  6 

G  1  P 

'  ,  G  2  0 

8  / 

8 

'O', 

T  5 

C2P 

‘  .  C20 

8/ 

8 

'O', 

T  5 

,  ' 

XO 

,  G20 

8/ 

8 

'O', 

T  5 

K  1 

'  .  G20 

8  ) 

WR  I  T  E  ( 

6 

.  500  ) 

OULT) 

. 0ULT2 , 0ULT3 , OULT 

FORMAT 

( 

8 

'O', 

T  5 

FA  I 

LURE  AT 

THE  SEAM' 

,  G20 

8 

/ 

a 

'O', 

T  5 

F  A  I 

LURE  AT 

THE  SIDE 

,  C  20 

8 

/ 

8 

'O', 

TS 

,  ' 

F  A  I 

LURE  AT 

THE  END 

.  G20 

8 

/ 

8 

'O', 

T  5 

F  A  : 

LURE  LOAD 

,  G  2  0 

8 

» 

STOP 

END 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

1  O 
1  1 
1  2 
1  3 
t  4 
1  5 
1  6 
1  7 
1  8 

1  9 
20 

2  1 
22 

23 

24 

25 

26 

27 

28 

29 

30 

3  t 

32 

33 

34 

35 

36 

37 

38 

39 

40 

4  t 
42 


'  , A . A  1  ,D,F1  ,  K  ,  I  ,  7 , 6  I 

',HIF2,NUIF3,P! 

',SS,SE$P,SSC) 

4  R  •  8  :  '  ,  NE  ,  NF  ,  NP  ,  NSM  ,  NSC  ,  NS  ) 


IMPLICIT  REAL  • 8  ( A -$ ) 

PEAL  « •  XI  lOOO)  ,  XI  (  ioooi 
INTECER«4  NE.NF.KOUNT, IV 
CALL  FREAO IS,  '  5R«  8 
CALL  FR  E AO (  5 .  *  SR  *  8 
CAL  L  FREAO (  5  ,  '  4R ■ 8 
CALL  FREAO  <  5  .  ' 2  I «  4 
NPsNP-2. 000 
KOUNT«NE 

I  F  (  NF  .  CT  .  KOUNT  )  FOUNT  *NF 
CALL  FREA0<5.'R»8  V:  '  ,X,NEI 
00  10  I  V*  1  ,  NE 

C 1  E*C 1 E*X I  IV) 

G2E*G2E+X (  I  V  )  *  *  2 
CONTI NUE 

CALL  FREAO (  S  ,  '  R  *  8  V  :  '  ,  X  t  ,  NF  ) 

DO  102  I  V  *  1  ,  N  F 
C  1  P*C  1  P  +  X 1  (  IV) 

C2P«C2P*X  1  (  I  V  > • «  2 
CONT I NUE 
C 1 P«NSH/A*G 1 P 
G 1 E  *  NS  H / A*  G  1  E 

XOs(2.000«NS/SS*2.000*C1E/SE*NP*C1P/SP>/(NSC/SSC«2  000«NS/SS*2.000 
l*NE/SE  +  NP*NF/SP  )>A/NSH 
IFIXO.CT.O.SOOI  KOUNTsI 
G2P*NSM/A/XO*G2P 
C2E*NSH/A/X0»C2E 

K  1  *  2 . 000*  NS •  I  A/NSH/XO*  1  .  ODO  I  /SS  +  2  , OOO*  I  C2E  ’  G  1  E  >  /  SE^NP*  (  GXP  *  C 1 P  >/SP 
C  1  1  sO  14400*A«0«*4*FHK/E/T»*3/8»0 

C  1  2  *2  000«A*F2*  (  1  000-+-NU  >  *  <  1  .  000+2  .  DO»H/0  I /  0/T/E 

ODO«A«»3>F3/3  O0O/B««2/A1/E 
.  0o*A*P»SP*F3/8* *  2 
000»SP«P*  (6  .  O0O/NP  +  A*s 2  *  F  3  /  8  *  *  2  ) /A 
OOOM  1  .  ODO+  (  NSH°  1  OOO  H  ( A/ NS H/ XO "  1  ODO)  )/KI 
A , 0 , F  1  ,K , E , T  ,  8 
H1F2,NU,F3,P 


C  1  3  s2 
C  2  1  A  8  2  . 

C  2  1  6  8  2  . 

C  2  2  *2. 
WRITE! 6 , 100) 
WRITEII,  200  > 


WRITEIG  ,3001  NP.NSH.X0.K1 

WR I  TE I  6  ,  400 )  Cl  1  ,  C  1  2  .  C  1  3  ,  C2  1  A.  C2 1  B  ,  C22 
C  C  A  e  C  1  1*C12^C13^C21A  +  C22 
C  C  8  8  C  1  1+C12  +  C13«-C21B<*C22 
C  C  A  s  1  ODO/CCA 
CCB*  1  ODO/CCB 


43 

WR  I  T  E  (  6 

,  500 

> 

CCA , CCB 

44 

1  OO 

FORMAT  < 

4  S 

1  '  0  *  . 

T  5  , 

’  A  '  . 

T 

1 0 , C20 . 7/ 

4  6 

ro', 

T  5  . 

'O'. 

T 

1 O . G20 . 7/ 

47 

1  '  o  '  , 

T  5  . 

'FI' 

t 

T 1 0 , G20 . 7/ 

48 

1  '  0  '  , 

T  5  . 

'K'  , 

T 

10.620.7/ 

4  9 

1  '  0  *  , 

T  5  , 

'  E  '  . 

T 

10.620.7/ 

50 

1  '  0  •  . 

T  5  . 

'  T  '  . 

T 

10.620.7/ 

5  1 

1  '  0  '  , 

T  5  . 

'  B  '  . 

T 

10.C20.7) 

52 

200 

FORMAT ( 

53 

1  '  0  '  , 
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T 
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S  4 
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T 1 0 , G20  7 / 

55 

ro1, 
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•  NU  ' 

. 
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56 
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# 
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57 
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'  P  '  . 

T 
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58 
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FORMAT ( 
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T  5  . 
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60 

1  '  0  '  . 

T  5  . 

'  NS  H 

' 

,T10(  G  2  0 . 7/ 

6  1 

1  '  0  *  . 

T  5  , 

‘  XO  ' 

, 

T 1 0 . C20 . 7/ 

6  2 

1  0  O  '  . 

T  5  . 

'K  1  ' 

T10.C20  7  ) 

6  3 

400 

FORMAT  < 

6  4 

1  ‘  0  *  , 

T  5  , 

*  C  1  1 

.  T 1 O ,  G  2  0 . 7/ 

6  S 

1  '  0  ' 

T  5  , 

'  C  1  2 
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1  '  0  ' 

T  5  , 

'  C  1  3 

' 

.TIO.  G  2  0  7/ 

6  7 

ro*, 
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‘  C  2  1  A 

' , T 10 , G20  7/ 

6  8 

1  '  0  '  , 

T  5  , 

0  C2  1  8 

'  .  T  1 0 . 620 . 7 / 

6  9 

ro(, 

T  5  , 
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70 
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FORMAT  < 
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73 

STOP 

7  4 

END 

EncS  of  File 
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1 

IMPL ICI 

T 

REAL 

*  6  (  A  -  H  . 

,  0  -  2  ) 

2 

REAL  *  6 

L  . 

K 

3 

RE AL • 4 

X  A 

.  X  6 

,  V  A .  YB  , 

X  V  (  soo 

4 

1 

SHEAR  < 

500,5) 

5 

INTEGER 

*  4 

NO  ( 

$  ) 

6 

CALL 

FR  E AO (  -2 

.  # ENOF  I  L  E  '  . 

1  ) 

7 

CALL 

FR  E AO (  -2 

,  '  ENO  L  I  NE  ' 

0  ) 

6 

CALL 

FTNCMO  <  ' 

ASSIGN 
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*  T  I 

9 

CALL 

FRE AO ( 2 . 

'  2  R 

X  A  1  , 
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1  0 

CALL 

FR  E AD (  2  , 

'  2  R 

Y  A  1  . 

YB  1 

1  1 

CALL 

FR  E AO (  2  . 

'  2  R 

X  A2  , 

XB2 

1  2 

CALL 

FREA0I2, 

*  2  R 

Y  A  2  . 

YB2 

1  3 

CALL 

FR  E AO  (  2  , 

'  2  R 

X  A3  . 

X  B  3 

1  4 

CALL 

FR  £  A  0  (  2  . 

'  2  R 

Y  A3  . 

YB3 

1  5 

NP  * 

0 

1  6 

lO  CALL 

FRE AO ( 5 . 

'  3R  *  6  : 

'  .  L  . 

S  L  , 

1  7 

CALL 

FR  E AO (  5  , 

'  SR  *  6  : 

'  .  EF  , 

EC 

3  1 

32 

33 

34 

35 

36 

37 
36 

3  9 
40 

4  1 
4  2 
43 

4  4 
4  5 


6  6 
6  9 


60 
6  1 
62 
6  3 
6  4 


9  4 

95 
9  6 


1  1  2 
1  1  3 


V  V  (  500 , 5  )  ,  RMOM (  5  OO  ,  5  )  , 


2sT!TLES  11*1001, 


1  9 

CALL 

FREAO ( 5 

.  '  I  . 

4  R  * 

20 

CALL 

FREAO ( S 

'  4  R 

*  6 

2  1 

CALL 

FR  E AD (  S 

.  '  3  R 

•  6  : 

22 

CALL 

FREAO  <  5 

,  '  R  * 

6  :  ' 

23 

NP 

s 

NP  ♦  1 

24 

AS 

z 
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25 

SW 

B 

T  /  (  4 
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H/L 

26 

S  B 

£ 

8  1  /  L 

27 
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8 

Cl  /  H 

26 

S  8 
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EF  • 

S  B 

29 
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S 
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S 

30 
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GW  ) 
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SBL  *  OSORT <  B*K«CF/2  .  OOO/CC/AS  > 
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DL*SD*SL 

I  F  <  S  B  L  -  C  )  12,13,11 

SF  *OSORT (  (SBL**2‘C**2)/2  OOO) 

FL  *  SF*  SL 

P  2  *  OS  I NH <  <  0 l ♦ F L  )  / 2  000  > 

T  2  *  DCOSH (  (  OL  • F  L  )/2  OOO) 

RH2*  OCOSM (  (OL*FL  1/2 .0001 
VI*  OS  I NH <  ( 01  - F l  )  /2  OOO  ) 

R  2  *  S  0  ♦  S  F 
RN 1  *  SO  -  SF 

A2*  SW*HL*R2**2*T2/(  (R2«*2*T2*  (R2**3«R2»2  .  ODOsm®  *  2«R2*P2 
I - RM  2  *  S W/ E  C / S 1 » > 


4  6 

l*RN1**2*RM2*  ( RN  1  *  »  3*  V  1  -  2  O0O*S8L**2BRN 

1 

*  V  1  * T2°SW/IC/S  I 

)  >  / 

4  7 

6  2 

ODO/GF/EC/S 1 

46 

B  2 

«-A2*RN1**2«RM2/R2**2*T2 

4  9 

GO 

TO  400 

SO 

1  2 

SF 

*  OSOR  T (  <C**2  -  SBL**2)/2  OOO) 

5  1 

F  L 

*  SF  *  SL 

52 

S  1 

s  OTAN (  F L /2  .  OOO  )  •  0 T A N M (  0 L / 2  0 0 0  > 

S3 

V 

*  OS  1  N  (  F L / 2  .  OOO  >  •  OS  I  NH (  0 L / 2  .  OOO  ) 

54 

T  1 

*  OCOSlFL/2  ODO)  *  0 C 0 S H ( 0 L / 2  0 0 O > 

55 

P  1 

s  OS  I  N (  F L /2  OOO )  *  OCOSH (  OL/2  OOO  > 

56 

RM 

s  DCOS  (  F L /2  .  OOO  >  *  OS  I  NH (  0 L / 2  .  OOO  1 

57 

R  1 

8  S 0  *  *  2  •  S F  *  *  2  •  2.000  *  SBL  * 

* 

2 

5  8 

0  1 

8  SO  **  3  -  3.000  *  SF  **  2  *  SO 

2 

ODO  *  SBL 

*  «  2 

*  SO 

S  9 

R  N 

*  SF  **  3  •  3.000  *  SO  •*  2  *  SF  ♦ 

2 

OOO  »  SBL 

*  8  2 

8  S  F 

6  0 

H  1 

s  2.000  •  SF  *  SO  *  ( P 1 *RN  ♦  RM*01 

♦ 

S  1  *  P  1  ■  0  1  ■ 

.  s  1  *  RM*  R  N  « 

SW8 

6  1 

1  T  1 

/EC/SI  •  S1*V*S W/EC/SI)  ♦  R1  *  <S1*P 

1 

*RN  *  s  1  8  R  M  *  0 1 

-  P  1  *  0  1 

•> 

6  2 

2RM 

«  R  N  ) 

6  3 

A  1 

8  SW  •  HL  «  ( 2  ODO  ®  S  F  *  S  0  ♦  S1*R1>  / 

2  OOO  /  GF 

/  EC 

/si/ 

H  1 

6  4 

B  1 

*  SW  *  HL  *  (  2  .  OOO* S F * S 0 • S  1  -  R1)  / 

2  OO©  /  GF 

/  EC 

/SI  / 

H  1 

6  5 

GO 

TO  400 

6  6 

1  3 

SF 

s  OSORT (  (  C  *  *  2  -  S B L *  *  2  )  /2  .  ODO  » 

6  7 

P  3 

s  OS  I  NH  (  OL / 2  OOO  ) 

RM3=DC0SH  <  OL / 7 
T3=  S0**2*RH3 


70 

V  2  8  2  OOO*  S  0  *  R M 3  - 

SO**  2  *  S  L / 2 

O 

Q 

O 

«  P  3 

7  1 

R3sP3*SD**3*RM3« 

SW/EC/S I 

72 

R  N  2  s  -SD**2*SL 

*  R  M  3 / 2  .  OOO  +  SO 

8*2* 

P  3  • 

P  3 

73 

A  3  8  H  L  «SW* V2/2 

ODO/CF/EC/S I 

/  (  T3 

*  R  N2  * 

7  4 

B  3  8  A  3  *  T  3  /  V  2 

75 

400  CONT I NUE 

7  6 

WRITE  (6.20) 

L  , 

S  L  ,  H  ,  EF  , 

EC  , 

CC 

77 

IS.  S  N 

7  8 

20  FORMAT  (  '  1  '  , 

TS  , 

'Width  of 

Bu  i 

1  d  i 

"0 

7  9 

1  C  20  .  7  , 

/ T  5  ,  'Length 

o  f 

8  u  i 

1  d 

'  SW/ EC/SI 


EW  , 


GW  ,  I! 


C  I  . 


S  l  , 


S  8  . 


I  in  I 

ng  (in) 

G207/T5,  ' Height  of  Building  (in) 

G  2  0  7 / T  5  ,  'Young'  's  Modulus  for  Frame  (Kips/in»«2i 
C20.7/T5.  'Young'  's  Modulus  for  Concrete  (K ips  /  in*®2 
G20.7/T5,  Shear  Modulus  for  Concrete  (K  ips  /  in«  »2  I 
C20  7 / T  5  .  'Young'  ' s  Modulus  for  the  wall  (K  ips/  in»«2 

(Kips  /  i  n»  »  2  1 


65 

7 

G  2  0 

7/TS  , 

'Shea 

r  Modulus  for  the 

wa  1  1 

6  6 

6 

C  20 

7  /  T  5  . 

'I  of 

Beam  (  in*»4  i 

6  7 

9 

G  2  0 

7/T5  . 

'  1  of 

Col umn  «  i  n  *  »  4  i 

8  6 

• 

G  2  0 

7/T5  . 

'  I  of 

S  1  ad  (  i  n  *  »  4  ) 

8  9 

1 

C  2  0 

7/TS  . 

'  S  t  i  f 

f ness  o  f  Beam  ( 

i  n 

«  »  3  ) 

90 

2 

G  2  0 

7  /  T  5  . 

'St  i  f 

fness  of  Column 

< 

i  n  •  «  3  ) 

9  1 

3 

G  20 

7  ,  /TS 

'  S  t 

i  f  f  ness  o  f  w a l  1 

i 

i  n  *  B  3  i 

9  2 

4 

G  20 

7  /  T  5  . 

'  S  t  i  f 

fness  of  Fr  ame 

(  K 

i  p  /  i  n  ) 

9  3 

5 

G  2  0 

7  /  T  5  , 

’  T  o  t  a 

1  S  t  i  f  fness  o  f 

F  r 

ame  (  K 

WRITE 

FORMAT 


G  2  0  7 

6,30) 
(  /T5  , 


AS  . 
Of  F  f 


0  ,  B 
a  me  s 


9  7 

1 

G  20 

•  7  /  T  5  , 

Spacing  Between  F 

9  6 

2 

C  2  O 

.  7  /  T  5  . 

'Thickness  of  Wall 

9  9 

3 

G  20 

7/T5  , 

'Area  of  Slab  (in) 

1  OO 

4 

C  20 

7  /  T  5  , 

'Thickness  of  Slap 

1  O  1 

5 

G  2  0 

7  /  T  5  , 

'Beta 

1  02 

6 

C  20 

.7  /  T  5  . 

'  K  appa 

1  03 

7 

C  20 

7  /  T  5  , 

'Hor  iiontal  Load  (1 

1  O  4 

6 

G  2  0 

7  > 

1  05 

WRITE 

i6,50»  G  F  , 

SBL,  C.  SO.  SF 

1  06 

SO  FORMAT 

(  '  0 

'  .  T  S  , 

'UNIT  ST  I  FFNESS  <  K 

1  07 

1 

C  20 

7/T5  . 

'  b  ( d i mens i ons 

1  O  8 

2 

C  20 

7  /  T  5  , 

' c  Iwird) 

1  09 

3 

G  2  0 

7/T5  . 

' Sma 1 1  d 

1  1  0 

4 

G  2  0 

7/TS  . 

' Sma 1 1  f 

1  1  1 

5 

G  20 

7  ) 

NOl  NP  ) 

X  =  X  I 


236 


End 


I  1  4 
I  1  S 


I  27 
i  28 


I  4  1 
I  42 
I  43 
I  44 
I  4  5 
I  46 
I  4  7 
1  4  8 
I  4  9 
I  SO 
1  5  1 
I  52 
!  5  3 
I  54 
1  5  5 
I  56 
I  57 
I  5  6 
I  59 
I  60 
I  6  1 
I  6  2 
!  6  3 

I  6  4 
o  f 


4  9  9 

500 


IMS6L-CI  500 . 502 . 50  1 

Y  *  ML  /  CF  «  (  1  .  OOO  +  SW/M 1 /EC/S  l  •  I  I  2  ODO«SF«SD 


1  6 

1  ■ OCOS  <SF«X) 

1  7 

2*  OCOSH  <  SO*  X  ) 

♦  <2 

OOO*  S  F 

«S0»  S  1  -  R  1  1  8 

1  8 

1  9 

80  FORMAT  <  T  5  , 

5020 . 7 

> 

20 

NO  (  N  P  1  s  NO  < 

NP  1  ♦ 

1 

2  1 

X  V  1  NO (  NP  1  ,NP 

)  *  Y 

22 

YVINOINP)  ,  N  P 

)  *  X 

23 

Y 1  *2  OOO  • 

OS  I  N  < 

SF»X  1 

*  OCOSH  <  SO*  X  ) 

24 

1 OCOS  <  SF  *  X  1  • 

OS  1  NH 

(  SO*  X  1 

«  <  A  1  •  S  0  ♦  8 

25 

Y  2  *  2.000  ■ 

OCOSI 

S  F  »  X  1 

•  OCOSH  <  SO*  X  > 

26 

1SF  -  A  1  ■ S  F  * • 

2  )  ♦  2 

OOO  • 

OS  I  N 1  SF  *  X  )  • 

#  OS  I  KM (  SO®  K  )  >  ) 


2000*  A  1  • SF  •  S  D 
Y  3  *  2  .  OOO  • 


-  8  1  • SF • • 2  ) 

OS  !N(SF«X  I  > 


OCOSH (  SO*  X  ) 


<  B  1  • SO 
8  1  *  SF  1 

*  ( A  I  ®  SO®  >  2 

DSIMM(SO®X) 

s  (81®S0®«3 


A  1  *  s  F  )  *  2  .  OOO 


2000*81*50* 

(  8  1  •  S  0  *  *  2  -  2 


3  OOO®  A  I  »SF  ®  SO 


29 

1  *  •  2  - 

3  OOO®  8  1  • SO 

•SF* 

•  2 

♦ 

A  1  •  SF  ® 

•  3  1  ♦  2  OOO 

8  OCOS  <  SF  *  X  1  •  OS  I  NH 

30 

2  <  S  0  8  X 

1  •  <  A  1  •  S  0  *  « 

3  ♦ 

3  . 

OOO*  8  1  *  SF 

•  SO*  •  2  *  3 

OOO « A  1  * S 0 « S F * • 2  •  Bis 

3  1 

3  S  F  *  «  3 

1 

32 

CO  TO 

504 

33 

50  1 

Y*  2  . 

OOO* A2*  OCOSH 

<  R  2* 

X  1 

♦  2 

000*82 

• OCOSH  <  RN  1  • 

X  1  *HL/GF 

34 

WRITE 

<6.601  X  .  Y 

35 

CO  TO 

SOS 

36 

502 

Y*  2  . 

OOO*  A3 • OCOSH 

<  S  0* 

X  1 

♦  2 

OOO*  83 

*  OS  I  NH <  SO®  X 

1  *  X+HL /GF 

37 

WRITE 

<6.801  X .  Y 

36 

CO  TO 

505 

39 

504 

WRITE 

<6.801  Y 1 . 

Y  2  . 

Y  3 

40 

RMOM ( NO  I  NP  )  ,  NP  1  * 

-  EC 

8 

S  1 

*  <  Y2 

•  2  OOO* SSL  *  *  2* Y  ♦  2  .  OOO* SBL • • 2* 

1HL/CF) 

SHEAR ( NO ( NP I 


NP  ) 


WRITE  <6,601  X, 

>5  X  *  X  ♦  H  I  NCR 

IF  <  X  IE  XF  >  CO  TO  499 

CO  TO  10 
90  CONTI  HUE 
STOP 

1999  00  100  IV  *  1.  NP 

CALL  PLOTITIXVI  1,  IV),  YV<  I 
1  5  0.  Y  A  1  ,  YB1,  9.0.  1) 

lOO  CONTINUE 

DO  110  IV  *  1.  NP 

CALL  PL OT  I  T  < SHEAR  <  1  ,  I  V  )  . 

1  XB2  .  5.0,  Y  A  2  .  YB2  .  9  O . 

1 1 O  CONTINUE 

00  120  IV  *  1 .  NF 

CALL  PLOTITIRMOMI  1.IV)  , 

1  X  8  3  ,  5  O  .  Y  A  3  ,  Y  8  3  ,  9 

120  CONTINUE 

CALL  PLOT  1  T  (XV.  YV,  ND( 

1  ,  9.0,  1  ) 

STOP 
ENO 


-  EC  «  SI  *  <  Y3  - 

SHEAR  <  NO  <  NP  )  ,  NP  ) 


2  .  OOO®  SB l *  ®  2®  Y  1  1 
,  RMOH <  NO  <  NP  1  .  NP  > 


Y  V  <  1  , 
1  1 


YVI  1  ,  IV)  ,  NO  <  I  V  )  , 


1 

IMPLICIT  R  E  A  L  •  8 

<  A  -  H 

,  0  *  Z  > 

2 

REAL 

•4  L  ,  K 

3 

CALL 

FREADIS  ,  '  4R 

«  8  ' 

•L.SLI. 

SL2 , Ml 

4 

CALL 

F  R  E  A  0  (  5  ,  '  5R 

•  8  : 

,  EF  .  EC  . 

GC . EW. 

GW  1 

5 

C 

NOTE 

*  W  =  >  HL  ( 

MOR  I 2  ONT AL 

LOAD) 

6 

CALL 

FRE AO  <  5  ,  '  1  , 

4R  8  t 

:  '  ,  N  .  E  , 

B  I  .  C  I  . 

ML  ) 

7 

CAL  L 

FR  E AD  <  5  ,  '  4R 

•  8  : 

.  0  .  T  ,  e 

K  1 

8 

CALL 

FREAO (  5  .  ' R  » 

8  '  , 

S  1  > 

9 

AS  * 

0  •  L 

1  0 

sw  * 

T  /  <4.000* 

<  H/  L 

>  •  8  3/EW 

♦  H/GW/L 

1  1 

SB  « 

B  I  /  L 

1  2 

SC  8 

Cl  /  M 

1  3 

S  8 

1 2 . OOO  «  EF 

•  SB 

/  H  8  * 

2  /  < 

1  Ol 

1  4 

SN  8 

OF  L  OAT  <  N  )  « 

S 

1  5 

GF  8 

SM  / <  S  L  1 *Sl 2  > 

1  6 

S  B  L 

S  OSORT(6«K« 

GF/2 

. OOO/CC/AS 1 

1  7 

C  * 

(CF/EC/SI)  • 

*  0 

2S  OO 

1  8 

SO  « 

OSORT  <  <  C« • 2 

*  S8L*e2)/2  OOO ) 

1  8 

Z  *  H  L  *  S W/  (  CF*EC*SI  ) 

2© 
2  1 
22 
23 


4  7 

4  t 

4  9 
50 

5  t 

52 

53 
5  4 
5  5 
56 
S  7 
56 

5  9 
60 

6  1 
6  2 
6  3 
6  4 
6  5 
6  6 
6  7 
6  6 

6  9 
70 

7  1 

72 

73 
7  4 

75 

76 
7  7 
7  6 

7  9 
80 

8  1 
8  2 
8  3 
8  4 


9  3 
9  4 
9  5 
9  6 
9  7 
9  8 
9  9 
1  00 


1  09 
1  1  0 
1  1  1 
1  1  2 
1  1  3 


000*2  O0O*SB/SC > 


0  l  USD'SL  1 
0L2*S0*SL2 

WRITE(«  ,2001  l.SLI  ,SL2.M.EF,EC.GC.EW,GW.8I  ,CI 
FORMAT (  '  1  ‘  , T5 ,  'Width  of  Building  (  in  ) 


,SB,SC,SW,S,SN 

'  ,  G20  7  . 


24 

4/  T5,  'Length  of  Fi 

PSt 

span 

' . G20  7  /  T 

5  .  ' 

L  eng  h 

25 

It  of  second  span 

' . G20  7 

/  TS  , 

'Height 

o  f 

du  i  1 

26 

4 d  i  ng 

'  ,  C20  7  / 

TS,  'Young'  's 

Mo  du 1 u  s 

for  F  r  a 

27 

4r»«  *  , 

G20  . 

7  / 

T5 , ' E  for 

Cone  ret© 

2  8 

4  * , G20  7  /  T5 , 

' 

Shear  E 

for  concrete 

•  .620 

28 

4  7  /  T  5  ,  '  E  f Of  w«  1 

1 

'  ,  G20 

7  / 

TS  ,  ' 

30 

lsh«ar  e  Fop  1  1 

'  .  C  20 

7  / 

T  S  ,  ’  I  f 

o  r 

beam 

3  1 

4 

* . C20  7 

/  T  5  .  '  I 

of  CO 

1  uran 

32 

4 

G  20 

.  7  / 

T  5  ,  '  I  OF 

s  1  ad 

33 

4  ' . C2©  7  /  T5 

t 

'  S  t 

i  f  fness  o  f  deant 

.  G  2  © 

34 

4.7. 

35 

4  /  T5,  'Stif fness  ©  f 

c  o 1 umn 

C  2  0  7  / 

T  5 

.  'Stiff 

36 

Iness  o  f  hi  1  1 

'  ,  C20  7 

/  T  5  ! 

'Stiff  ness 

©  f 

37 

4  ©  f  F  p  am©  <  i  n  *  *  3 ) 

' 

,6207/ 

38 

4  4  total  stiffness 

o 

f  f  p  ame 

'  ,  G20 . 7  ) 

39 

WR  IYE<6 . 400 1 N  ,  E,  T. 

AS  ,  0 

.  8  .  K 

.  H  L 

4© 

400 

FORMAT ( /  T  S  .  ' Number 

o  f 

F r am©  s 

' 

4  1 

£  5 . G20 

7 

/  TS 

'  S p ac  i  ng 

Bet  we©n 

F  p am©  s  <  i n  1 

4  2 

4  *  .  G  2  0  7  /  T 

4  3 

45,  'Thickness  ©f  W* 1 

1  < 

i  n  1 

'  .  G20 

.  7  /  TS 

.  'Area  ©  f 

4  4 

4  S  1  ad  (  in) 

'  , 

G20.7  /  TS.'Th 

i  ckness 

o  f 

S  l  ad 

4  5 

4 

'.620.7  /  T  5 , 

'Beta 

4  6 

6  "  ,  G 

4  2  0.7  /  T  5  ,  ‘K«ppi> 

I'Hor  iiontil  Lo«d  (hips/in) 
l F ( SBl - C  >  t  2  ,  13,11 
SFsOSORTI (C8»2*SBL®82)/2. 000) 
FD2*2*SF«S0 
FL2*SF*SL2 

R * S 0  * s 2  -  S F *  *  2 • 2  O0O*SBL8«2 
FI  1  *  S  F  *  S  L  1 

S  1  *  0  E  X  P (DL  1  I*0SIN(FL1  ) 

T 1  *  0  E  X  P  (  OL  1  >  » DCOS  <  FL  1  > 

U  1  s  0  E  X  P  <  *  0  L  1  I  *  OS  1  N  (  F  L  1  > 

V 1 s  OEXP  <  • OL  1  I  • 0  C  OS  (  F  L  1  1 
S2*0EXP(DL2l*0SIN(FL2) 
T2*0EXP!0L2t«0C0S(FL2l 
U  2  *  0  E  X  P (  -0L2>*0SIN<FL2> 

V  2  *  0  E  X  P <  °  0  L  2  >  *  OCOS  ( FL2  > 
P1*T»«R»2  00O»SF*SD«S1 
OOO*  S  F *  S  D  8  T  1 
OOO*  SF  * S  0*  U 1 
ODO*SF»SO*V1 
O0O*SF*S0*S2 
O0O*SF«S08 T  2 
OOO®  SF»S0*U2 
ODO*SF*SDI  V  2 


,  G20  7  /  T  5  „ 


OlsSI»R*2 
R  1  *  V  I  *  R  *  2 
WUll  I  *  R  °  2 
P2*  T2*R  •  2 
02=S2*R*2 
R2* V  2  »  R  ♦  2 
W2sU2*R»2 


E  1  *  X  1  • 
F  1  *  Y  1  ■ 
G  t  *  X  1  • 
H  1  *SW, 
T  3  =  S  0  - 


*  3  • 

3  OOO 

V 

SO* 

S  F 

8 

8 

2 

-  2  ODO* 

S  8  L  *  *  2 

« 

S  0 

d  3  - 

3  000 

a 

SO  8 

a 

2 

a 

SF 

♦  2  ODO« 

S  8  L  «  *  2 

8 

SF 

S0» 

GC  *  AS 

/ 

< 

EC 

a 

S 

I 

8 

B 

*  K  > 

S  F  8 

GC  «  AS 

/ 

< 

EC 

8 

s 

I 

a 

B 

»  K  > 

\ 

SW/ 

<  EC  «  S 

I 

) 

<  EC 

•SI  1- 

X 

1 

«  3  * 

T  1  ♦  S  F 

■ 

8 

3  » 

s 

1 

- 

3 

0008 

so* 

*  2 

a  S  F  * 

s 

1  - 

3 

OOO  » 

S  D  • 

SF  » 

«  2  »  T  1 

»  3  8 

St-SF 

a 

a 

3  * 

T 

1 

4 

3 

0008 

SO  * 

»  2 

a  SF  - 

T 

1  - 

3 

ODO  • 

SO* 

SF* 

a  2  »  S  1 

8  3  8 

V  1  -  S  F 

a 

a 

3  « 

U 

1 

♦ 

3 

0  0  0  8 

so® 

«  2 

«  S  F  a 

U 

1  - 

3 

OOO  * 

S  0  8 

SF  a 

a  2  *  V  1 

•  3  “ 

U  1  +  SF 

8 

a 

3  8 

V 

1 

• 

3 

000« 

SD*SF 

*  »  2  » 

u 

1  * 

3 

ODO" 

so* 

»  2  * 

SF«V1 

T  1  - 

SF»S1 

S  1  * 

SF  «  T  1 

V  1  ♦ 

SF  *  U  1 

V  1  • 

S  0  *  U  1 

85 

T  5 

s 

T  3 

- 

2 

OOO 

a 

S  B  L 

* 

8 

2 

a  T  4 

- 

T  1  *  SW/  ( 

EC* 

S  I  > 

8  6 

S  5 

= 

S  3 

- 

2 

ODO 

a 

S  B  L 

« 

8 

2 

*  S  4 

• 

S  1  *  SW/  ( 

EC* 

S  I  > 

87 

V  5 

s 

V  3 

- 

2 

OOO 

a 

S  B  L 

a 

8 

2 

8  V  4 

♦ 

V  1  8  SW/  ( 

EC  » 

S  1  1 

8  6 

US 

= 

U  3 

♦ 

2 

OOO 

a 

S  B  L 

a 

8 

2 

a  U  4 

*• 

U  1  *  SW/  ( 

EC* 

S  I  1 

8  9 

T  6 

r 

S  0 

a 

a 

3 

«  T  2 

♦ 

SF  » 

a 

3 

8 

S  2  - 

3 

OOO  *  S  0 

a  a  2 

•  S  F  *  S  2  ' 

3 

OOO  a 

S  0  *  S  F  a 

-  2  * 

T  2 

9  0 

S  6 

c 

S  0 

a 

m 

3 

-  S  2 

- 

S  F  » 

a 

3 

a 

T  2* 

3 

ODO •  S  D 

a  a  2 

*  S  F  *  T  2  • 

■  3 

OOO* 

SO«SF« 

»  2  * 

S  2 

9  1 

V  6 

s 

S  0 

a 

a 

3 

a  V  2 

- 

S  F  » 

o 

3 

a 

U2  + 

3 

ODO  »  S  D 

»  •  2 

*  S  F  •  U  2  • 

■  3 

ODO*SD»SF» 

«  2  « 

V  2 

9  2 

U  6 

= 

SO 

a 

a 

3 

*  U  2 

♦ 

SF  - 

» 

3 

a 

V  2  * 

3 

OOO* SO 

•  SF 

*  a  2  *  U2  • 

3 

OOO* 

S  0  «  •  2  » 

S  F  • 

V  2 

T8*S0*T2-SF  »  S  2 
SA  =  S0*S2*SF*  T  2 

V  8  8  S  0  * V2*SF»U2 
U8sSF*  V2  •  SO*  U2 

T  9  2  T  6  °  2  O0O»SBL*»2*T8T2*SW/<EC»S!  > 

S  9  =  S  6  -  2  O0O*SBL**2*S8-S2*SW/  (  EC*S  I  ) 

V  9  =  V  6  °  2  O0O*SBL*  *2* V8*V2*SW/  I  EC'S  I  > 

U9:U6«2  ODO*SBL • »2*U8*U2*SW/ I EC»S I 1 


1  0  1 

WR  I  T  E  <  6 

.223  1  1 

T5  .  S  5  ,  VS 

1  02 

223  1 

FORMAT ( 

T  5  =  • 

.  G  20  7 /  ' 

1  03 

WR  I  T  E  (  6  , 

,  2  2  3  2  ) 

P  1  ,  0  1  .  R  1 

1  04 

2232 

F  ORMA  T  < 

Pis' 

.6207/' 

1  05 

WR  I  T  E  «  6  , 

,  2233  1 

G  1  .  X  1  ,  Y  1 

1  06 

2233 

F  ORMAT  <  ‘ 

'  G  1  =  ' 

.  C  20  7  /  ' 

1  07 

WR1TEI6  , 

,  2234  1 

E  1  .  F  1 

1  08 

2234 

F  ORMA  T  <  1 

E1=' 

, C20  7 / 

WR  I  T  E  (  6  ,  2  2  3  5  )  R.FD2 

FORMAT  <  ’  R  =  ’  ,  G  2  O  7/'  F 

WR  I  T  E  (  6  ,  2  2  3  7  >  T9,S9,V9 

FORMAT ('  T9* ' , C20  7/ ' 

WR  I  TE  (  6  ,  2238  )  P2.02.R2.W2 


S  5  s 

,  G20 

7/  ' 

VS  =  ' 

,  G20 

7/  ' 

US  *  ' 

G  2  0 

7/ 

.  W  1 

0  1=' 

,  C  20 

7/  ' 

R  1  s  ' 

.  G20 

7  /  ' 

Wl:1 

.  G  2  0 

7  / 

,  H  1 

IMS' 

,  G  20 

7/  ' 

Y  1  *  ' 

.  C  2  0 

7/  ' 

HI:’ 

.  C  20 

7/ 

Fir' 

.  C  20 

7/  1 

02-  • 

.  G  2  0 

7  /  . 

i  .  U  9 

S  9  = 

.  G  2  O 

7  /  ' 

V  9  s  • 

.  C  20 

7  /  ' 

U9  =  ' 

.  G  20 

7  / 

1  1 4 

1  1  5 

1  1  6 

1  t  7 

i  i  e 

1  1  9 

1  20 

1  2  1 

1  22 

1  23 

1  2  4 

1  25 

1  26 

1  27 

1  28 

1  29 

1  30 

1  3  1 

1  32 

1  33 

1  34 

135 

1  36 

1  37 

1  38 

1  39 

1  40 

1  4  1 

1  42 

1  4  3 

1  44 

1  4  S 

1  46 

1  47 

1  4  6 

1  49 

1  SO 

1  5  1 

1  5  2 

1  S3 

1  54 

1  55 

1  56 

1  57 

1  5  8 

1  S  9 

1  6  O 

1  6  1 

1  6  2 

1  6  3 

1  6  4 

1  6  5 

1  6  6 

?  6  7 

1  6  6 

1  6  9 

1  70 

1  7  1 

1  72 

1  73 

1  74 

1  75 

1  76 

1  77 

1  7  8 

1  79 

1  80 

1  8  1 

1  82 

1  83 

1  8  4 

1  85 

1  86 

1  8  7 

1  8  8 

1  8  9 

1  9  O 

1  9  1 

1  9  2 

1  9  3 

1  9  4 

1  9  5 

1  9  6 

1  9  7 

1  9  8 

1  9  9 

200 

20  1 

202 

203 

204 

205 

206 

207 

208 

209 

2  1  O 

2  1  1 

2  1  2 

2  1  3 

2  1  4 

2  1  5 

2  1  6 

2  1  7 

2  1  8 

2  1  9 

2  20 

2  2  1 

2  22 

2  2  3 

2  2  4 

2  25 

2  26 
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2238  FORMAT ( *  P2« ’ ( G20  7/4 

WRITE(6,2239)  2 

2239  FORMAT ( '  Z«  '.G2Q  7/) 

GO  TO  500 

13  SPiO 

S2«0EXP  <  DL  1  > 

T  2  *  0  E  X  P <  «  OL  1  ) 

S  3  *  0  E  X  P  <  012  ) 

T  3  *  0  E  X  P  (  -  0L2  ) 

R  1  1  a  -SD«*3*S2“S2»SW/EC/S  I 

SI  1tSDt*2*S2*S0*O«Sl  1«S2-S2»SL  1*SW/|C/SI 
T 1 1»S0»»3#T2-T2*SW/EC/S! 

U  1  1*SD«»2*T2*SD*«3*SL  1»T2°T2«SL  1®SW/EC/SI 

R  1  3  a  S  0  » »2«S2 

S  1  3«S0«S2»  (OU-21 

T  1 3  a  S  0  s a  2»  T2 

U13*S0*T2»I0LH2I 

R  1  4aS0«  «  3  - SW/EC/S  I 

S  1  4*S0« • 3 

T14*S0«*2 

U 1 4  «S0*  *  3  +  SW/EC/S I 
R  1  5«S0*  *  2 
S  1  5  a  S  0  •  (  0  L  1  -  2  > 

T  1  S*S0«  <  012  -  2  ) 

UlSsSO*  t  DL  1  +  2  ) 

R 18«SD*B*K»EC»SI/CC/AS»S08*3 
S  1 6  a  1  ,O0O  +  B*K»EC»SI/CC/AS*S0**2 
R  2 1  * *S0«»3*S3-S3«SW/EC/S I 

S21*S0«*2»S3-S0»* 3«SL2»S3-S3»S12«SW/EC/S 1 
T21*SD*«3«T3-T3*SW/EC/SJ 

U2 1 sS0*»2*T3+SD*«3*SL2*T3* T  3  <  S  L  2*  SW/ I  C/S  I 

R  3  1  *  S  0  »  •  2  *  S  3 

S3  1  «S0aS3*  <  DL  2  -  2  > 

T3  1  *S0®  *  2*  T3 
U31*S0»T3*  (  0  12*2  ) 


02«',G20  7/'  R2*'(C20 .7/'  W2*4,G20  7/) 


WRITEI6 , 

,1231) 

R  1  1  .  S  1  1  , 

,  T  1  1  .  U  1  1 

123  1 

FORMAT  <  ‘ 

1  /  > 

R  1  1  a 

4  , C  2  0  7/ 

'  S  1  1  *  4  . 

G  2  0 

7/  4 

T11*‘ 

,  G  2  0 

7/  ' 

U  1  1  *  4 

,  G20 . 7 

WR  I T  E (  6  , 

,  1  232) 

R13.S13, 

.  T  1  3  .  U  1  3 

1232 

FORMAT ( 1 

1  /  > 

’  R  1  3* 

*  ,  G  2  O  7/ 

'  S  1  3  =  '  . 

G20 

7/  4 

T  1  3  a  4 

.  G  2  0 

7/  4 

U  1  3  * 

.  G  2  0  7 

WR1TEIS, 

,  1  233) 

R  1  4  ,  S  1  4  , 

.  T  1  4  .  U  1  4 

1233 

FORMAT  <  ‘ 

1  /  ) 

R  1  4  a 

.620  7  / 

■  S  1  4a  *  . 

G  2  O 

7/  4 

THs' 

.62© 

7/  4 

U  1  4  a  • 

.  G  2  O  7 

WR  I  T  E  <  6  . 

,  1234) 

R15.S1S. 

,  T  1  5  ,  U  1  5 

1234 

FORMAT ( ‘ 

1  /  > 

’  R  1  5  * 

4  .620.7/' 

’  S  1  5*  '  . 

G  2  O 

7/  ‘ 

T  1  5  a  ' 

,  G20 

7/  4 

U  1  S  *  ' 

,620.7 

WR I  T  E  (  6  . 

1  23S  ) 

R  1  6  .  S  1  6 

1  2  3  S 

FORMAT ( ‘ 

R  1  6  a 

'  .  G  2  0 . 7/  ‘ 

S1l«' 

G  2  O 

7/  ) 

WR  I T  E (  6  . 

1  236) 

R  2  1  .  S2  1  . 

,  T2  1  ,  U2  1 

1236 

FORMAT ( * 

1  /  > 

WRITEI6, 

R  2  1  a 

* , C20 . 7/ 1 

'  S  2  1  *  '  . 

G  2  0 

7/  4 

T  2  1  *  4 

.  G20 

7/  4 

U2  1  *  4 

,620.7 

1  237) 

R3  1  .  S3  1  . 

,  T3  1  .  U3  1 

1237 

FORMAT ( ‘ 

1  /  > 

R  3  1  a 

'  .  C  2  0 . 7/  • 

S31e', 

G  2  0 

7/  4 

T3  1  *  ' 

.  G  2  0 

7/  4 

U  3  1  *  ' 

, G20 . 7 

WR  1 T  E (  6  . 

1  238) 

2 

1238 

FORMAT ( * 

2  *  '  . 

G20 . 7/  ) 

2  13  1 
2  132 
2  133 

2  135 
2  136 
2  13  7 
2  13  8 
2  13  9 


GO  TO  500 

I  SF * DSOR T  (  < S 8 L * « 2  -  C « *  2  > / 2  00©  > 

FL  1  *  SF  «S  L  1 
FL2sSF*SL2 

R  1  la  (SD^SF>**2»2  OOOaS8L**2 
R22*  <  SO-SF  )  *  *  2  -  2  000«SBL«*2 
S22eOEXP  (  0 t  1  ♦ F L  1  ) 

T22=0EXP  <  - DL  1  -FL  1  ) 

U22*0EXP  <  “OL  1  I  *  OS  I N (  F  L  1  ) 

V22*DEXP  <  -DL  1 ♦FL  1  ) 

S33*0EXP  (  0  L  2 ♦ F  L  2  ) 

T33*0EXP<  -DL2-FL2  > 

V33*0EXP (  • OL  2  +  F  L  2  > 

U33*0EXP(  DL2-FL22  t 

XI  1*  (  S  D  ♦  S  F  )  »  *  3  -  2  ODO*SBL**2*  (  S  0  *  S  F  ) 
X  2  2  a  (SO-SF  >**3-2  .O0O»SBL*'2*  (  SO-SF  ) 
G  1  1  a  X  1  1  ♦SW/  (  E  C  *  S  I  > 

Mils  SW/  <  EC  «  S  I  >  -  X  1  1 
G  2  2  a  X  2  2  ♦  S W/  (  E  C  *  S  I  > 
C33=S0+SF*B«K*EC*SI«X1 1/GC/AS 
M  2  2  *  S W/  (  EC*  S  I  >  -  X22 
H33*S0-SF+B*K«EC*S]*X22/CC/AS 
R4:S22*R  1  1 
T4*T22*R  1  1 
U4sU22*R22 
V4=  V22«R22 
R5  =  S33aR  1  1 
T  5  =  T  3  3  »  R  1  1 
U5=U33*R22 
V5aV33«R22 
S6=S22*H1 1 
T6eT22*G1  1 
U6sU22»H22 
V6=V22*G22 
S7=S33«H1 1 
T7e  T33  «  G  1  1 
U7=U33»H22 
V  7  a  V  3  3  •  G  2  2 
WR  I  T  E  <  6  ,  2  1  3  1  > 

FORMAT  (  '  S  6  *  ' 

WRITE(6.2132> 

FORMAT  (  '  R  4  *  ' 

WR  I  TE  I  6  ,  2  1  33  I 
FORMAT  (  '  Gila 
1/'  X  1  1  *  '  ,  C20 
WR  I  T  E  <  6  ,  2  1  3  5  > 

FORMAT  <  4  R  1  1  a 
WR  I  T  E  <  6  ,  2  1  3  6  > 

FORMAT  (  '  G  3  3  a 
WR  I  T  E  <  6  ,  2  1  3  7  ) 

F  ORMA  T  (  '  S  7  a  ' 

WR  I  T  E  <  6  ,  2  1  3  8  > 

FORMAT  (  '  R  5  a  ' 

WR  I  TE  •  6  ,  2  1  39  1 
F  ORMA  T  (  '  Z =  ' 

WR  1  T  E  «  6  .  4  9  9  > 


S6jT6.U6.V6 


C  2  O  7  /  ' 

T  6  a 

. G20  7/ 

4  U  6  *  ' 

.  G20  7  / 

V  6  a 

.  G  20 

7/  > 

R  4  ,  T  4  .  U  4 

.  V  4 

G  20  7  / 

T  4  a  ' 

.  G  20  7  / 

'  U  4  =  ' 

.  G20  7  /  4 

V  4  a 

.  C  20 

7/  > 

G  1  1  ,  H  1  1  . 

C  22  , 

M  2  2  ,  X  1  1 

.  X  22 

.  G20 . 7/  4 

M  1  1 

s  '  ,  G  20 

7  /  4  C  2  2  *  '  ,  C  2  O 

7  /  4 

H  2  2  a  ' 

G  2  O 

/  4  X  22=  ' 

.  G20 

7  /  > 

R  1  1  ,  R  2  2 

, G20  7/ 4 

R  2  2  =  '  .  G  20 

7  /  > 

G  3  3  ,  H  3  3 

.  C20  7 /  ' 

M  3  3 

=  4  .  G  20 

7  /  > 

S7.T7.U7.V7 
. G20  7/ 4  T  7  a  '  _  C 

R  5  .  T  5  .  US  ,  VS 
,  C20  7/'  T  5  a  '  .  C 

2 

G20 . 7/  > 

GF , SBL , C , SO , SF 


/  4 

U  7  =  ' 

.  G  2  O 

7/  4 

V  7  =  ' 

,  G20 

7/ 

/  ' 

U  5  * 

.  G  2  O 

7/  ' 

VS  a  * 

.  G  2  O 

7/ 
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227 

WRITE(6 

.  300 

) 

2  2  A 

300 

FORMAT < 

T15,  'X'  ,  1  OK , 

Y 

in » 

229 

4  9  9 

FORMAT ( 

'O', 

T  5  ,  'UNIT  STIFFNESS  (K1PS/IM) 

'  ,  C20 . 7 . 

230 

4/  T  5 

'  O 

1 d i mens  ions 

' , C20  7 

TS  .  '  C 

23  1 

A 

'  ,  C20  7 

/ 

T  5  ,  '  Snti  1  1 

d 

232 

A 

.  C20  7  /  T  5  .  '  Snu 

1 

233 

A 

'  , C20  7 

/ 

T  5  ,  ' 

23  4 

4 

C  2  0 

7  /  T  5  ,  '  rt  2 

,  C20 

235 

4.7  / 

T  5 

. 

*  .  G20 

7  ) 

2  3  6 

500 

S  T  OF 

237 

E  N  0 

End  of  f 

1  e 

240 


1 

REAL 

I  NC 

2 

CALL 

F  R  E  A  0 

S  . 

'  5R  *  4 

.  R  L 

3 

CALL 

FR  E  AD 

s  , 

'  6  R  *  4 

.  so 

4 

CALL 

FREAO 

5  . 

'  4  R  «  4 

.  A  1 

5 

CALL 

FR  E  AO 

S  . 

'  4R  •  4 

.  A  2 

WR  I  TE (  6  ,  200 ) 
FORMAT  <  TS  ,  ' 
N«R  L  1  / I NC ♦ 2 
K»* INC 


,  R  L  t  , R  L  2  ,  INC,  ALPHA,  SETA) 
,  S  0  I  ,  S  F  ,  S  F  1  ,  ML  ,  CF  ) 
.B1,C1,0U 
,  B2 .  C2 , 02  ) 

Y  1  *  ) 


1  0 

00  to  I ■  1  .  N 

1  t 

X  *  X  ♦  I  N  C 

1  2 

IF  (  ALPHA-BETA  >30, 35  ,  40 

1  3 

30 

Y* A  1  *EXP  |  (S01+SF1  )  •  X  1  *  B  J  •  EKP  (  »(S01*SF1 

>  *  X  )  *C  1  «  EXP  < (SOI 

"SF1  )  *  X 

t  4 

1  )  ♦  0  t  •  E  X  P  <  •  ( SO  t  *SF 1  )»S  )«HL/GF 

1  5 

CO  TO  50 

1  6 

35 

Y*A1*CXP(S0*X  >^B1*X*EXP(SD8X  )  ♦  C  t  8  E  X  P  (  - 

SO*  X  >  ♦  D  1  »  £  X  P  (  •  SO  *  X  >  <*HL  / 

1  7 

1  CF 

1  A 

CO  TO  50 

1  9 

40 

Y*EXP(SO*X)*  (A1*C0S(SF*X  >*B1*SIN(SF*X  ) 

)*EXP(  -SO*  X)e(Cl 

*COS(SF* 

20 

1X)401«SINISF*X)  > ♦ M  L / C  F 

2  1 

SO 

WR  I TE  <  A  .  1 00  )  X  .  Y 

22 

1  00 

FORMAT (  T5  .  2E  1  5  .  7  ) 

23 

1  0 

CONTINUE 

24 

WRITEIS,  300  > 

2  S 

300 

FORMAT  <  TS  .  '  X2  V2  ' 

> 

26 

N  *  R  L  2  /  I  N  C  <•  2 

27 

X  *  °  INC 

2  A 

00  20  I  *  t  ,  N 

29 

X  s  X ♦ I NC 

30 

IF ( ALFHA-BITA  )  (O ,  65  4  70 

3  t 

60 

Y*A2*EXP  (  <  SO  1 ♦SF  1  I  *  X  t ♦B2«IXP (  •  1  SO  HSF  1 

)»X)*C2*EXP( (SOI 

-  S  F  1  )  •  X 

32 

1 >  *028  EXP  < • (SOI -SF1 >»X>*HL/GP 

33 

CO  TO  AO 

34 

6  5 

Y*A2*EXP(S0*X  )*I2*X*|XR(S0*X  1  *C2*  E  XF (  ° 

SO*  X  > ♦02«  X«EXP  (  • 

SO*  X  > ♦ML 

3  5 

1  /CF 

36 

CO  TO  AO 

37 

70 

Y*EXP($0sX  >  «  <A2*C0S<SF«X  1 +B2 » S  I  N ( S F » X  ) 

(♦EXP  (  •  SO  »  X  1  *  (  C2 

*  COS  (  SF  » 

3  A 

1XU02>S1N(SF>X)  >  oHL /CF 

39 

AO 

WR  !  TE  (  6  ,  100  1  X.V 

4© 

20 

CONTI NUE 

4  1 

STOP 

42 

END 

End  ©  f  f  i 

i  1  « 

